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1. Introduction 


We start with a brief discussion of the physical motivation behind the 
mathematical considerations to be presented in Part I of this paper. 


Often in theoretical physics one is concerned with a functional relationship, 


o=8(9(5) (1.1) 


which states that the present value o of a physical quantity is determined by 
the values g(s) of a second quantity at all times s in the past. In (1.1), g(s) 
is the value of the second quantity s time units ago. If the functional § is given, 
the relation (1.1) may be used to predict the present value o of the first quantity 
from a knowledge of the “‘history” g of the second quantity. The relation (1.1) 
may be interpreted as expressing the causal nature of a class of physical processes. 
For definiteness we assume that the values of o and g(s) le in appropriate 
normed vector spaces, not necessarily finite-dimensional. 

In many physical situations, the value o in (1.1) will be, in some sense, 
more sensitive to the values of g for small s than for large s. Intuitively speaking, 
the “memory” of the system described by (1.1) will “fade away” in time. In 
order to make this idea precise, we shall introduce a norm ||g|| in the function 
space of the histories g for which (1.1) is meaningful. We first choose a number p, 
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1<p<ov, which will be kept fixed, and then define the norm of g to be 


pia ee 2 aa a 
= V/ F (gis) ms) Pas if 15 p<00 
| sup |g(s)| A(s) if p=oo, 


(1.2) 


lI | 


where |g(s)| is the magnitude of g(s) and where h(s) isa real-valued function 
which approaches zero rapidly as sco. Thus, in computing the norm lg || 
we assign a greater influence to the values of g for small s (recent past) than 
for very large s (distant past). We call the function / in (1.2) an influence 
function. The physical idea that the memory of the system is fading corresponds 
to the mathematical assumption that the functional § in (1.1) is continuous 
with respect to convergence in the norm (1.2) of the function space of histories. 
The influence function / characterizes the rapidity with which the memory is 
fading. 

It is always possible to reduce the relation (1.1) to one which is normalized 
in the sense that the possible histories g have the value g(0)=0 at the present 
instant and that the value of § for the zero history g(s)=0 is zero. Now, 
physical experience indicates that phenomena which one would expect to be 


described exactly by a normalized functional relation (1.1) often follow a simpler. 


relation of the form ‘ai dy d 


c= os 9a (1.3) 


s=0 


where { is a linear transformation. In particular, it appears that (1.3) accounts 
well for the observed phenomena in the case of slow processes. This observation 
leads to the conjecture that, in some mathematically precise sense, the relation 
(1.3) approximates the normalized relation (1.1) for slow processes. It is the 
purpose of the present paper to prove an approximation theorem of this kind. 
Theorem 2 of §5 asserts not only that (1.3) is the first-order approximation to 
' the normalized relation (1.1) for slow processes but shows also the form of the 
approximations of higher order. The theorem is based on the assumption that 
the functional § of (1.1) is not only continuous but also Fréchet differentiable 
at the zero history g (0) =90 in the function space with norm (1.2). A normalized 
functional § satisfying this assumption will be called a memory functional. The 
continuity of § for histories g different from the zero history is not needed to 
prove the approximation theorem. 

In Part II of the paper we apply the approximation theorem to constitutive 
equations of continuum mechanics. Our main interest there is in the logical 
status of the theory of Newtonian fluids within the framework of a recently 
proposed general theory of fluids with memory effects *. 


The theory of compressible Newtonian fluids is based on the following con- 
stitutive equation for the stress tensor S: 


S=—pI+2nD+A(trD) I; (1.4) 


here the rate of deformation tensor D is the symmetric part of the velocity 


gradient tensor; = (@) is the hydrostatic pressure the fluid would be supporting 
* Cf. [1], [2] and [3]. 
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if it were at rest at its present density 9; 7 and / are functions of @ alone and 
are called coefficients of viscosity. 

Experience shows that for some substances, such as water, the theory of 
Newtonian fluids accounts remarkably well for experimental measurements over 
a very wide range of conditions. Other substances, such as molten plastics, 
definitely do not obey (1.4) exactly, but yet their behavior appears to approxi- 
mate that of Newtonian fluids in the limit of slow motions. 

In §6 we review the mathematical definition of a simple fluid. As we have 
frequently remarked in the past, we believe this definition is capable of covering 
the behavior of nearly all real fluids including such substances as molten plastics 
which exhibit “‘hereditary’’, ““non-Newtonian”’ and ‘“‘viscoelastic” effects. Here 
we add to the definition of a simple fluid a new requirement: we require that 
the functionals occurring in the definition be memory functionals in the sense 
of the definition used in this paper. As we have indicated, this requirement 
is related to the physical assumption that simple fluids have a fading memory. If 
this assumption were not satisfied, the term “‘fluid’”’ would hardly be appropriate. 

In §7 we conclude that the theory of Newtonian fluids is indeed the complete 
first-order approximation to the theory of simple fluids for slow flows. We also 
indicate what an experimenter should expect to find as second-order corrections 
to the constitutive equation of a Newtonian fluid. We point out that several 
special flow problems for incompressible second-order fluids lead to third-order 
linear partial differential equations. 

In §8 we apply our approximation theorem to the theory of the general 
simple materials defined in [7], Part III. 


I. The Approximation Theorem 
2. Influence Functions and Histories 
An influence function h of order r>0 is a real-valued function of a real 
variable with the following properties: 
(i) h(s) is defined and continuous for 0Ss<o. 
(ii) h(s) is positive, h(s)>0. 
(iii) For each o>0, there is a constant M,, independent of «, such that 


A (s/o) Ze —< 
sup iG HOt 020 S1.. (2.1) 


The last condition (iii) means that /#(s) must decay to zero at a fast enough 
rate as s—>oo. In fact, we have 
sup s’ h(s)"==N-= oo, (2.2) 
s>0 


which follows from (2.1) by taking o=1 and observing that 


sS2o 


v / 
(2) ACE) SA) wp apy SAO Me 
Let a real Banach space Y with norm || be given. We then define a story 
g to be a measurable function defined for 0Ss<oo with values g(s) in %. Two 
such functions will be regarded as the same if they differ only on a set of 
measure zero. 
25* 
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For a given influence function # and a given number ~, 1<p<o~, we define 
the Y, ,-norm, ||g||,,p. or simply ||g||, of a history g by 


Ba ae ea ae 
lako= VF ligolayras it t3p<e, (23) 
gle =suplg (sas) if poe. (24) 


In (2.4) and subsequently sup f(s) stands for the essential supremum of f(s), 
i.e. for the greatest lower bound of the suprema of all functions which differ 
from /(s) only on a set of measure zero. 

The set of all histories with finite SF, ,-norm forms a Banach space, which 
we denote by % ». 


We remark that (2.2) is a necessary but not a sufficient condition for the decay 
relation (2.1). However, if a function f(s) satisfies (i), (ii), and the limit relation 


lim s’(s)=0O monotonically for large s, (2.5) 
S—>0O 


then it also satisfies the decay relation (2.1) and hence is an influence function of 
order v. To prove this fact we consider the expression 


h (s/o) _ (s/a)" h(s/a) 
GWE CIN ec Ce 


Since s’f(s) is monotonically decreasing for s larger than some value s,, it follows 
that (2.6) is not greater than 1 for all s>s,. For 0<oXsSs, the denominator s’h(s) 
of (2.6) has a positive minimum because / satisfies the conditions (i) and (ii). The 
numerator of (2.6) is bounded by the maximum of s’/(s), which exists and is finite 
because of (2.5). It follows that (2.6) is bounded, for oXs<oo, by a constant M, 
independent of «, which is the content of (2.1). 
The function 
1 


rir 


is an influence function of order y, An exponential, 


h(s)= 


PU ((S) =— Cit Ss Ole 


is an influence function of any order. 


3. Retardation 


The retardation I}, with retardation factor x,0<a<1 , 1s the linear transformation 
g —g,, defined, for all histories g, by 


(12.9) (8) = Gals) =g (as). (3.4) 
We show that J} maps the space ¥,» into itself. For p+00, we have* 
gall? =f (Ig (es)| A(s))Pds 


ae (s/ 1 Pall 
= af Mig ais? (Gas s lla + (sup Ay? 


h(s) 


* Since no confusion can arise, we omit the indices h, p on the norm, 
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Since /; is continuous and positive, 


o2s=0 h (s) 


for any o>0. The decay condition (2.1) implies that 


sup =e <o’M,<o. 


Hence, ace 
A(s/o) __ K 
DirG gO (3.3) 
Combining (3.2) and (3.3) gives 
IIgellSe ? Kallgll- (3.4) 


It is easily verified that this inequality remains valid for =co when we put 
‘Va=1. It follows from (3.4) that when g has a finite norm g, has a finite norm 
and is, therefore, in ¥%,,. Intuitively, retardation replaces a given history by 
one which is essentially the same, but slower. 

If, possibly after a suitable alteration of g(s) on a set of measure zero, the 
limits 


0 : h era) 
G=limg(s), #=lim 5 (gs) — > 9) (3.5) 
j=0 


exist for k=0,1,...,”, then we shall say that the history g has m generalized 
derivatives at s=0. Of course, the existence of ordinary derivatives implies the 
existence of the corresponding generalized derivatives. Here we shall use the 
term “‘derivative” always in the generalized sense of (3.5). Generalized derivatives 
at s=0 of any order may exist even though g(s) is not continuous near s=0. 

The Taylor transformation IT,, is the linear transformation g—J/T/,,g defined, 
for all histories g which are n times differentiable at s=0, by 


UIT, g) ( TT a gy (3.6) 


where the g are the derivatives (3.5). 

This Taylor transformation replaces g by its 2‘ “Taylor approximation” 
IT,,g. The history //,, g is a polynomial of degree <n. 

We note that retardation /} preserves the differentiability of a history and 
that //,, and [} commute; 1.e., for all histories g having m derivatives at s=0, 


n 


arey® II,.g=D.0,g- 3.7) 


na 


The set of all histories which have derivatives at s=0 and which also 
belong to the Banach space Y, , will be denoted by Y,. This set G, is a linear 
subspace of Y,,, but it is not closed in Y%,, and hence not a Banach space. 

The following theorem is an analogue of the classical Taylor approximation 
theorem. 

Theorem 1. Assume that n, p, and the order r of the influence function h satisfy 
the inequality 


atiesy: (5 =0 if p=). (3.8) 


360 BERNARD D. COLEMAN & WALTER NOLL: 


Then the Taylor transformation IT, maps the subspace Z, of L, , into itself, and, 
for all g in G,, 
lim = lta(g — 12, 9)|h,2 = 9- (3.9) 


We can also write (3.9), using (3.7), in the form 
Lng — 1, [9 = Ga— Un Ga = 9%"), (3.10) 


where the order symbol 0 must be understood in terms of the norm of the 
function space Y,,. Roughly speaking, the theorem states that in the space 
G,-L,,, a slow history is close to its Taylor approximation and that the 
distance between them is 0 (”). 


Proof of Theorem 1. We consider only the case when # is finite. The case 


when p=co can be treated analogously. 
First, we show that every polynomial of degree <” has a tinite norm. For 


this purpose it is sufficient to prove that 


J (sth(s))?ds=L,< 0 for kn. (311) 
0 
It follows from (2.2) that 


[e,2) 


f(t h(s))?ds = / (s" h(s))? sO ds < NP f s-P0-® ds. (3.12) 
1 1 Z 


The inequalities (3.8) and ku imply that #(y—k)>1 and hence that the 
integral 


is finite. Since s* h(s) is continuous, it then follows from (3.12) that the integral 
AA ris, finite: 

We have shown that any polynomial of degree < belongs to FL, , and 
hence to Y,, because it has m derivatives at s=0. Since a Taylor approximation 
IT,,g is a polynomial of degree <n, it follows that JT, maps J, into itself. 


The definitions (3.5) and (3.6) imply that the history 


satisfies the limit relation 
ee PS) 
tim == (ps (3.14) 


The definitions (2.3), (3.1), and (3.13) show that the assertion (3.9) of the theorem 
is equivalent to 
oh 
jim (|f(«s)| h(s))Pds=0. (3.15) 


We observe that 


Tt fre s)| h(s))? ds = ‘oe f (16 n(<\)'as (3.16) 
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and investigate the latter expression. Let e>0 be given. It follows from (3.14) 
that we can find a o(¢e)>0 such that 


ES for O0<sSXo(e). 


Hence, ole) Fe) 


ae ik (I7(o) n(\Fas <= pose ef (sn(2))Pas 


a"? 
ef ( = h( =\Pa() = Eos 
0 


where L,, is given by (3.11). On the other hand, we have 


ats [ (\to1A(2))fas—or f (IM 409))"(Gaey) a G18) 


a(e) a(e) 


IIA 


where w= f(r — nm) —1. Applying the decay condition (2.1) to (3.18), we find that 


gin f (toa Sass are, tae (3.49) 


o(e) 


On combining (3.17) and (3.19), we see that 


af (6 a(S) Pass PL, +o Mf, lll xii 


holds for all e>0 and all O0<a<1. The assumed inequality (3.8) insures that 
w=p(r—n) —1 is positive. Therefore, by choosing first ¢ and then « sufficiently 
small, we see that the right side of (3.20) can be made as small as desired. It 
then follows from (3.16) that the limit relation (3.15) holds, which completes 
the proof of the theorem. 


In the special case n=O, oe g=0 Theorem 1 states the following. If 


g ©,» is continuous at s=0 with value 0, then 
lim |Z g|| =o. (3.24) 
If ie and g—0, instead of (3.21), we have the stronger result 
9. =1,g =O (a). (3.22) 


4. Memory Functionals 
We first recall some definitions* from the theory of functions defined on 
a real vector space # with norm || || and having values in another real vector 
space 7 with norm | |. 
(1) Suppose g;,.--,g, are variables in #. Then a function $(gj, ..., YJ), 
derined for all values i the variables g; in # and having values in 7, is called 


* These demations are analogous to those given in Chap. XXVI of [4]. 
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a bounded k-linear form if it is linear in each variable g; separately and if there 
is a constant M, independent of g;, such that 


B(g,,--.,9x)| SM ||grll --- ll gall- (4.1) 


The form %(g,, -.-, gy) is said to be symmetric if any permutation of the variables 
leaves the value unchanged. 

(2) A function $(g), defined for all gc # and having values in 7, is called 
a bounded homogeneous polynomial of degree k if there is a bounded symmetric 
k-linear form $(g,, ..., g,) Such that 


Bg) =¥@,---.9)- (4.2) 
The symmetric n-linear form ¥(g,, -.-, Js) is uniquely determined by the homo- 


geneous polynomial %(g) and is called the polar form of the polynomial. 


(3) A function %, defined on a neighborhood of 0¢.# and having values in 
T is said to be n times Fréchet-differentiable at 0€ # if there are bounded homo- 
geneous polynomials 6°§(g) of degree k=0,1,..., such that 


3g) => BG) + Ig | RV), (4.3) 
R=0" 3 

WEL lim |9(g)| =0. (4.4) 
Ilg|l° 


The polynomial 6*§(g) or its polar form 6°§(g,,..., g,) is called the k' Fréchet- 
differential or the k variation of § at OCH. The differentials may be obtained 


recursively by 
Rl yj 
YFG) =80),  F(g) =A! lim, |B(Ag) — Y= 'G(G)|. (4.5) 
j=0-" 

We here consider the case in which # is the function space of all histories 
g with the following properties: 

(«) g has a finite Y, ,-norm, 

(B) g has m generalized derivatives at s=0, 

(vy) g has a zero limit at s=0: 

(0) 


hina. g'(s) == g-== QO. (4.6) 


s—>0 
(0) 2, p, and the order 7 of the influence function h obey the inequality (3.8). 
The conditions («) and (f) state that #< G,. Condition (6) insures that the 
conclusion (3.9) of Theorem 1 is valid for all g c%. 
A function § defined on a neighborhood in” of the zero function 0 CH CF, 
and having values in a real Banach space 7 will be called a memory jubchona 
of type (h,n) if it is m times Fréchet-differentiable at 0€ Y and if it is normalized by 


°F (g) =8(0) =0. (4.7) 


In some applications it may be more natural to assume that § is defined and 
Fréchet-differentiable on a neighborhood of zero in the entire space Y,,. How- 


ever, the approximation theorem of the following section applies only to histories 
which belong to the subspace ¥ of Y, ,. 
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5. The Approximation Theorem 


Theorem 2. Let § be a memory functional of type (h, n); the histories g in 
the domain space H of § have values in the space S, and § itself has values in J. 
Then, for each k-tuple of indices (j,, 75, ...,7,) such that 


$275 7.5°°354,5", Atiet---+1,=%, (5.1) 


there exists a bounded k-linear form (; 


with variables in SF and values in F 


such that, for all gc HX, it 
Big.) = D hy. DL vos Ga) -0(08); (5.2) 
where 
G9=1,g9, i€., ga(s)=g(«s), (5.3) 
and where 
pee (5.4) 


is the generalized 7 derivative at s=0 of g,, defined according to (3.5). The sum 
nm (5.2) ts to be extended over all sets (j,,...,7,) of indices satisfying (5.1), and 
the order symbol 0(«") has the usual sense. The multilinear forms {,_;, are uniquely 
determined by &. 

This theorem permits the asymptotic approximation of a memory functional, 
for “‘slow”’ histories, by a polynomial function of the derivatives at s=0 of the 
argument function of the functional. It is worth noting that the approximating 
expression 

=D halGe 9) (5.6) 

ae sii 
regarded as a functional of g, is not a memory functional in the sense of the 
previous section because it is not even continuous at the zero function g=0C*%. 


For »=1 and n=2, (5.6) reduces to 


(2) 


u(g)=49), lg) = L(g Veh (9. 9 9) +h(g ) 7) 


respectively, where {, and f, are linear and f,, is bilinear. 


Proof of Theorem 2. For 0, the theorem is a trivial consequence of 
(3.21). For n=1, the proof is based on a combination of the Fréchet-differen- 
tiability assumption (4.3) and equation (3.10) of Theorem 1 which, in the notation 
of (5.3), may be written as 


Gu =f, Ga + 9(a"). (5.8) 
Since g—0 by (4.6), the definition (3.6) gives 
a 9s) (s) — gee g! (5.9) 
The result (3.22) applies to both g¢# and ll, gc#H: 


Tl, $= 0(@), (5.10) 
J,=0(a). (5.11) 
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Combining (5.11) with (4.3), (4.4), and (4.7), we obtain 


(Gn) = 5 FB (Ga) + 0a"). (5.12) 
 *! 

Consider now a bounded homogeneous polynomial $(g) of degree k, 1 SkSn, 
in the sense of (2) of §4. The differentials 6*%, k>1 are such polynomials. 
The polar form of ¥ will also be denoted by ¥. Using the multilinearity of this 
polar form %, we obtain by (4.2), (5.8), and (5.10) 


= 


¥(g.) = 8 (HGa+0(a"),---, Ly, Ja + 0(a")) 
(I Oy, ce bie Oa) ) + DF (O( (ol <2), 5 OR, cos Us 


where each term of the sum contains at least one variable o(«”). It follows 
from the boundedness (4.1) of $% that the terms in the sum are all o(« aii 
Since k=1, we have 


¥B(Ga) = BL Ga) + 0 (2"). (5.13) 


We now investigate $ (//,, g,). Using (5.9) and the multilinearity of the polar 
form %, we obtain 
is ae ® : 
Sh ye ani sig, ..., iG). (5.14) 
ji=l k= 1J e: 

Due to the symmetry of $, all terms of (5.14) which differ only in the order of 
the indices 7,,...,7, are the same. Collecting these terms and separating all 

terms of order higher than 1 in «, we get 
BLL Ih OT OP Me pecs, (sir gd pat g) + O(a"*2) (5.15) 


(Jay «+63 7k) 
k fixed 


where the sum is to be extended over all k-tuples (j,, ...,7,) satisfying (5.1) and 


where the m,,_;, are positive rational numbers. 


Now, the function 
Vir. (Gis ++. &) = Bs q, -.., sa), (5.16) 


with variables a;¢ Y and values in 7 is clearly multilinear, because 93 is. Also, 
Y;,...;, 18 bounded. Indeed, application of (4.1) to (5.16) yields 


Wee (Cyt eae a,)| < M||sh al]... || ”* a, ||, 
and (3.11) shows 


[| s” a)|\; = | a,| {ae 


hence 
Bae Cay wake a,)|< M|a,|...|a,|. (5417) 
Substituting (5.16) into (5.15) and Be (5.4) yields 
BUT Fa) = Yin Gs <--> Ga) +O ("49 (5.18) 
+ Ik) 
ies 


where the f,,_;, are bounded k-linear forms. 


Functionals and Continuum Mechanics 365 


Applying the results (5.13) and (5.18) to ¥= 5, 8 in (5.12), we obtain 
formula (5.2). 

The uniqueness of the J, ;, follows easily from the observation that they 
are linearly independent and of order O («”) or lower. 


We remark that the multilinear forms [,__;, of (5.2) are not necessarily sym- 
metric. 


II. Applications 
6. The Concept of a Simple Fluid 
The notion of a simple fluid has been given a definition within the frame- 
work of a general theory of the mechanical behavior of materials*. This defini- 
tion is based on the following two physical assumptions **: 


(a) The present stress depends on the past history of the first spatial gradient 
of the displacement function. 


(b) A fluid has no preferred configurations. 


Using the principle of material objectivity ([7], §41), it was shown in re- 
ference [/] that the constitutive equation of a simple fluid can be written in 
the form (22.12) of [Z]. This functional relation, in a slightly different notation, 


reads 
le.) 


S(t) = (C,(t—s); @(é). (6.1) 


s=0 
Here S(t) is the stress and g(t) the density at time ¢. C,(t) is called the right 
Cauchy-Green tensor at time t relative to the configuration at time ¢. This 
tensor is defined by 
C,(t) = B'(t) F(x) (6.2) 
where 


F, (7), = Vos ys (22,:2) (6.3) 


is the gradient the displacement function §=y,(#, t) which gives the position 
at time t of the material point having the position x at time ¢. The stress tensor 
S(t) is symmetric. The Cauchy-Green tensor C;(t) is positive definite and 
symmetric; for t=t¢ it reduces to the unit tensor I: 


C)mT. (6.4) 
The functional § in (6.1) is isotropic; 7.e., § obeys 
Q 9 (C(9) 0) 07 9 (0 C(s) Q"; @) (6.5) 


identically in the history C(s)=C,(t—s) and the orthogonal tensor Q; here 
Q’ is the transpose of Q. It follows from (6.5) that the value of § for the “rest 


ANTS 24 
x 3 we give a survey of the theory of simple fluids with emphasis on physical 
applications. In that paper we anticipate some of the results rigorously derived here. 
Although the language and the definitions of the present paper are slightly different 
from those used in §§ 6 and 7 of [3], the arguments presented here prove also the 


theorems stated there. 
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history” €(s)=TI is a scalar multiple —p(o) of the unit tensor I. Defining 


B(G0).0) =P) F+9(F+ G),0). (6.6) 
we may rewrite (6.1) in the form 
S(t) =—p(o(t))1 +3 (Gl —s) —Te(t)). (6.7) 
The functional § (G(s); o) is defined for functions G(s) with the property 
me G(0) = 9, (0.8) 
and it has the value 0 for the zero function G(s) =0; 7.e., 
(0, 0) 0. (6.9) 


It is also isotropic in the sense of (6.5). 

We now assume that for each simple fluid defined by a constitutive equation 
(6.7) there exists an influence function / of an order 7 such that the functional 
% of (6.7) is a memory functional of type (h, m) in the sense of §4. The domain 
of $ is a class # of histories G whose values G(s) are in the space ¥ of all 
symmetric tensors. For the norm | | in Y we use 


| A| = \tr 42 for ACS. (6.10) 


a 


The range space 7 of § is the same as the range space of the histories G, 
7.e. the space Y =-¥ of all symmetric tensors with norm (6.10). Equation (6.9) 
insures that § has the normalization (4.7) required for a memory functional. 
The density @ enters into (6.7) only as a real parameter. The assumption that 
wis a memory functional implies that its domain of definition contains a neighbor- 
hood of zero in a function space # which is defined by the conditions («)—(6) 
of §4. The condition (y) expresses the assumption that all histories GCW 
correspond to motions which are continuous at the present instant s=0 (cf. (6.8)). 

If the simple fluid under consideration is incompressible, we must make some 
alterations in our starting assumptions. For every possible motion in such a 
fluid the density @ is constant and the tensor €,(t) is unimodular. In addition, 
the stress is determined by the history of the motion only up to a hydrostatic 
pressure p. Consequently, the equation (6.7) must be replaced by 


S@)=-p1+5 (C,(¢— s) — 1) (6.11) 


in which the indeterminate pressure p and the functional § may be normalized by 


a 


0=trS(t) +36 =tr§ (Ces) — 0), (6.12) 


If we were to limit the domain of the functional $ in (6.11) to kinematically 
possible histories, this domain would not contain a neighborhood of the zero 
function in an appropriate function space, %; therefore, § could not then be 
a memory functional. We assume, however, that § becomes a memory functional 
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when its domain is extended by putting 


co co 


& (C(s) — 1) =§ ((det CG) C(s)\— 1) (6.13) 


s=0 s=0 
when C(s) is not unimodular. 


Aside from the added properties (6.12) and (6.13), the functional $ occurring 
in (6.11) is assumed to be of the same type as that in (6.7) with respect to both 
isotropy and memory. 


7. Approximations of Order n for Simple Fluids 
We now apply Theorem 2, §5, to the memory functional occurring in the 
constitutive equations for simple fluids (6.7) or (6.11). 
Suppose a motion with Cauchy-Green tensor C,(t) is given. The correspond- 
ing history G is defined by 


G(s) =C,(t—s) —I. (7.1) 
If C,(t) is m times differentiable with respect to 7 at t=#, the k“ Rivlin-Ericksen 
tensor A,, kR=1, 2,..., mis defined as follows: 
dk - (k) 
Aa ae Cr(7) wes (=A)" G., (7.2) 
(k) 


where G is the k™ derivative of G(s) at s=0, as in §2. 

We now consider histories G, obtained from GC¥# by retardation as in (5.3). 
The corresponding Rivlin-Ericksen tensors are 

(k) 
Aj = a" A, = (— 1)" Gy. (7.3) 
(k) 

They differ only by the inessential factor (—1)* from the tensors G, to be sub- 
stituted for Je in the approximation formula (5.2). This formula, applied to 
the constitutive equations (6.7) or (6.11), yields the following expression for the 
stress tensor S, corresponding to the retarded history G,: 


Se=— PTLD thy (Aen AB) +010"), (7.4) 
Tyo Tk 
where the summation is to be extended over all sets of indices (7,, ...,7,) obey- 
ing (5.1). The terms m,,_;,(A%,..., 4%) are linear in each of the variables. 
For compressible fluids, it is understood that # and m,,__ ;,, depend on the density @. 
The equation (7.4) remains valid even when the derivatives shown in (7.2) 
exist only in the generalized sense of (3.5). 
The multilinear forms m,,__;, in (7.4) are isotropic functions, which means 
that they obey the identities 


Ont jo (Ui 55 U7) OF =m, ,,(QU, OF, 2%, 0.0;,,Q") (7.5) 


for all orthogonal Q and all symmetric tensors U,,..., U,. This proposition 
follows easily from the fact that the memory functional § occurring in (6.7) 
or (6.11) is isotropic in the sense of (6.5), from the observation that the conjugation 
G(s) +>QG(s)Q7 leaves the norm ||G|| unchanged, and from the uniqueness 
of the multilinear forms m,,__,;, asserted in Theorem 2, §5. 
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It follows from known theorems on isotropic functions that each form 


Wy... (Aj, ++», Ay), because it is isotropic and multilinear, may be expressed 
as a sum 
Mj, ip (Aj +1 Ap) = 2 Si (7.6) 
of products of the form 
Ci = Mi Pi Po = Pru (A;, Aj, --- A), = A,,, aii: Fa (7.7) 


where the g;'s are traces of products of some of the tensors A;, and are such 
that each A,, 7=1, ...,k, occurs precisely once in each term C;- In the case 
of an incompressible simple fluid the coefficients “; are constants, whereas for 
simple fluids in general the j; are functions of the density. Thus, for simple 
fluids, a finite number of scalar material functions (@), “4;(@) suffices to deter- 
mine the stress S, to within terms of order m in «. 

The case »=1 of (7.4) is of particular interest. With use of (7.6), (7.7), we 
obtain 


S.=— pI+ 1 At+ 4 (tr AQ I+ oa). (7.8) 


A 
a 2 
Now, the first Rivlin-Ericksen tensor 4,=2D differs from the rate of defor- 
mation tensor D only by the factor 2*. It follows that (7.8) is, to within terms 
of order o(«), simply the constitutive equation (1.4) of a Newtonian fluid. 
When the fluid is incompressible, the case n= 2 of (7.4) takes a remarkably 
simple form. The observation that tr A4,=0 for isochoric motions and use of 
(7.6), (7.7) yield 
S,=— pl +n At + B(Ai)? + y A$ + 0(07), (7.9) 


where 7, 6, and y are material constants and where # is an indeterminate 
pressure. This pressure # differs, in general, from the mean pressure p defined 
by the normalization (6.12), because it is obtained from p by absorbing all scalar 
multiples of J arising from m,, and m, through use of (7.6) and (7.7). 


Motivated by (7.9), we can define an incompressible second-order fluid by the 
constitutive equation 


S=—p1+74A,+BAji+y Ag. (7.10) 
Incompressible Newtonian fluids correspond to the special case B=y=0 of (7.10). 


In some dynamical situations equation (7.10) leads to a linear partial differ- 
ential equation for the velocity, just as in the Newtonian case. For example, 
consider a rectilinear shearing flow which, in Cartesian coordinates 5 My Soules 
a velocity field with components 


{v’} = {0, v(x, 2), O}. (7.11) 


If the body forces are conservative, substitution of (7.10) and (7.41) into 
CAUCHY’s dynamical equations leads to the following third-order partial differ- 
ential equation for v(x, t): 
GWT MONO 15. OM 
at CA On a One 
or O:7): 


ee pectin . (7.12) 
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A non-steady flow of the Couette type, in cylindrical coordinates 7, 3, z, 
has a velocity field with contravariant components 


{v} = {0, w(r, t), O}. (7.13) 
For such a flow, instead of (7.12), we get 


(7.14) 


ow 0 Ow roeD) 
7 = eke 
G 1 or ale 


In the Newtonian case, since b=0, equations (7.12) and (7.14) reduce to 
diffusion equations. As for the diffusion equations, many physically interesting 
solutions of (7.12) and (7.14) may be obtained by separation of variables*. It 
would be desirable to develop a mathematical theory of third-order partial 
differential equations of the type (7.12), (7.14). 


8. Simple Materials in General 


The considerations of the previous two sections may easily be extended to 
the general simple materials defined in Reference [7], Part III. 

The form of the constitutive equation of an tsotropic simple material may 
be obtained from (6.1) by replacing the scalar density o(¢) by the left Cauchy- 
Green tensor B(t), taken relative to an undistorted reference state (cf. [1] (22.10)). 
A consideration analogous to the one which led to (6.7) shows that the consti- 
tutive equation of an isotropic simple material may be written in the form 


S (2) = h (B(d)) +8 (Gilt—9 —I; Bit). (8.1) 


Here, the functional § depends on a tensor parameter B, instead of on a scalar 
parameter @ as in (6.7). We assume again that there is an influence function 
h of order » such that the functional § in (8.1) is a memory functional of type 
(h, n). In place of (7.4), we then obtain the following approximation formula 
for the stress S, se by a slow motion of an isotropic simple material: 


S,=h(B) + DW, ..;,(A4.-.-» 43; B) 0 (a*), (8.2) 
(1) «+> 7h) 
where m,,__;, is linear in each of its first # tensor variables but not necessarily 
in the last es: variable B. 

The formulae corresponding to (8.1) and (8.2) in the case of anisotropic 
simple materials are obtained from (8.1) and (8.2) simply by replacing all tensors 
T occurring in these formulae by their conjugates R7 TR with the rotation 
tensor R=R(t) of the displacement from the reference state (cf. [Z], (22. 8). 

In the case of isotropic materials, the function and the Heirs ein 
(8.1) are isotropic in the sense that they obey the identities 


Q(B) Q*=H(QBQ"), (8.3) 
8 (O G(s) Q7, 9BQ"’) = QB (G(s), B) Q* (8.4) 


* In particular, the sinusoidal vibration problems discussed for Newtonian fluids 
in §§ 345—346 of Lamp’s treatise [5] are easily solved for second-order fluids. Also, 
special solutions of (7.14) corresponding to sinusoidal vibrations of a fluid between 
coaxial cylinders can readily be found. 
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for all orthogonal tensors Q. As in §7, it follows that the My in in (8.2) are 
isotropic functions of all their variables; 7.e., they obey the identities 


O mW. ie (UL sey U,; B) o¢ = mM, 4.(Q U, Q*, meets Q U, Os OB Q*) (8.5) 


for all orthogonal Q and all symmetric tensors U,,...,U,, B for which m,, _ ;, 
is defined. The methods developed by SPENCER & Riviin ([6], [7] and (8]) 
may be used to derive explicit representations for the ™,, of a type analogous 
to but more complicated than (7.6), (7.7). Using such a representation, one 
can show that, in the case »=1, the approximation formula (8.2) reduces to 


S, = §(B) + Att, (B) + t,(B) Aj + tr (AZ f,(B)) f,(B), + 0(a), (8.6) 


where §, f,, f, and f, are isotropic functions of the one variable B and hence 
have representations of the form 


§(B) =6,1+ 6, B+ B, B’, (8.7) 


in which By, B; and f, are scalar functions of the three principal invariants of B. 


The first term §(B) in the expression (8.2) for the stress S, corresponds to 
purely elastic response. The sum in (8.2) may be interpreted as representing 
the internal friction for slow motions. 
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An Energy Theorem for Viscous Fluid Motions 


ROBERT FINN 


1. Introduction 


It is intuitively clear that an infinite amount of work is necessary in order 
to accelerate a body %, initially at rest and surrounded by an incompressible 
viscous fluid which occupies the (three-dimensional) region exterior to &, into 
a state of uniform motion with constant velocity, —w,. For if W(t) is the work 
done in accelerating the body from the initial time ¢ =0 till the time #, E(¢) the 
kinetic energy of the fluid particles at time ¢, and Q(t) the rate at which energy 
is converted into heat in the motion, one has always the relation (see, e.g., [J, 
p- 268]), 


(1) Wi ET Ot) at 


Letting u(x, ¢) denote the velocity field of the fluid particles at time ¢ and adopt- 
ing a coordinate frame fixed to Z, then the new velocity field, 26 (a, t) =w6,+ 
u(x, t) will tend to a steady motion w(a#) which vanishes on the boundary Y 
of Z but does not vanish identically. Hence we must have Q(t) > Q) > 0 for all 
sufficiently large ¢, and therefore 


W(00) = f Qodt = ov. 


This reasoning is in two senses incomplete. First, it is not known whether 
every (or even any) steady motion can be achieved by such a continuous defor- 
mation in which the velocity field satisfies at all times the hydrodynamical 
(Navier-Stokes) equations. The second, and presumably more lasting objection, 
is that no information is obtained about the relative magnitudes of the two 
terms on the right side of (1). 

We propose in this paper to answer the second question and circumvent the 
first by studying directly the kinetic energy content of a steady disturbance 
field defined exterior to Z. That is, let w(a) define a flow in the exterior & of 
B, and let w(x) +, at infinity, (a) =0 on &. We shall show that then the 
kinetic energy of disturbance, E = {|26(a) —26,|?4V, is necessarily infinite. 

é 


This theorem is a relatively easy consequence of results we have proved in [2] 

if suitable assumptions are made on the rate at which (a) tends to w,). The 

reader who is content with such a discussion, or who wishes only to see the 

motivating idea, can refer directly to Theorems 1 and 2 in $3. Our principal 

effort will be directed to proving the result under the simplest and most natural 

assumptions possible (Theorem 6). In this form, the theorem is applicable to 
Arch. Rational Mech. Anal., Vol. 6 26 
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the solutions whose existence is demonstrated in [3] and in [4], for which little 
information about asymptotic behavior is available, and shows that the kinetic 
energy of disturbance associated with such solutions is necessarily infinite. As 
additional reward for this effort, we obtain as corollaries of our method several 
results of independent interest, notably a general theorem on asymptotic behavior 
of solutions with finite kinetic energy of disturbance (Theorem 4), and a theorem 
of Liouville type (Theorem 7). 


2. General remarks 
We shall study time-independent solutions of the Navier-Stokes equations 


Aw — w-Vw—Vp=0, 


2) V-w=0. 


The notation is the usual one of vector analysis. The velocity field w (a) and 
pressure p(x), #@=(x,, X2, X3), are assumed to possess continuous derivatives of 
all orders entering in the equations. In (2) we have normalized the density @ 
and viscosity coefficient « to be unity. This can always be arranged by multi- . 
plication of w& and by suitable constants, and we assume this done. 

Corresponding to a given constant vector w#,, we introduce the linearized 
system 

Aw — w,:-Vw—Vp=0, 
(3) 
V-w=—O0O, 


which are the equations satisfied by a perturbation of a solution of (2) about 
the uniform flow w= w,. A fundamental solution tensor x(x, y) = (Xi j), P= (y) 
corresponding to (3) has been given explicitly by OSEEN [5, p. 34]. For each 


fixed 7 (y;;, p;) defines a solution of (3) in # and a solution of the adjoint system 
4 Aw +w,:-Vw—Vp—=0, 
V-w=0 


in y. Using this tensor, any solution of (3) in a region Y with boundary Y can 
be represented in the form 


(5) (a) = By + $ [H- Ty —Y- TH+ (y-W) (y-n)] dS + fy: VaAV 
z a 


where we have set w(a) = (ax) — wy, and terms of the form u-T® are to be 


understood in the sense w-Tv=u;(T@),,n;, with 
(TB) =— 95.4 (54 a sl, 6;,= 1,0=] 
0x; Ox; } 0, t=b7, 


Ve (n;) = unit exterior directed normal, and summation extended over repeated 
indices. A similar representation is available for the pressure, 


(6) P=Pot$lu-TH—b- Tu + (h- wu) (wWy-n)]dS + fp-w-Vuidv 
p3; ‘Sa 


where we have introduced a ‘‘pressure’”’ y*= 2, -V 


1 
——] corresponding to the 
vector (a, ¢). (rer) : 
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In the case that the region is a sphere V and w#,=—0, we shall need the Green’s 
tensor G(x, y), P(x, y) associated with (3). This tensor has also been given 
explicitly by OsEEN [5, p. 105]. Using it, we find 

=¢u-TGdS aA -VwdV, 
(7) i : 

=$w-TPdS+f[P-w-Vwav. 

B3 V 
The formulas of OSEEN yield the important estimates 
|x, y)|=max|xi]<Crey, [Vxl<Creg, 

(8) Ip(x,y)|<C roy, 

[Vp (x, w)|<CrZ3, [Ppl ee 77 


and for any fixed sphere, 


(9) (Gris; gi = O72 |[P(@,y|<Cre,, |VGa,y)|<Crz3 


ay’ vy’ 
uniformly in x and y, for some fixed constant C. 


Throughout this paper, the symbol C will be used to denote a constant, the 


value of which may change in the context. Thus, from the inequality « <C/1+ p? 
we may conclude « <C|f| as B varies over some range. Analogously, we denote 
by e« such a quantity which is known also to be small in magnitude. 


3. Preliminary Estimates 
Lemma 1. Let w(x), p(x) define a solution of (2) in the exterior of a smooth 
closed surface X, and suppose that for some constant vector wy, U(x) = w(x) — Wo 
satisfies the inequalities, 
fluPavV<co, l|u(e)|<M< o, 
é 


uniformly for all x in & Then u(a) +0 as #—->oo. 


Proof. Let « be an arbitrary positive quantity. Choose Ry, so large that in 
the region & exterior to a sphere about the origin of radius Ro, f |ul?dV<e. 


ERo 
Let # denote an arbitrary point such that |#| > R)+4, and let V, be a sphere 
centered at x and of radius a to be determined in the interval 1=a<2. 
We ai 
—/{G-u-VudV+ by Wy: Vudv + $u- TGdS, 


Va Xa 
+, being the bounding surface of v An integration by parts, using the fact 
that G(a, y) vanishes on 2, yields 
x) =fu-u-VGdV —fu-w)-VGdV+u-TGdS. 
Va Va xa 
From the estimates (9) and ScHwarz’ Inequality we find 


Ns Lay av,+ |i f fe pees cea 


F 
Va 


26* 
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> 


2qS<e, hence we may choose an a in the given 


By the choice of Ro, farsi 


range such that §japaSce. Using the second hypothesis, we obtain 
Xa 


(10) jwte) |< f Pav, +e. 
Similarly, for some 6 in the range 3Sb=4 we obtain 
lr@l<c f Plate. 
Vk om 


This latter inequality is valid uniformly for all x interior to V,. Multiplying by 
~2 and integrating over V, gives 


A F re 4 
* [uw (a)| |e (y)| 777 a 2 77 
PR alg O [es aVpamaeat f ley)! dV,\ +e<e, 
V, V, , 


by ScHwarz’ Inequality. Inserting this result in (10) completes the proof of 
the lemma. 


Lemma 2. Let (a), p(x) define a solution of (2) in a neighborhood & of infinity, 
and suppose W(x) >, as e>oo. Then V0 and p =o0(logr) as x>~. 


These results are proved in Theorems 1 and 2 of [2]. 
Lemma 3. Under the hypotheses of Lemma 2, suppose in addition that { |u|? dV 
& 


<o, U=W—w,. Then there is a constant py such that p(a)>p) as >. 
If X' is any smooth closed surface which together with its exterior lies in &, then 


(14) (a) Sees Tb —- TU+ (W- &) (w,-n)] a> te u-Vudv, 


& denoting now the exterior of &. 


Proof. By ScHwaArz’ Inequality, the estimates (8), and Lemma 2, we see 
that the volume integral in (11) is absolutely convergent. Further, this term 
tends to zero as #->oo. To see this, let “2 be the annular region bounded by Y 
and by a sphere 2, about the origin of (large) radius R, and let & be the exterior 
of Lp. Let V be a unit sphere centered at x. Then for || >R+4, ve agah f+f. 


V @r-V 
For given e > 0, we may choose R so large that ju |?dV<e. ice ee of the 


[o-u-VudV\<e 


ae uViiav|< é. For this choice of R and x, we have 


estimates (8) and Lemma 2, we may choose ay 50 large that 


and also 


fv U- rar =e | iileay f parce 


Ep—V Er 


and hence | fp-u-VudV|<e, the desired result. 
E 
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In the (finite) region “,, we have by (6) 
p(x) =$[e-Tp—P- Tu + (P-t) (wy-n)JdS+$[ JdS+fp-u-Vudv. 
g ZR Pp 


Fixing # and letting R->ov, we see that of the two terms which depend on R, 
the latter tends to a limit, hence so must the former. Set 


F(x)=lim $[u-Tp—w- T+ (b- a) (u-n)] dS. 


Roo ie 
In this integral, all terms clearly tend to zero with the possible exception of the 
middle one, which involves the pressure. Let 
H(x)=lim fb-Tuds, 


ROR 


and denote by 6”H(«) a difference quotient corresponding to a displacement h 
in an arbitrary direction. A simple application of the mean value theorem shows 
that 

6" H(a) = lim 6d'b- TudS=0 


in #, by (8) and Lemma 2. Since / is arbitrary, we conclude that H(x) = F(a) = 
const =/,, and Lemma 3 follows directly. 

Theorem 1. Let 26 (a) be a solution of (2) exterior to a smooth closed surface X, 
and suppose that for some constant vector Wy == 0 and some «> 0, | 6 (ax) — Wo| < 
Cr-2-* as £00. 

Il, pes => >); > > > 
: EA EID STATE 10S NEST 


Then if U(x) = ww (x) — Wy, we shall have f |u(a)|2dV <co if and only if a=0 
Proof. By Theorem 8 of [2], 
U (a) = Uy(a) + Ue (x), 
with 
Uy (x) = ${-Ty—y- TUF (X- Ww) (Wn) + (X- w) (Wen) dS, 


and U,(a), as is easily verified from the estimate (27) of [2], has finite square 
integral over &. A simple application of the mean value theorem, using the 
estimates (8), 


iy (x) = — a- x(x, 0) + tt, (x), 


where Uj; (a) is square integrable over &. (The point y=0O is assumed interior 
to »'.) Thus, it suffices to show that a-y(#, 0) has infinite square integral 
whenever a+ 0. Set @=(d,, a, a3). Then |a- x (x, 0)|? =; ; Xin 4 &, a quadratic 
form in the components of a. We integrate this fon termwise over a sphere 2p 
about the origin, and observe that by the symmetry properties of 7;; (see, ¢.g., 
the defining relations [5, p. 34] of OsEEN) all terms for which 7 + # integrate to 
zero. Thus, we need only study the sum of squares 47, a7 for each of the indices 


* The vector a has a physical significance as the sum of the force exerted on 2 
by the fluid, the flux of momentum across 2, and a momentum flux arising from 
the motion at infinity. 
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i =1, 2,3 (no summation on 7). Tedious but formal estimation yields 
$yiidS>CR4, $y3.dS>OR4 6y33dS >CR4, 
ZR ZR =R 
for some positive C, and from this the theorem follows. 
Theorem 2. Let ae be a solution of (2) exterior to a smooth closed surface &, 
and suppose that |W (a) —w |< Cr-4~* as x00, Wy + 0. If w(x) =0 for x on 2, 
then necessarily aA: tee ibe 


Proof. If q is the force on 2, then by Theorem 13 of [2], 


4 [ Ou; Ow; ray > Of 
ass - 0 71 . 
Z ( OX; se OX; } . 


It follows that 4+ 0, and the result is obtained directly from Theorem 1. 


Lemma 4. Let y (a) be a vector valued function defined and piecewise continuously 
differentiable* throughout space. Then there exists a constant vector ¥ such that, 
for any choice of the origin of coordinates, 


(essa iV 4f \vypar, 
é 


the integration being taken over the entire space €. 
This lemma is proved as Lemma 3.4 in [6]. 


4. The Main Theorems 


Theorem 3. Let w(x) be a solution of (2) in the exterior € of a smooth closed 
surface X, and suppose that for some constant vector Wy, U(x) = W(x) — Wy satis- 
fies the inequalities 

f\aPdV<ow, |e(x)| <M<oo, 
6 


uniformly for all xin &. Then w(a) has finite Dirichlet Integral over &, [ \Vw|?dV 
é 


= 00, Denoting the value of this integral by A and the contribution from the annular 
region SF, bounded by X and by a sphere Xp of (large) radius R by D(R), we have 


As Ps = “ap (u-n) — Jai)? (wy-m) — p(i-n)} as 


and 


|A — D(R)| < CR. 


Proof. Lemmas 1, 2, and 3 show that 2 (a) >2,, V26 (x) +0, and p(x) +f, 
as #->co. We shall assume, without loss of generality, that fo=0. Formal 
integration by parts yields the relation 


(12) —A= Once oH aap (ue -) — |e? (wy -n) — piti-n)bas. 


ZR 


x 
That is, continuous throughout space and possessing bounded derivatives which 
are continuous except perhaps on a finite number of smooth surfaces. 
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set Hir ae jaj2dS, K(r) = =H vy) dy, and let us denote an interval of the form 
i, ne by Te, There must bist an 7) such that for every integer 7> 7, at least 
half of the vale Kr i) 1,=V,¥—1,..., 27 —1 which occur in any J, satisfy 
the inequality K(7,) <—. For otherwise, there would be an infinity of disjoint 
intervals {J,}, on a 453 which at least half the es K(7;) would 
satisfy K(r ee. ar . We would conclude fH ") dr >= “- for each such 


I,, and hence ffipav = oo, a contradiction. 7 
é 


Denote by {£;} the set of such disjoint unit intervals E, =[7,,7;,1] for which 
K(7;) cs On Si such interval there must be a set S; of measure wv (S;) — 
on weet (rv) <—. It follows that there is a fixed 6) > 0 such that on each E; 
there are points a v,, and a set S? of measure = 6, contained between these 
points and of distance = 06, from each of them such that on So. H<—. 

Consider now $ p(u-n) dS for vin some S?. Since 0, we have by some? 
Inequality, oe 


(25 DIES OSA 


€ 
<= Ollas. 


We proceed to estimate the integral on the right, using the representation (11) 


for p i: : 


(13) 


ao Typ —p-Tu— (p-u) (@-n)]dS+ fp u-Vaidv 


é 


i -u-VudV + o(-.| : 


In the interval E;, we consider the subinterval |[7,, 7,|, and deform the volume 
integral by an integration by parts: 


(14) as i-VudV + ae w- u) (w-n)dS— [w-é-VpdV +0(5), 


Xr +2, 6—Ay 


where .o% is the annular region bounded by 2), and 2,,. We evaluate the integral 
of |p| over 2, by integrating (14) under the sign. We must therefore study 


$ |p (a, y)|aS, and if [Vep|dS,.. 


We distinguish cases (a), (b), (c), according as y lies in 7, on 2,, or &,,, or in 
€ — A. 

Case (a). Let 9 be the distance from y to 2,, @ the distance from y to &. 
Using (8), we obtain 


G 
|p (x, y)| < 02 +(v+ 5) —27 (7+ eo) COS@ ’ 
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being the polar angle measured from the origin, with the line joining the origin 
to y as axis. Thus, 
hip. y)|4S2< Clog ” 
0 


when 7 is large compared to 99. The contribution to $ |p| dS arising from ¥% 


y 
Saat 


[w-vaavpyas|se]/ [toe atel, [ |ujpav 
. 0 . 
ar Sy 


Ly Ay 


is thus 


= erlogr- |/4 =e |rlogr 
by the choice of the interval E;. 
Case (b). We have 


oe (u-n) ya S| = Coe 5 ple Pas < | t08 5 
0 
with a similar estimate for 2’. 


Case (c). Since |Vep| < C/r 


|7.y» we find for y in oes SAPD OS 


Caer u- ei we 


— A, 


Collecting these results and inserting them in (13) yields 


TAS 


€ , 
y 


Sls jy 087 


— 


for 7 in some Sy, Hence we obtain from (12) that for all such 7, 


|D(r 


and therefore 


1 e Cf adi) 
(15) | D7) SAP ae ae tee 


whenever 7 is in an SP. 


Suppose now that for some R,, D(R,) >A. Then since D (7) is non-decreasing, 
we obtain from (15) 


(16) tO Wap) Rea re S? 
DG) A | dee 0 elsewhere. 


Now in every interval I,=[r, 27], at least half the unit subintervals, 7.e., 
the intervals E;, contain a subset S? of measure at least 59, do independent 
of vy. Integration of (16) from R, to R thus yields 


1 
DiRyoa ~ Dara = © 4o(R? — Ri). 
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This is utter nonsense, as we see by letting Roo, so it must be that Di <A; 
all y, and therefore 


1 dD(y) r(1—e) if réS? 
A= DY) s-d7 0 elsewhere. 


Integration now yields 


¥\ 3 = 050,(h eR?) 


from which 
A 5D) e) CR 


for all sufficiently large R. This completes the proof of Theorem 3. 


Theorem 4. Under the hypotheses of Theorem 3, if Wy +0, then | w (a) —w,| 
<Cr-3+* as w%->oo, for any e>0. 


Proof. For x in &, the representation (5) gives 
u(x) = — $[u-Ty—y%- Tu + (x -U) (Wy-n)]dS— fy uVuad. 


(A discussion analogous to the one at the end of the proof of Lemma 3 shows that 
the other boundary integral vanishes in the limit.) The estimates (8) show that 
the surface integral over 2 has the behavior O(r) at infinity. To study 
the volume integral, we let |#]=2R and introduce the region Yp, Bee 


by 2'and by a sphere about the origin of radius R. For ally in %, |x (a, y)| <= 
by (8). Hence, by hypothesis and by Theorem 3, 

fy u-Vu av| iO Rh V4 i apav |) f / f \Vulzav 

Sr 


2CR?: 
Finally, setting &=& — Zp, 


Hix a. benef av|/ [\Pupar 
ER ER 


AGI 


by Lemma 4 and Theorem 3. Thus |u(a)|<Cr?. The method of proof of 
Theorem 1 now shows that the vector @ must vanish, and the stated result 
then follows from Theorem 8 of [2] *. 

From Theorems 1, 2, 3, and 4 we obtain the following results: 

Theorem 5. Under the hypotheses of Theorem 3, if there 1s no flux of fluid 
or momentum flux across 2’, then the force exerted on & by the fluid ts zero. 

If 7%, 0, Theorem 5 follows from Theorems 1 and 4. But if W,=0, we have 
the following lemma: 


* One sees easily that if in addition there is no net flux of fluid across 2, then 
u(x) satisfies the stronger estimate (27) of [2]. For the application of Theorem 4 
to the later results of this paper, the estimate |u| < Cyr! already suffices. 
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Lemma 5. Under the hypotheses of Theorem 3 af wy=0, then 
(17) w(x) = 6 {y-TwH—W- Ty —(X-#) («6 -vi)}dS + O(). | 
Proof. As in is demonstration of Theorem 4, we have 
(a) = am [(wo-Ty—-y-Tw]ds getty Se Be Bo 
An integration by parts changes this expressions to 
Ww (x) =—$[w6-Ty—x- TW +(x) RUB TONNE ot) AA te 


Again we evaluate the volume integral separately in the two regions -Y, and 
&. Since* for ®)=0, |Vx| < By we have 
Voy 

[6-w. VxaVv 
Er 

One sees easily that a6 (a), considered as defined only in é, can be extended 
to the entire space (say, by invariance at inverse points with respect to 2p) 
to have Dirichlet Integral not exceeding twice the integral over &. Thus, 


oa hey 

C| Pek gy 
y2 yu 
xy 


ER 


fBias C [ |VePav < CR 
Voy 
ER Er 


by Lemma 4 and Theorem 3. In -%g we have uniformly |Vx| Pee hence 


AdV eC aR S 


fw-w- Vxav| =< OR Alen 
Pp é 


and the lemma follows. 


Using Lemma 5, we may complete the proof of Theorem 5. For (17) shows 
that w(x) =${}dS+W,(x), where w(x) has finite square integral over &. 
2 


This is the decomposition needed for the proof of Theorem 1, which now shows 
that under the assumptions, F =0, g.e.d. 


Theorem 6. Let w(x) be a bounded solution of (2) defined exterior to a smooth 
closed surface &, and suppose W(x) =0 on X. Then for every choice of by, we 
shall have f |v — wW,|?dV =o, or else w(x) =0 in &. 

& 


Proof. Suppose that for some WW, f|%—w,/?dV<co. If w&)=0, then 
4 


v6 (x) =ti(a), and hence the constant A of Theorem 3 is zero. Hence, by Theo- 
rem 3, f|Vw|?dV =0, and therefore w(x) =0 in &. If %&)--0, we conclude 
i 


first from Theorem 4 that |26 (ax) —W,|<Cr+ and then from Theorem 2 that 
f |w — wy dV =o, a contradiction. 
oe 


Liens: W, + 0, the best uniform estimate available is |V'z|<Cyv~3, which is not square 
integrable over &. This is the reason the elaborate apparatus of [2] is necessary 
in the general case. 
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Finally, we may obtain a theorem of Liouville type: 


Theorem 7. Let wW(a) be a solution of (2) which is defined throughout three- 
dimensional Euclidean space &, and which satisfies the conditions, 


i) | (x)|< M <0 throughout € , 
ii) J | w(x) — 6, )?dV < 
é 


for some constant vector W,. Then w(x) = Wb, in €. 


Proof. By Theorem 3, f |Vw|?dV=0, hence (a) =const. By ii), 6 (a) = 
Wo, 9.e.4. é 


Note added in proof: After submitting this paper I learned that a similar theorem 
has been proved for motions in two dimensions by W. WoLIBNER [Studia Math. 
12, 279—285 (1951)], who obtained it as a corollary of a method for determining the 
force exerted on 2 in the motion. Important extensions of WoLIBNER’s theorems to 
time-dependent motions of compressible fluids have been given by A. KrzyWIcxI 
[Studia Math. 15, 114—122, 174—181, 252—266; 16, 48—55 (1955—1957)]. 
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Self-Superposable Magnetohydrodynamic Motions 


RICHARD R. GOLD 


Communicated by C, TRUESDELL 


1. Introduction 


The formal concepts of superposability and self-superposability of nonlinear 
partial differential equations, in particular those equations of interest in hydro- 
dynamics, have been examined in [/] and [2]. Following a brief review of the 
literature and some remarks regarding the philosophy of the subject, super- 
posability and self-superposability conditions were established for the equations 
governing motion of viscous fluids, compressible and incompressible, Newtonian 
and non-Newtonian. Corresponding conditions were determined for the equations 
of motion of magnetohydrodynamics and for the potential equation of gas 
dynamics. 

The present paper is an extension of the basic ideas set forth in Chapter 8 
of [2]. Following these formal results, we now seek solutions of the resultant 
system of equations, 7.e., we wish to construct self-superposable magnetohydrc- 
dynamic motions. The appropriate equations are given in the next section. 
Solutions are obtained in Sections 3 and 4 for two and three-dimensional motions, 
respectively. Both the cases of steady and nonsteady flow are considered. 

The mathematical complexity arises, in part, from the nonlinearity of the 
basic equations resulting from the terms (U- V) U and (B- V) Bin the momen- 
tum equation and curl (Ux B) in the induction equation. As a result of this 
difficulty a large number of the magnetohydrodynamic problems which have 
been treated with regard to such questions as stability, shock waves, magnetic 
fields, and aerodynamic applications involved linearizing techniques of one sort 
or another. In the astronomical literature one finds that very often U is either 
identically zero or is assumed to be known, in which cases the equations are 
linear in the remaining electromagnetic variables. Several authors have begun 
the study of the full equations under simplifying assumptions which make the 
system more amenable to solution, as is noted in the discussion of Case 2 in Sec- 
tion 4. It is at this point that the introduction of the concept of superposability 
may be of definite use in exploring the character of the basic equations by pro- 
viding additional exact solutions and general properties. In general this is our 
primary purpose, and solution of specific boundary-value problems is not at- 
tempted in this work. 

Obtaining exact solutions of a particular system of nonlinear differential 
equations is but one of the potential applications of the concept of superposability. 
The pertinent conditions would allow one to examine general classes of functions 
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which may have this property of additivity (see [1] and [83], for example). In 
view of the role that the additivity or superposability property plays in almost 
all methods for solving linear differential equations, an examination of these 
methods applied to particular nonlinear equations in the light of generalized 
superposability results seems warranted. Finally, it is worthwhile to study the 
several formal mathematical questions which may be raised with regard to a 
subclass of solutions of a nonlinear equation or system of equations (or, more 
generally, to nonlinear differential operators of interest) which exhibit this 
linearity. 
2. Magnetohydrodynamic Equations. 
Superposability and Self-Superposability Conditions 
The motion of an electrically conducting fluid in the presence of a magnetic 
field obeys the equations of magnetohydrodynamics. The fluid is treated as a 
continuum, and the classical results of fluid dynamics and electrodynamics are 
combined to express the phenomenon. For a viscous incompressible fluid with 
constant properties (7.é., 6, l., €, @, v) the equations of momentum and of in- 
duction constitute the full magnetohydrodynamic system: 


eu/at+ (U-V) U—4(B-V) B=—Vs+rP?2U, 


(2.1 a) 
A=(420u,)1=constant, s=p/o+2+1B?/2, 


6B/ot =V x (Ux B) +V2B, 


(2.1b 
” = (4a0u,) 1 =constant, 


where the external force on the fluid other than the Lorentz force is expressed 
by the potential function 2, U is the velocity vector, # the pressure, @ the den- 
sity of the fluid, uw, the magnetic permeability, ¢« the dielectric constant, » the 
kinematic viscosity, o the electrical conductivity, and B the magnetic induction 
(the magnetic field strength or intensity is H=Buw;"). Inherent in these equations 
are the following assumptions, often referred to as the magnetohydrodynamic 
approximations. The electrical component, e, FE, of the ponderomotive force 
in the momentum equation is small compared to the magnetic component, 
j x B. The displacement current in Ampére’s law and the convection current, 
o, U, in Ohm’s law are negligible compared to the conduction current, o E. 
A detailed discussion of these equations and assumptions is available in a number 
of references ([4]—[6]). 
Both U and B are solenoidal vectors; hence 


V-U=0, V-B=0. (2.1¢, d) 


System (2.1) must be solved for U, B, and ~. The electric current density, j, 
is then readily obtained from Ampére’s law, j =(4% u,)+V xB. The electric 
field strength, E, is obtained from Ohm’s law, E=jo!—UxB=7Vx B— 
Ux B. Clearly V-j =0, while the divergence relation for EF, V-E=4z 7" @,, 
determines the excess charge density @,, which need not be zero in general. 
The superposability conditions can now be derived in the manner outlined 
in detail in [7] and [2]. Briefly, if (U,, B,, s,) and (U,, By, s,) are two solutions 
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of (2.1), what is desired is to determine s,+ 31+ SO that (G+ U,, B, +B, 
s,-+8.+2) is also a solution. Introducing the three solutions into (2.1a) and 
subtracting the latter result from the sum of the first two, one obtains 


_ (U,-V) U, — (U,V) U, + AL (By: V) By + (Ba-V) By] =Va.  (2.2a) 


Similarly, from equation (2.1b), we have 


—Vx(U, x B,) —V x (U, xX B,) = 0. (2.2b) 

For self-superposable flow, U,= U,= U, B,= B,=B and equations (2.2) become 
—(U-V)U+A(B-V) B= 3Axz, (2.3.a) 

VY x(Ux B)=0; of ©U xX B=Vyz. (2.3 b) 


In order for the system to be self-superposable, 7.e., if (U, B) is a solution, (” U, 
nB) is a solution, there must exist a pair of single-valued scalar functions z 
and y defined by (2.3). The conditions of integrability for (2.3) are 


curl [(U-V) U] =Acurl [(B-V) B], (2.4a) 

cunl (CU << Bye 0. (2.4b) 

On the assumption that the flow is self-superposable, system (2.1) may be 
simplified by replacing the nonlinear quantities with a and y via equations (2.3). 


Eliminating z from the resultant momentum equation (2.1a) by taking the curl 
and noting (2.1b), one obtains 


Obict==p V7 baw nro aur, (2.5 a) 
éB/at =nV?B. (2.5 b) 


Equations (2.1c, d), (2.4) and (2.5) comprise the self-superposable system which 
determines U and B. The pressure would then be determined from (2.1a). For 
steady flow (2.5) become 


Vee =O (y= 0); (2.6a) 
ViB=0 (40). (2.6b) 


Similar results were obtained in [2] as well as in [7]. For B=0 the conditions 
and equations noted above reduce to the nonconducting case discussed in [1]. 
It is of interest to note that for steady flow equation (2.5a) is satisfied for a 
perfect fluid (» =0) and § need not be harmonic; thus (2.6a) no longer applies. 
Similarly, for a steady flow in which the electrical conductivity is infinite (7 =0) 
equation (2.5b) is satisfied, and B need not be harmonic. 


It may be remarked that the self-superposability condition for equation 
(2.1a) expressed as equation (8.2.2) of [2] (equations 2.3a or 2.4a above) can 
be written in the alternate form 


curl (U x curl U) = Acurl (B x curl B) (2:7) 
by using equations (2.1c, d) and the vector identity 


(F-V)F=17P— Fx cul BF. 
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Force-free fields which are of considerable interest in the study of static equi- 
librium (U=0) in astronomical problems ([8] and [9]) are characterized by the 
condition Bx curl B=0. For U=0 (2.7) reduces to Bx curl B=Vw, which 
contains this result. The self-superposability relation is in fact a generalization 
of this force-free field condition to the case where U0. Since B satisfies the 
diffusion equation (2.5b); in the present system there is an obvious analogy 
with the several electromagnetic and astronomical analyses of the magneto- 
hydrostatic problem, where putting U=0O in the induction equation (2.1b) 
results in this same diffusion equation ([4, pp. 4 and 5]; [3] and [14]). 


3. Plane Flow (U independent of z, and U, = 0) 

Let us assume in this section that the magnetic field is restricted to the plane 
of motion, 7.e., Bis independent of z, and B,=0. Introducing the stream function 
w(x, y, t) along with a comparable function ®(x, y, ¢) for the oe field so 
as to satisfy (2.1c,d) identically (YU=y,, U=—y,; B,=@,, Bk 
where the comma in the subscript denotes partial Higaeanatony a cyeiae 
equations (2.4, 5) to the form 


(Ay) =oVy, (3.4) 

O(p, Vi p)/o(x, y) = 10(@, Hi ®)/A(x, y), (3.1) 
D ,=Vi ® + fi(d), (3.4) 

O(p, P)/A(x, y) =0. (3.4) 


The vorticity §, which is normal to the plane of motion (§,=&,=0), is €,= 
= xx—V,yy=— Vip, where Vi =V7 (V7). The Jacobian of the functions h(x, y,/) 
and g(x, y, t) is given by O(h, g)/0(x, vy) =h ,&y—h, g,~, and f(t) is an arbitrary 
function of time. Equations (3.1) must be solved for y and ®, while (2.14) 
reduces to the following equations for the determination of the pressure, p/@ +2: 


(plo+ 2) -=—Vyi— Fly te yl ty, Wy—18, “WG+r(MWy],,, (3-2a) 
(pio+ 2) y=¥, 21 3 ys $y bea CP y* Vip—Aa Dy. Vi O—y[Viy] > (3.2b) 


(ble+ 2) .=9. (3.2¢) 
In plane polar Se ihe vr, %, t) and ®(r, #, t) are defined by the conditions 
U=r'y 5, U=—y,,; B,=11® 5, Bsx=—@,. Equations (3.1) have the same 


form with the operator V2 now defined by Vi=@/072+7r10/0r +r 72/00, and 
the Jacobian relations now given with respect to y and @. Equations (3.2) apply 
using this definition for V7 and replacing 0/0 by 6/07, o/dy by 7190/00. 

In order to solve system (3.1) for y and @, let us introduce general orthogonal 
curvilinear coordinates, « and f, in which case by setting «=y the Jacobian 
relations (3.1b, d) may be simply satisfied. The functions h,(«, 8) and hg (x, B) 
are the usual metric coefficients, (ds)?=hj(da)?+A3(dB)?. Equation (3.1d) 
implies immediately, since y,,=1 and py ,=0, that 


D —— 0) . ies a) 
The vorticity is given by 
E,= — Vip = — (hy he)? (hg ht’) a (3.3 b) 
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while, by use of (3.3a), ¢,=—Vi® becomes 


C= co (Ay Np) > (hg hi” ® .) .= > hy” Daa at os De. (3.3¢) 
The Jacobian relation (3.1b) reduces to E, g=A®_.E,,, or, by use of (3.3, ¢), 
(1 —A@®*,) E, p= —A(h’) .p P ~Digp. (3.3 d) 


Using (3.3a—c) in (3.1c) reduces it to the form 
@ , = (hy h,)+ (hy hy* ae aia f(t) = UT] hy? OR Uf] ? ® , + f(t) ‘ (3.3) 


Applying 0/08 to (3.3) and using (3.3a) gives (hi?) p® aa=&s,pP, 4, and substitut- 
ing this result into (3.3d) yields 


E,,=0, (3.3 f) 
from which it follows that 
(hi?) ® ..= 0. (3.4a) 
By use of (3.3 f), (3.1a) becomes 
ae =e (Ay hy) = (hy hy* ee aa, hy” orae ag 67% ‘ (3.5 g) 
Applying 0/08 to this result and using (3.3 f) yields 
(11”) €, aa 0. (3.4b) 


In view of (3.3a) equation (3.1c) would immediately require that ¢,,=0, in 
which case (3.3f) follows directly from (3.1b). In addition to being useful in a 
subsequent application, the somewhat lengthier discussion given above emphasizes 
the fact that €, ,=0 is required for all solutions ®, including, for example, the 
case where the nonlinear terms are identically zero (w=0), 7.e., ®B=A-*y. The 
significance of (3.3f), namely, that €,=&,(«, ¢), is quite apparent. First it follows 
that (3.1b) can be satisfied only if both sides vanish identically, so that 0(y,V7y)/ 
d(x, y) =0. (For steady flow V?®=constant, and (3.1b) reduces immediately 
to this result.) This implies that the vorticity (in addition to the magnetic field 
in view of (3.1 d)) is constant along the streamlines, where y =constant. Finally, 
equation (3.1a) along with (yw, Vi w)/@(x, vy) =0 defines p(x, y, ¢) in the corre- 
sponding nonconducting case, obtained by KAmP£ DE FErrET (see the discussion 
in [1]). To this is now added the related magnetic field for a conducting fluid 
and the correction to the nonconducting pressure distribution. 

All the solutions y, ® of (3.1) may now be obtained from equations (3.3), 
(3.4). The more interesting ones are discussed in the following outline. 

(a) For €,,=-0 equation (3.4b) is satisfied first by the condition (hy") s=0; 
One is thus led to the class of flows where the streamlines are parallel lines, 


y =y(x, t), or concentric circles, y =y(r, t). In the former case (3.14) is satisfied 
by 


wy (x,t) =c, x8 + 6, x? 4 m, (t) x + F(x, t), (3.5) 


where c, and c, are arbitrary constants, m,(é) is an arbitrary function of time, 
and F is any solution of the one-dimensional heat conduction equation 


Bai Ps (3.6) 
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for A,=y». It should be remarked at this point that the reduction of our problem 
to the classical equations of LAPLACE, Poisson, and HELMHOLTz or to the well 
known heat conduction or wave equations, et al., completes the solution. In 
any particular problem the form of the solution will depend on the domain 
of the independent variables and the initial and boundary conditions to be 
satisfied, while the general solution may be obtained by separation of variables, 
the use of the fundamental solution, or the employment of infinite series or various 
integral techniques. Since these will not be considered here, the reader is referred 
to the vast literature on such initial value and boundary value problems. Equation 
(3.1d) shows that @ is also independent of y (but not necessarily a function 
of time). The corresponding solutions for ® from (3.1c) are, therefore, 


D(x, t) =, x2 + c,x% + ma,(t) + F (x, t) (3.5 a) 
where m3 (t)— 2C,n=f(t) and F is any solution of (3.6) for K,=7, and 
Di ha hn eer + gy iO) HS =r (3.5b) 


The (é) and (k) components of equation (2.1a) show that s ,=s,—0 while the 
(j) component reduces to (py ,—vy ,x),,=1(t) —6q» =s ,; hence 


plo + Q=sly,t)— $A B= [mi (t) — br] y + mg() FAB, (3.50) 

where m,(¢) is an arbitrary function of time and @ is given by (3.5a) or (3.5b). 
For y=yp(r, t), Vi=rd/ér(rd/dr), and equation (3.1a) is satisfied by 

w (r,t) = [Rk 7? + 1, (t)] Inv + ky 7? + G(r, 2) (3.7) 


where G is any solution of 
Gok, (G4 9iG.,)5 (3.8) 


for K,=y. From (3.1d) it follows that ® ,=0, and the corresponding solutions 
for ® are 


@(r,t) =kalnry + kyr? +n9(t) + G(z,t), (3.74) 
where 13(t)—4yk,=/(t) and G is again a solution of (3.8) (K,=7), and 
O(r) =k(4y) 17? +k, Inr+ke, f(t) =—Rk. (3.7b) 


Using (2.1a) in polar coordinates, one obtains, first, 


(pe Vi'y),.= r[ni(t)—4kv] =v71s 5 
or 
s(7, 0, t) = [nj (t) — 4k, 7] 84+ q(7,2). 


Applying the result to the e,-component equation, one obtains 
== tn 1), 
and the pressure is therefore given by 
plo +2 =[m() — 4h] 0+ Sr t(yl,— 1G!) dr—ZA@’,. (3.70) 


In view of the fact that the pressure must be uniform in any region, it follows 
either that j(t)—4k,»-0, in which case the multi-valued potential function Q 
Arch, Rational Mech. Anal., Vol. 6 27 
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must contain a compensating external force, or, if no such force exists, 1.¢., if 
Q is single-valued, 1, (t)= 4%, vt+ m3 (é)- In equation (3.7c) y is the solution (3.7), 
while @ is given by (3.7a) or (3.7b). 

(b) For &,,-+0 the alternate condition specified by (3.4b) is &, 4.=0, 1.e., 
£,=—V?y isa linear function of y, 


Vip = Kgy. (3.9) 
Equation (3.1a) reduces to 
pray, (3.9a) 
and it follows that 
p(x, y,t) =e"! (x,y), (3.9b) 
where H is a solution of 
ZGH=KH. (3.9c) 


G. J. TAYLOR [10] introduced this solution for the nonconducting problem in 
his study of the decay of vortices. Corresponding to this solution for y, since 
(h7") g==0, equation (3.7a) shows that @ is also linear in y, say, 


D(x, y, t) =m, (t) yp + me (ft). (3.9d) 
Equation (3.1c) then becomes 
myrtmy+m=nmVWy +f, 
or, by use of (3.9), 
y= (4 Ky— my, my") p + (f — ms). 


Comparing this result with (3.9a), one obtains 


My (t) = eRe om, (8) = fF 8), 
and 


O (x, y, t) =c et ®! A(x, y) + mg(t). (3.9e) 


® decays according to the magnetic viscosity 7, while y decays according to ». 
In both cases the space variation H (x, y) is a solution of (3.9c). If y=v=k, 
m, (t)=constant =c, m3 (t)=f(t), and 


D=cy+m,(t) =ce*" A(x, y) + m(t). (3.9f) 


Generally, electromagnetic dissipation is more significant than viscous dissipation, 
1.€.,»<y. Estimates based on a simplified kinetic theory approach for an ionized 
gas such as hydrogen (reference [5]) show that »/7=2X10-7a/0, where « is the 
degree of ionization and 9 is expressed in cgs units. For sufficiently low values 
of the density, which will occur in several stellar applications, it is observed 
that »/7->1. From equations (3.2) and (3.9a, d) one readily obtains 


plo +2 = 3 (y?, + yy) — £ Ky (Am? — 1) y? + mg (t). (3.98) 


For y=v=k, m,(t)=c, and if c= 4-4, the nonlinear terms in (2.1.a) vanish identi- 
cally (7=0), and (3.9g) reduces to 


Plo +2 = — 4} (w?, + y*,) + mg(t). (3.9h) 


Self-Superposable Magnetohydrodynamic Motions 389 


(c) If €.,=0, equation (3.3g) becomes 
Cia Orem, (3.10) 


from which it follows, using (3.4b), that (47°) , i &,&, ,=0. This is satisfied first 
by the condition (/;°) g =0, which now refers to the steady flow counterpart 
of (a), namely p=yp(x) and p=y(r). For the former, equation (3.1a) yields 


p(x) = cy x8 +c, x*7 + Cex + Cy. (3.11) 

The corresponding solutions for ® are given by (3.5a, b). Equation (2.1a) now 
reduces to Sy=—6c¥, so that s(x, t)=—6c,yy+m,(t), and 

plo +2 =m, (t) — 6e,y—4A@’, (3.11) 


gives the pressure distribution corresponding to ®(x, t) as given by (3.5a), while 
if B= @(x), equation (3.5b) may be used, giving 

plot Q=c—bavy—dA(kytx +¢;)2. (3.41 b) 
Notice that in the last case the well known Poiseuille flow solution given by (3.11) 
(parabolic velocity distribution) has the usual linear pressure distribution in 


y modified in this self-superposable example by the magnetic pressure, which 
appears as a quadratic term in ~%. 


For y=y(r) equation (3.1a) gives 
DW sky 7 In? Re nr -P hy 7*. (3.12) 


The corresponding expressions for ®(r, t) and ®(r) are given by equations (3.7a,b). 
Note that from (3.12) the vorticity is given by Vi y=4(k, +h.) +4, Inv. If 
P(r, t) is expressed by equation (3.7a), 


riss=—7Vy),=—4hrr', s=—4k,v0+ (7,1). 
S,=4,=1 (y?,—4@*,), and the pressure is given by 
plot Q=—4khvI—341G,4 fr (y?,—Aa®*,) dr. (3.12a) 
If ®(7) is given by (3.7b), this becomes 
plo +2 =— 4hyv0 — 2 A[R(2n) 47 + hy rt}? + f(y? — 4B) dr. (3.126) 


By discussion similar to that following (3.7c), we see that 2 must be multi- 
valued to assure the uniformity of the pressure in (3.12a, b). If this is not so, 
it is necessary to take k,=0, in which case the solution for p(7) reduces to the 
result V7 y=4k,=constant, which will be considered next. 

(d) The second possibility under (c) is €, ,=0, which says that the vorticity 
is constant throughout the fluid, 7.e., 

Vip = 4. (3.13) 

This is a solution of both the steady and unsteady flow problem such that 


p(x, y,t) =a x + a,y? + b(t) xy +V(x,¥,2), (3.13) 
27* 
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where 0(é) is an arbitrary function of time, 2(a,-+a,)=c,, and V is an arbitrary 
harmonic function (V7V= 0), while in the steady flow case 


p(x, ¥) =a + a,y? + a,xy + V(x, 9). (3.13 ¢) 


Equations (3.3e) and (3.4a) give the corresponding solutions for ®. From 
(3.4a) the condition (f;”) ,=0 refers to the previously obtained solutions for 
w(x, t), p(7,t), p(x), and p(7), which now must be modified in view of (3.13a) to 
Ne w (x,t) =c, 27/2 +m (x +m (t), yy.t)=m (t) Inv + n(t), (3.14a, b) 

w(x) = Cy 47/2 + Cy % + Cs, w{r)i=Cgln ¢ + Cg. (3.14¢, d) 


The solutions for ® are given by equations (3.5a, b) and (3.7a, b), while the 
necessary adjustments in (3.5c), (3.7c) and (3.11a,b), (3.12a, b) for the pressure 
are apparent. If ® ,=0, equation (3.3e) becomes ® = (¢y*—G® ,) hi where 
é,=—V2y=—«, and /(t)=—c. Upon substituting this result into (3.4a), we 
obtain as the only new possibility ® ,=c(c,7) 7, or 


DP =c(ey) ++ Ce. (3.15) 


This represents the combination V7? y=c,, Vi B=cy + where p(x, y) and ®(x, y) 
may be expressed by (3.13c) and (3.15) respectively. Determination of the 
pressure depends upon specific forms of the arbitrary harmonic function V(x, y). 
If one specifies c?A=cj7?, the nonlinear terms vanish identically (~=0). 

An alternate assumption for the magnetic field is B=B,. Equation (2.1d) 
requires that B, does not depend upon z. In this case (B-V) B=0, and system 
(3.1) applies with ® replaced by B,, f(¢)=0, and (3.1b) becomes 4 (yp, Vz y)/2(x, y) 
=0. Both the vorticity and the magnetic field are normal to the plane of motion 
and are constant along streamlines. The solutions for y determined above are 
still applicable, while the corresponding magnetic field (as well as the resultant 
values of E and j) now take on a new physical interpretation. This will have 
additional significance in subsequent investigations of particular boundary-value 
problems. Note that an immediate consequence of this assumption is that 
corresponding to a constant transverse magnetic field, B,, one obtains the general 
class of motions wherein the vorticity is everywhere constant or proportional 
to w [equations (3.13) and (3.9)—(3.9c), respectively] in addition to motions 
corresponding to p(x, t), p(7, t), p(x), and p(r), previously defined. 

The solutions outlined above for the steady flow case are obviously valid for 
a perfect fluid as well as one in which the electrical conductivity is infinite. For 
a perfect fluid (v0) in steady flow, p(x, y) need not be biharmonic, as is seen 
from (3.1a). We require only that Vi y=/(y), and the equivalence of these 
self-superposability conditions for plane flow with SToKEs condition for yp to 
be the stream function for steady flow of a perfect incompressible fluid has been 
previously noted. For a viscous fluid with infinite electrical conductivity (7 = 0) 
in steady flow B, (or ®, if f(t) =0) need not be harmonic. The stream functions 
derived above are applicable, while the magnetic field is specified only by the 
relation B, (or ®)=g(y). For an inviscid fluid with infinite conductivity, the 


only requirement in steady flow is that the vorticity and magnetic field be con- 
stant along streamlines. 
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4. Three-dimensional Flow 


The assumption of plane flow in Section 3 requires that both conditions 
U = U(z) and U=0 be satisfied. R. BERKER [11] observed that either restriction 
may be deleted, resulting in two classes of ‘“‘pseudo” two-dimensional motions. 
To study further the character of the present system of equations we shall now 
examine one of these possibilities with regard to self-superposable magnetohydro- 
dynamic flows. Because of the presence of the magnetic field, this type of motion, 
which is characterized in detail by BERKER for the nonconducting case, is found 
to have additional physical significance. 

Cased, U=-U(z), U=20: 

For the same form of magnetic field, it is easily shown that the potential 
functions p(x, y,t) and ®(x, y,¢) can be introduced in a manner similar to 
that of Section 3. The vorticity components, §=V x U, and the components 
of C=V x B become 


§=iU,—jU,—kViy, C=1B,,—Jj B,,--kVi®, 


where Vj is the two-dimensional Laplacian operator. The self-superposable 
system (2.4, 5) can be written as follows: 


(Vip =p, Us=rVU+hAl, (4.1.a, b) 

O(y, Vi p)/8(x, y) =18(®, Hi ®)/a(x, y), (4.2a) 

a(y, U)/e(x, y) =20(®, B,)/e(x, y), (4.2b) 
®,=1KWO +f), Bi=nVi B,, (4.3.a, b) 

aly, P)/A(x,y) =0,  a(y, B,)/(x, y) = 6(®, U)/e(x, y), — (4.4a, b) 


where /,(¢) and /,(¢) are arbitrary functions of time. Equations (4.1) are obtained 
from the reduced momentum equation (2.5a); (4.2) represent the first self-super- 
posability condition (2.4) ; (4.3) are obtained from the reduced induction equation 
(2.5b) and (4.4) from the second self-superposability condition (2.4b). Observe 
that equations (4.1) —(4.4a) are just equations (3.1), so that the functions w(x, y,¢) 
and @(x, y, t) are the two-dimensional solutions derived in Section 3. Equations 
(4.1)—(4.4b) define the corresponding z-components, U, and B,, for each two- 
dimensional solution. The (i) and (j) component equations of (2.1a) are unchanged, 
so that in terms of the obvious 2D and 3D subscript notation it follows that 
S3p,x=Sap,z and S3p y=Szp,y- The (k) component equation reduces to 


U,; + Wy Oe Olea A(@, Be P , Bs) SSeS Seis! Vz U,. 


By (4.1b), (4.2b), this becomes s ,=—/,(f) which, along with the results noted 
above, implies that 


$3 (%,y, 2, t) = — zh (4) + Sen(*, 9 t) + f(t). (4.5 a) 
Since $3 p= (p/o-+2)3 p+ (A/2) (@?, + B, + Bi), it follows that 
(ple + 2)3n = (Ble + Q)on — (A/2) Bi —zf,() +A), (4.5b) 


where (pf/o +2), p is the pressure distribution determined in Section 3. 
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Introducing « and f in the manner previously described, we obtain rate 
(4.1)—(4.4b) the following additional equations: using (3.3a), from (4.2b) w 


have 
U, 


4, 


paA®, Bp, (4.6a) 
and from (4.4b) 
B, p= Di, ig (4.6b) 
Combining equations (4.6), we obtain (1—1®*,) U, ,=0. 
(i) U,,=0. This implies B, ,=0 from (4.6). Thus U, and B, are also functions 
of y, and (4.1, 3b) become 
U, -= 0 (Ay hy) (hg MU, 2), tA) = VM Ua — VEU a tht), (47a) 


BLy= n (ry hy) (he hy* B,a),o= UI hy* Byes mt H&s Be (4.7) 
Applying 0/08 to (4.7) and using (3.6f), we see that 
(11°) U, 42 = 0, Oi) s Pudi 0- (4.8a, b) 


Following the scheme of Section 3, for a representative number of cases the 
corresponding solutions for U, and B, can now be obtained. 


(a) &,,+£0, (27”) ,=0. Corresponding to the two-dimensional solutions given 
by equations (3.5), (3.7) are the following results for U, and B,: 


U,(%, t) =e, x? + 6, x +m, ()) +E (4, f), (4.9a) 
where my; (t)—2¢,v=/,(t), and F is any solution of (3.6) for k,=v. 
Bx = 7 By OR Blt) Cs 4 + Cela). (4.9b, c) 
From (3.5c) and (4.5b), 

plo + Q = [mi (t) — 6e,v] y + s(t) — fd) 2 — (0/2) (@, + BY, (4.94) 
where ® and B, are given by (3.5a) or (3.5b), (4.9b), and (4.9c). 

U,(7,t) =k, lnr +k, +n, (t) + G(r, 2), (4.10a) 
where nj (f)— 4vk,=f,(t), and G is a solution of (3.8) for k,=». 

B,(v,t) =7(B,,,+7"B,,) or B,(r) =kgn7r+k,(n +0). (4.10b, c) 
Equations (3.7c) and (4.5b) provide the corresponding pressure distribution. 
pe rans (11°) +0, €,a.—=0. Equation (4.8a) indicates that U, is also 
U.=(p, (4.11 a) 
in which case for f,(#)=0 (4.1b) reduces to (3.9a). The vorticity components 
 &=4y=4U, &=-ay.=4U, &=—yp=—koU, 


If c.=\V—kg, ks<0, we note that this solution satisfies the condition 


SV Ce On (4.41 b) 
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The vorticity vector has the same direction as the velocity vector at every point 
of the fluid and is proportional to it. Motions of this type (Betrami flows) have 
been treated by several authors [1], [15]. The solutions for p(x, y, t) and ® (x, ¥, 2) 
are given by equations (3.9b,c) and (3.9d—f) respectively. Since B, is also 
hi in y, it foll that 

inear in wy, it follows tha ot Feng (4.410) 


where the analogy with the previously mentioned force-free field condition of 
CHANDRASEKHAR is apparent. 


The remaining cases are readily derived from the results of Section 3. Let 
us refer now to the second possible solution. 


(i) U,,+0, @,=24. Thus, 


OD) = A-*y + mg (t), (4.12a) 
and from (4.6) in 

U=YAB, + g(a,1), (4.12b) 
i.e., U,—VAB,= function of y. Equations (4.1, 3b) become 


U, =o (hh 
Bri= (yh 


2) 1 [(hg ha Uy, x), + (a he” Us) sp) HA), (4.13 a) 

2)? [he hy* Be) lie (Ay hy" B, ~),] : (4.13 b) 
Substituting (4.12b) into (4.13a) and using (4.13b), we obtain 

) 


VAton)*(v — 9) Beet Fa” 8 aa 80+ 7A) =7 8. (4.14) 
Applying 0/08 to (4.14) and using (3.3f) yields 
VA)» — 1) Baar + (i), p80 0- (4.15) 


Representative solutions for U, and 5B, can once again be presented, using the 
scheme of Section 3. 


(a) &,,=+0, (47°) ,=0. We shall restrict ourselves to the case p=y(x, t), O= 
@ (x, t) in the following discussion. Equations (4.1, 3a) become 


YP t= VW axe t pi (t) x + pelt), © .=7© ,,+ h(t), (4.16a, b) 


where #, and #, are functions of time only and where the Jacobian relations 
(4.2, 4a) are identically satisfied. We now require, however, the following addi- 
tional relations: 

U=VaB,+2(x,t), O=A-tyt mill). (4.17a, b) 


Substituting (4.17b) into (4.16b) we have yp ,=7y + VA (f.—ms3), and comparing 
this result with (4.16a) yields 


(7 —») Y2x= Prlt) * + Pelt) — VA Tf () — mad], (4.18) 
(n —») pr=n br(t) x +7 bo (t) — VA [fe(t) — ma(0))- (4.18b) 

Substituting equation (4.17a) into (4.1b) and using (4.3b) gives 
VA (y — ») V2 B= VA? (9 — 9) Br t= 822-8 th. (4.19a) 


Or 
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Applying 0/dy to this result, we obtain 
Van (ny —v) B,y:=0- (4.19b) 
Now if 7=v=hk, equations (4.18) require that b, (t) =0, po(t) =VA(f.— ma), 


that w(x, t) is now a solution of 


YPr=RkYrxt Pot), (4.20a) 
and O=/-!y+m,(t), where 
= VA(fo— ms) (4.20b) 


defines ®(x,t). Equation (4.19b) is satisfied for B, ,,==0, and B,(x, y,t) 1s thus 
a solution of the two-dimensional conduction equation, 


B, »= RV? B,. (4.21 a) 


Equation (4.17a) defines U,(x, y,t) where g(x, t) is any solution of the one- 
dimensional equation 


WE Aa See Op (4.21 b) 


If v--y, integrating (4.18a) twice with respect to x, differentiating the result with 
respect to ¢, and then comparing with (4.18b), we see that 


, (t) = constant = 4, 
— VA[fo(t) — mg (t)] = constant = ag, 
M(t) = aynt+ ag, 
m, (t) =a, t+ VA(n—v) [f f, (d) dt— me (t)]+-a,=ay vt+ (n—v) f po(t)dt+a,. 
Therefore, w(x, t) is defined by 


y= a, x8 + a,x? + ms (t) x + me (t). (4.22) 


This is the special case V7 y=c, to be considered. The corresponding solution for 
(x,t) is given by (4.20b). Equation (4.19b) now requires that B,y1=0; hence 


B,= Hy, (x, t) + Hy (x, y)- (4.23) 


Putting (4.23a) into (4.3b), we see that A, .—1 A, ,.~=nVi H,=H,(x), where 
H, is an arbitrary function of x. Letting H,=H,(x, t)-++h,(x) in the equation 
H, ,— A, ,,=H 3, we obtain (for H;=—ynhi) H, 1= nM, ,, pune es 
H,= H,(%, y) +he(% ) in n V3 Hy = H,, so that for Hn (=O), we havea 


Note that H;=nhz—=—mnhi implies that hg (x )+Ay (x) =cyx-+c¢,. Thus ae 5 
as defined by (4.23a) becomes 


. A, = H, (x, t) +h, (x), (4.23 b) 
where HT, is a solution of 


Hi, 1 = My x2, (4.23 c) 
h, is an arbitrary function of x, and 


H, = H,(%, y) + 1% + ¢,— h, (x), (4.23 d) 
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where H, is harmonic, 7.¢., 


U,(x, y, t) is now defined by equation (4.17a), where g(x, ¢) is a solution of 


VR on — 2 42 Itt) =VAn(n —») HT, ;, (4.24) 
H, (x,t) having been determined from (4.23c). 


(b) €,.-0, (A1") p= 0, & a.=0=@ 4,4. The last equation requires that ® be 
linear in wy. In Section 3 equation (3.9a) permitted the required compatibility of 
the solutions for 7==y» by proper choice of the arbitraty function m, (¢). It should 
be noted that the condition &, ,,—0, implying that the vorticity be linear in y, 
can be satisfied only by (3.9) —(3.9c) since the choice of K, equal to a function of 
time would ultimately lead to a contradiction. In view of (4.12a) we now require 
m, (t)=A~* in (3.94). Substituting this into (4.3a) gives p,=yVi pt V/A (fe—ma). 
For n=v=R, p(x, y,t) is given by 


Wate Baa Ne VE = KA, (4.25) 
while (4.12a) defines ®(x, y, t) with mg(t)=f,(t). For 7=-y the only compatible 


solution possible from &, ,,,= 0 is 
Vi=¢. (4.26) 


Thus, y ,=¢4 VA(fo—m3) and 
p(x, ¥,t) = mg(t) + f(%, ¥), (4.27) 


4 
velocity components are thus independent of time. For y=v=hk, B,(x, y, t) 
is any solution of 


B, ,= kV? B,. (4.28a) 


Equation (4.12b) defines U(x, y,4) where, from (4.15), 24.—=0 implies that 
g=n,(t)p+n(t). In this case £,=—Viy=—Ksy, and y is given by (4.25). 
From (4.14), therefore, for ;=/,, it follows that (¢)=0, so that 


g=cop+ melt), malt) =h(t). (4.28b) 


For v==y equations (4.14), (4.15), which define g, show that for (h7") g 0 it is 
necessary that g ,,—0 and, further, that B, ,,=0. Thus B,, must be a function 
of y, which merely expresses an integrability condition for g in (4.14); while 
B,(x, y, t) is in general a solution of (4.3b), in order that (4.14) define g it is 
further required that B,, be a function of y. Since £,=— Vi p=— ¢, is required, 
equations (4.26) —(4.28) and (4.14) yield, for m,(¢)=c,, the formula 


Vay —n) 49 Bir=lgmg (2 = C9 + fi). 
Thus 
VA(v — 9) 19 B,= cy g(t) + h(x, y), 
but 


oe nin Vi B,=Vih= VA(v —)2n B,,= Con Ms - 
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Hence m3 (¢) constant = ¢s, 1.€., 
ms (t) =e, + VAlfa(t) — ma(t)] = 6s. 


As a result of the restrictions imposed on the solution y, ® for »--7 by equations 
(4.12a) and (4.26) —(4.28), B, is now given by 


B,=nA-4(v — n)* [eo (cgt + ea) + h(x, ¥)], Vih=fo0gee (4.28c) 


(c) €,,=0. For (hj°) ,=0 the solutions corresponding to w(x) and ®(x) are 
readily obtained from (a). The second possibility, called (d) in the corresponding 
discussion of Section 3, is ¢, ,=0, 1.€., 


Vig cy. (4.29) 


The two-dimensional solutions for y and ® are given by (3.13)—(3.15). In 
view of (4.12a) we now require that V?®=A3° V7 p=c, A+}, 1.€., ®@A=cin? in 
equation (3.15). For y=y=h, g is linear in y, which is defined by (4.29) rather 
than by (4.25). B,is given by (4.28a) and U, by (4.12b). For »==y, the discussion 
in (b) applies. 

Case 2. U+U (8), U,=0. 

The analogous problem of a “‘pseudo’”’ two-dimensional motion with axial 
symmetry is now considered. Magnetohydrodynamic flows of this kind have 
already been studied somwhat, apart from questions of superposability. The 
basic equations will therefore be introduced in this section but, for the present, 
the discussion of the solutions of these equations will be restricted to several 
immediate results. 

Potential functions p(y, z,t) and ®(r, z, t) can be introduced so as to satisfy 
the solenoidal condition upon U(r, z,¢) and B(v, z, t). Thus 


U=e,r*y ,.+eU,—e,71yp,, B=e,r1®,4e,B,—e,7'@,, 
and 
§=—e,U,,+e1* D'y +e,77 (7G) ,, 
where 
D2 = 8/87 — 1 afar + 2/022. 


System (2.4, 5) becomes 


(Dy) .="D'y,  (rU).=r DU) +A), (4.30a, b) 
(yp, 7-2 D¥y)/0(r, 2) +27 U, U, = A[2(, 2 D? O)/d(r, 2) +27 By By ,], (4.300) 
a(p,7 U(r, 2) =2.0(®,r By)/O(r, 2), (4.304) 

® =n DD +f,(f), (7 By) = D*(r By), (4.30¢, f) 

a(y, ®)/2(r,2)=0, 8(y,7 2 B,)/8(r, 2) = 0(®, 1 Y)/a(r, 2). (4.308, h) 


For the magnetohydrostatic problem (U=0) CHANDRASEKHAR (([8] and [9]) 
examined the case of a force-free field which, in general, is characterized by the 
equation Bx curlB=0. As previously observed, this relation is a specialization 
of the general self-superposability equation derived earlier. It is examined in 
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some detail for the axially symmetric problem presently being considered, and 
the reader is referred to these papers for the solutions when U=0. A non-zero 
velocity is accounted for in [12], where, under the additional restriction that 
the fluid be inviscid, a comprehensive treatment of the full axisymmetric equa- 
tions is given. Included are several interesting solutions of the full system, the 
derivation of a number of general expressions and integral relations, and a dis- 
cussion of boundary condition considerations. In view of the obvious need for 
examining system (4.30) in its own right a detailed analysis will not be attempted 
at this time. Simplifications for the steady flow case are apparent, while some 
rather general remarks relevant to the corresponding two-dimensional case 
(U,=0=B,) and to the immediate solution U=|/AB (note [12]) are better 
left for a detailed analysis. For the purpose of illustration, however, let us 
consider the following particular examples. 

(i) =@(r). Equation (4.30g) shows that y=y(r). The operator D? now 
assumes the form D?=r0/dr(r1d/ér), so that [for y=0, (4.30a) is identically 
satisfied and the following is not required] D+y=0 provides, after several inte- 
grations, the well known Poiseuille motion: 


y= Inv + cor + c37?, (4.31) 
where the c; are arbitrary constants. From (4.30e), taking k;=—/f,(¢), we have 
D? @=k, 74, so that 

® =k, (2y) 477 Inv + ¢,7. (4.31 b) 


From (4.30d), y,U;,=A‘®, Bs, while, from (4.30h), py,By.=@,U,,. It 
follows that 
Py (iis = A Byte Us, = U5 Bs’. 


Clearly, there are two possibilities involved, namely, U,,=By,=0 or U; ,= 
V7 Bo. ¢- In the former case, all the Jacobian relations (4.30c, d, g, h) are identi- 
cally satisfied, and complementing equations (4.31a, b) are similar expressions 
for rU, and 7B, as obtained from the relations D?(7U,) =k, (ky= — f(t), 
D? (7 By) =0 (n+0) using D?=rd/ér(r4d/07). Alternatively, U;, ,= VABy. and, 
therefore, yp p=VA ®,. From the latter it follows that p= @-+- constant, 
whence it is now required that c,=h, Va(2n)2, Co=0, Cg=VAcy. In addition, 

Uy = VA By + g(r). (4.32) 


Substituting (4.32) into D?(rU,) =k, 1, we obtain 
D2(r VA By tre) =VAD?(r By) + D2(r g) = D? (rg) = hav, 
since, from (4.30f) for 7-40, D? (ry By) =0. Substituting (4.32) into (4.30c) yields 
271U, Us; ,= 274 (VA Bo +8) VA Be; 
=2A7r1 By By, +2VArig By = 2dr By Bos. 


Since By ,--0 for this case, we must have g=0 (and therefore k,=0), so that 
U,=VAB, where By(r, z) is a solution of D?(r By) =0. 
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(ii) O= P(r, t). The previous discussion pertaining to the equations contain- 
ing Jacobians still applies. We require, therefore, that y=y(r,t) and that 
either U; ,=By,—=0, in which case all the Jacobian relations are identically 
satisfied, or Uy .=£ 0, By ,=-0. In the former case the solutions for y, ®, U, and 
B, are readily obtained from (4.30a, b, e, f), using D?=rd/0r(r44/6r) to express 
the resultant one-dimensional nonsteady flow equations. For U, ,=-0, By ,=0 
we obtain once more CAE, and y,=Va @ ,. From the latter, y=VAG+ 
h(t), and the solution for y and ® from (4.30a) and (4.30e) is thus restricted in 
a manner already experienced, 7.e., depending on the constants 7 and y. In 
addition, U,=VAB,-+g(r, t) is specified by (4.30c) in that g=O0 is required 
for Bs ,= 0. Referring to equations (4.30b, f), a consistant solution is obtained 
only when »=7 and f, (¢)=0. 
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Monotoniesatze bei Differentialgleichungen 2.Ordnung 


H. H. GLOISTEHN 


Vorgelegt von L. COLLATZ 


§ 1. Randwertaufgaben bei gewéhnlichen und partiellen 
Differentialgleichungen 
B sei ein beschranktes offenes Gebiet des m-dimensionalen Raumes der Ver- 
anderlichen (%,, X), ..., %,)=(,) mit dem Rand J’ und B= B+T. u(x,)€K soll 
heiBen: w besitzt in B stetige Ableitungen bis zur 1. und in B stetige Ableitungen bis 
Oru 


zur 2. Ordnung einschlieBlich. Zur Abkiirzung setzen wir u;= oe und u,;,= 


In B sei der Operator j OX; OX, 
(14.1) Tu=—a,,-u;,+F(x;,u,u,) und auf IP: Uw =u 0 oe (u =1,2,...,m) 


definiert. Dabei wurde iiber doppelt auftretende Indizes von 1 bis m summiert. 
Die symmetrische Matrix a,,(%;)=4,;(%,)* sei fiir jeden Punkt aus B positiv 
definit, «20 auf J’ und 6/do0 die Ableitung auf dem Rande in Richtung der 
Konormalen (s. CotLatz [1], S. 23). Dann ist bekannt (Corratz [1] und [2], 
MEYER [3]), daB unter der Voraussetzung: F(x;, , u;) monoton nicht-fallend 
hinsichtlich w in einem Bereich S des (%;, 4, u;)-Raumes mit (x;)¢ B der Mono- 
toniesatz gilt, d.h. fiir zwei Funktionen u(x,), v(x,;)CK folgt aus Tu< Tw in 
B und U,uswu,v auf I stets wv in B. Die Funktion F war dabei als stetige 
Funktion ihrer Argumente vorausgesetzt worden. 

In diesem Paragraphen soll gezeigt werden, da der Monotoniesatz auch fiir 
gewisse hinsichtlich w monoton fallende Funktionen F Giiltigkeit behalt. Dabei 
wird der Beweis des bisherigen Monotoniesatzes als bekannt vorausgesetzt, so 
da8 wir uns an vielen Stellen im folgenden kurz fassen k6nnen. 

Wir treffen folgende Voraussetzungen: 

(I) F(x;, 4, u;) besitze in S stetige partielle Ableitungen nach w und u,. 

(II) Fiir fest gewahlte stetige Funktionen B,(x;)* und C(x,) existiere eine 
Funktion o(x,;)€ AK mit 


(1.2) p90 in B, L*p=— 4,,-9;.7 B;- 9, +0920 inB, U,p>0 auf I. 
(III) Es sei die Funktion 
(1.3) G(%,,2) =F (4,,9 +. 2, 0;-+ 9, 2) —F(%;, v, v;) — B.-g, -2—C. 9-220 
fiir z=0 mit den in (II) gewahlten Funktionen 9g, B;,C und eine fest 
gewahlte Funktion v(x,) CK. 


* Hier dienen die Indizes nur zur Unterscheidung der Funktionen. 
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Satz 1. Unter den Voraussetzungen (1) bis (III) und Tu< Tvin B, U,usu,v 
auf I’ fiir zwei Funktionen u, v€ K, wobei v die in (1.3) gewahlte Funktion bedeutet, 


gilt uv in B. 
Beweis. Fiir die Funktion z=" definieren wir den linearen Operator 
PY 
(1.4) L2=— Q-Ajy hp — 2° Gn Qj HELD % in B. 
Andererseits erhalten wir mit u=v+-2, U;=U;+ 9° 24+; °4 und 
= 4;,(4— v) 5h =— D:4jn 2jp— 2° Aj Pj Ae — Gro Pj. 2 
—Tu— Tv—F(x,, u, uj) + F(%;, 2, ¥;) 


Lz=Tu—To F(x;,0 + 2, 0, + 92; + 9; 2) 4 

+ F(x;,0,0;) + (By- pj +C- @)-2. 
Auf I" wird U,(u—v)=U,p-2—-«@ Pn 0: Sei nun a= Max aia), das im 
Punkte P,=(%);) angenommen werden moge. Ist 7<0, so ist der Satz 1 richtig. 
Nehmen wir also an, es gelte 7>0. Dann kann Fj nicht auf J’ legen, denn dort 


be Oz . a:p dz 
pee ee S -— 
ware (3-|,s 0 und somit (z)p< & o 86 


iiber p und «. Fir A,B besitzt z ein differenzierbares Maximum mit (z,;)p=0, 
d.h. (4) p= (vj+9;-2)p, und — (a;,-2;,)p,20. Dann gilt fir den Operator Lz 
nach (1.4) mit Z>0 und der Voraussetzung (1.2): (Lz)p=0 und nach (1.5) mit 
Tu—Tv<0in Bund der durch (1.3) eingefiihrten Funktion: (Lz) p= —G(%p;,4) 20 
wegen Z>0. Diese beiden Ergebnisse widersprechen sich, d.h. die Annahme 


Z>0 ist falsch, es gilt z= Z%<0 in B. 


(1.5) 


le <0 auf Grund der Voraussetzungen 


Wie aus dem Beweisgang hervorgeht, gilt der 


Zusatz. Der Satz 1 bleibt unter der Voraussetzung Tux<Tv in B, U,us U,v 
auf I’ richtig, wenn statt (1.2) 


(ig2*) p>0 und L*p>0im B, U,p>0 auf l 
oder statt (1.3) 
(beta) G(a,,.2) = Oui oz 0) 


gefordert wird. 


Die Hinzunahme des Gleichheitszeichens in Tu< Tv in B kann bei Erfillung 
von (1.2) und (1.3) auch unter gewissen Voraussetzungen iiber die Funktion 
F(x;, 4, u;) bewiesen werden (s. z.B. den Zusatz zum Satz 3), doch soll dieses 
hier nicht weiter durchgefiihrt werden, da man in den meisten Fallen mit (4.2%) 
oder (1.3*) auskommen diirfte. 


Bemerkungen zum Satz1. Die Anwendbarkeit des Monotoniesatzes hangt 
davon ab, ob es gelingt, Funktionen 9, B;, C so zu bestimmen, daB die Voraus- 


setzungen (II) und (III) erfiillt sind. Fiir einige Sonderfalle soll dieses im folgenden 
untersucht werden. 


| 4. Sei F,(x;, 4, u;)20 in dem hinsichtlich « konvexen Bereich S, so lassen 
sich die Voraussetzungen (II) und (III) mit p=1, B,=C=0 erfiillen, denn 
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dann wird G(x;,z)=F(x;,v+2, vj) —F(x;, v, v) =E,(x;,0, v,)-z20 fir z=0 
(v bedeutet hier eine Zwischenstelle zwischen v und v+z). Wir erhalten also 
den Monotoniesatz in der bekannten Fassung. 


2. Fiir lineares F=—b;-u;+¢ -u (wobei der Fall c=O bereits durch 1. 
erledigt ist, wir lassen also auch negative c-Werte zu) wahlen wir B,=— b; 
und C=c und erhalten damit G(x,;, z)=0. Wenn sich zeigen la8t, daB irgend- 
eine Funktion y(x,) mit den Eigenschaften —a,,-;,—0;-y;+c-p20 und 
y>0in B, U,~>0 auf I’ existiert, so 1aBt sich der Satz 1 anwenden. In vielen 
Fallen gelingt es sehr einfach, solch eine Funktion m anzugeben (s. hierzu die 
Beispiele). Besitzt die Eigenwertaufgabe —a;,-w;,—b;-w;+4-c-w=0 in B 
und U,w=0 auf J’ den kleinsten Eigenwert A, und ist 1,>1, so gibt es eine 
zugehérige Eigenfunktion w,(x;)>0O in B. Diese Funktion erfiillt dann im 
wesentlichen (man mu8 wegen der Forderung U,y>0 auf I’ noch eine Grenz- 
betrachtung durchfiihren) die Voraussetzung (1.2). Aber der Nachweis j,> 1 
wird wohl haufig nicht einfacher zu erbringen sein, als die direkte Konstruktion 
einer Funktion . 

3. Fir die in w nicht-lineare Funktion F=—b,;-u,;+/(x;,u) setzen wir 
B,=— 6, und C= F,(x;, v)=f,(%;, v) und erhalten 

G(4;,2) =f(%,,0+ 92) — F(X.) — ful %i Y) 9° 2 =F fu (9) » (YZ), 
wenn wir die 2. partielle Ableitung von f nach w in einem beziiglich w konvexen 
Bereich S voraussetzen. Ist dort f,,=f,, 20, so wird die Voraussetzung (III) 
G(x,;, z)20 fiir z=0 erfiillt. Zum Beispiel wird fiir f=—v® stets f,<0, der 
Monotoniesatz in der bisherigen Fassung ware nicht anwendbar. Guilt jedoch 
u=0 in S, so ist f,, 20. Die Voraussetzung (II) verlangt 


(1.6) = Asn Pix — 93° Gi + E(%;, 0) “G20 


in B mit U,m>O0 auf J’ fiir eine in B positive Funktion gy. ZweckmaBig wird 
man hier eine Minorante bestimmen: F,(x;, v)=—w*=const und statt (1.6) 


(1.6') — 45° Dix — 0; Yj — O*- PZO 
in B untersuchen. 

4. Ist w eine Funktion der einen Veranderlichen x,=%, so haben wir den 
Operator 


(1.7) Tu=—u"'+F(x,u,u') fir a<x< a, (=+) 


und die linearen Randbedingungen 
(1.8) U,u=u(a,) +(—1)*a,-w'(a,) mit «,20 (V—ie2)s 

Die folgenden einfachen Beispiele mégen zur Erlauterung des Satzes dienen. 
Bei der Durchrechnung kam es nicht auf sehr groBe Genauigkeit an — die lieBe 
sich leicht durch Hinzunahme weiterer Glieder in den Ansatzen steigern —, 
es sollte vor allen Dingen die Einfachheit des Verfahrens gezeigt werden. Alle 
Beispiele lassen sich nach dem bisherigen Monotoniesatz nicht behandeln, man 
war bei ihnen auf andere Methoden angewiesen, z.B. Iterationsverfahren, Ritz- 
sches Verfahren usw., bei denen Fehlerabschétzungen miihevoller sind und 
meistens auch schlechter ausfallen. 
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Beispiel 1. Lineare gewohnliche Differentialgleichung (s. CoLLatTz [J], S. 183) 


(1.9) T= Ul VE) eA fie en u(+1)=0. 


Hier ist c=—(1-+*2)<0, um den Satz 1 anwenden zu kénnen, miissen wir 
zunachst die Existenz einer Funktion g(x)>0 mit den Eigenschaften L* p= 
—p" —(14+ 22) p=0 fiir |x|<4 und p(+1)>0 nachweisen. Fur p=1-—3-# 
wird L*¥p=4(1—x#+2%4)>0 (die Funktion g=cos 2x wiirde ebenfalls die 
verlangten Forderungen erfiillen). Der Satz 1 kann damit auf den Operator Tu 
aus (1.9) angewendet werden. — Als Naherungslosung fiir « machen wir einen 
Potenzreihenansatz, der die Randbedingungen bereits erfiillt: 


v(x) =a+bx?+cxt— (a+b+ 0c) x8. 
Es wird dann 
Tu—Tv=(4+a-+ 20) + (a+ 6-4 12¢) — 
(30a + 29b + 29c) x4— (a+ b) x8— (a+b+¢) x8. 


Aus der Forderung, daB die ersten Potenzen bis x4 einschlieBlich verschwinden 


sollen, folgt eine 1. Naherung w= (319 — 331 x2+ «4 —11x°) mit Tu— Tv, 
=e (124+ 11x4)20, also vy; fiir |w|<1. Aus der Forderung Tu—Tv=0 


fir, x=0 und *=1, (Tuy—Tv)"=0 fir, x=@.ergibt, sich eine 2, Nahering 


PE) Saker, = Bile 4 Soe’ Peed ie os 2 4) < 
= 366 (342 354 «2+ x4 —11%°) mit Tu— Tv, a 12x? —11x4)<0, 


also uXv, fiir | x| <1. Insgesamt erhalten wir somit die Schranken 


Ug 


sag Ol BNE 41%8) S u(x) Se (342 — 354x2-+ x4 — 11x) 


0) fara 4" 


Zum Beispiel fiir x=0: 0,9300 S u(0) S$ 0,9345. 
Beispiel 2. Nichtlineare gewohnliche Differentialgleichung 


(1.41) Ta=—' —u4—4eV=F% tir |x) <1, u(+1)=0. 


Den Bereich S des (x, u)-Raumes legen wir fest durch |x|<1, 0S w<U=const. 
Dann wird in S: f= —1— Zu=—1—2U=—@? und ¢(x)>0 ist nach (1.6’) 
aus der Bedingung L* y= — y"’ —w?y>0 zu bestimmen. Fir p=1—t- x? mit 
tT<1 wird L* p=2t—w*(1—1x?)=21t—w?>0 fiir w2< 217< 2, dh. der Satz1 
ist anwendbar, wenn uw U<# in S bleibt. Diese Bedingung muB hinterher 
aus den errechneten Schranken tiberpriift werden. Der Ansatz v—a (1— x?) 
liefert To=alt 1. 42 ; (1 2), Die Forderung Tv<Tu= + liefert den 
Wert a,=7 und Tv=>Tu=} den Wert a,=0,6340. Somit haben wir also zu- 
nachst die groben Schranken 


0,25 (1 — x?) < u(x) < 0,6340(1 — x4) <0,6340 =U <8. 


Um bessere Schranken zu erhalten, setzen wir v= (a+ 6x2) (1— x*). Verlangen 
wir Tu—Tv=0 fiir x=0 und x=1, so erhalten wir a,=0,4845, b,= — 0,0469 
und mit diesen Werten Tu— Tv<o0. Aus Tu— Tv— (Lu— Tv)” =0 fir x=0 
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folgen a,=0,4580, 6,= —0,0560 und Tu — Tv=0, also insgesamt die Schranken 
{0,4580 — 0,0560 x”) (1 — x?) < u(x) < (0,4845 — 0,0469 x?) (1 — x?) fir |x| <1. 


Bei Hinzunahme weiterer Glieder steigt der Rechenaufwand auBerordentlich 
rasch an, da dann nichtlineare Gleichungssysteme zu lésen sind. 

Bei gewohnlichen Differentialgleichungen wird man im allgemeinen die Nahe- 
rungslésung v(x) so wahlen, daB die Randbedingungen erfiillt sind. Im Gegensatz 
dazu ist es bei partiellen Differentialgleichungen haufig méglich, von Naherungs- 
lésungen auszugehen, die die Differentialgleichung exakt erfiillen. Bei den linearen 
Differentialgleichungen gelingt es, eine einfache Fehlerabschatzung aus den 
Defekten auf dem Rande anzugeben. Wir betrachten die 1. Randwertaufgabe 


Cf 12) Lu=— a;,:4;,—6;-u,tc-u=r(x,) inB, u=f auf I. 
Die Funktion g(x,) hat dann die Bedingungen 


(1.13) i o=—Lloeoin Bb wid. o> 0 aut 
zu erfiillen. 


Satz 2. Ser v(x;)\CK eme die Differentialgleichung (1.12) exakt erfiillende 
Naherungslésung fiir u(x,) und nimmt (u—v) auf IX den Wert Null an, so gelten 
in B die Schranken 


(1.14) p(%) Min “—* < u(x) — v(x) S p(x) Max “—" 


mit einer Funktion py, die den Bedingungen (1.13) gentigt. Ist Lp=0 in B, so 
gelten die Ungleichungen (1.14) unabhdngig davon, 0b (u—v) auf I’ den Wert Null 
annimmt oder nicht. 


; is : 5 _ es uUu—v a 6 
Zum Beweise definieren wir wieder den Operator Lz fiir z= —— wie in 


(1.4) und erhalten mit Tu — Tv=L(u—v)=0 und G(x,, z)=0 nach (1.5) Lz=0. 
Auf diesen Operator (der linear ist mit dem Koeffizienten L@ von z) kénnen 
wir wegen Ly=O in B die bekannte Fehlerabschatzung fiir die 1. Randwert- 
aufgabe anwenden (s. CoLLatz [1], S. 375), d.h. es gilt Min 22S Max: in B, 
woraus sofort (1.14) folgt. 

Zu beachten ist, da& wir bei der Anwendung des Satzes 2 die Funktion p 
numerisch verwenden, wahrend beim Satz 1 die Existenz einer solchen Funktion 
ausreichte. Man wird daher bei der Aufstellung von Fehlerschranken méglichst 
sorgfaltig bei der Wahl der Funktion » vorgehen. 


Beispiel 3. Lineare partielle Differentialgleichung in zwei Veranderlichen 


LO=—A4 -u—1 tur |xl<4, |yj<1 
(1.15) : 
und #=0 fur |[x[/=1 und |yf=4. 
Wir wahlen (x, y)=cos 3 - COS 3 >0 in B mit Lgp=0 (es lassen sich leicht 


weitere Funktionen gy angeben, doch die hier gewahlte diirfte die giinstigste 
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sein). Spezielle Lésungen der Differentialgleichung (1.15) lassen sich leicht 
bestimmen. Wir setzen die Naiherungslésung v(x, y) in der Form 


B I Tt 

v(x, y) = — 1+ A(cos x + cos y) + oe [cosh Edin COS 29 + cos Sta cosh ty) 
mit r=Ji- n?—1=1,2113 an, wobei die Symmetrieeigenschaften der Aufgabe 
bereits ausgenutzt wurden. Die Konstanten A und B bestimmen wir durch 


Kollokation so, daB v auf I” ,,méglichst gut Null‘ wird. Die Funktion 


peel nolo yet) ee ge 4 x) — Bcos = x 
w (x) p(x,1) p(x, 1) : i og 2 


ist dann fiir die Fehlerabschatzung ausschlaggebend. Da (x, 1) den kleinsten 
Wert bei x=1 annimmt, werden wir die Kollokationspunkte ebenfalls nahe an 
den Rand x=1 heranlegen. Aus der Forderung w (0,5) =w (0,9) =0 erhalten wir 
A=0,91814, B= — 0,42682 und fiir diese Werte w,,3x< 0,01660, Win > — 0,00784, 
also in |x| <1, |y|<1 die Fehlerabschatzung 


— 0,0079 - p(x, y) S u(x, y) — v(x, y) S$ 0,0166 - p(x, y) 
mit 


v(x, y)=—1+0,91814 (cos x+cos y) —0,23351 (cosh r x- cos 5 teos 3 *coshty). 


Speziell hat man fiir ~(0, 0) mit v(0, 0) =0,36925 und (0, 0)=1 die Schranken 
|w(0, 0) — 0,3736| <0,0122, d.h. der Naherungswert 0,3736 ist auf etwa 3% 
genau. — Durch eine andere Wahl der Kollokationspunkte lheBe sich dieses 
Ergebnis wahrscheinlich noch etwas verbessern. Bei dem hier gemachten Ansatz 
mit nur zwei freien Parametern A und B 1aBt sich natiirlich der Satz 1 direkt 
anwenden, indem man A und B so bestimmt, daB einmal v, (x, 1) <0, ein andermal 
V,(%, 1) 20 wird. Wir wiirden dann A,=0,88098, B,= —0,35701, 4,=0,92544, 
B,=—0,42541 erhalten und damit |w(0, 0) —0,3784| <0,0070, also fiir den 
Naherungswert 0,3784 einen Fehler von etwa 2%. Die direkte Anwendung des 
Satzes 1 macht aber bei einem Ansatz mit mehr als 2 freien Parametern gréBere 
Schwierigkeiten, so daB wir dann doch bequemer (1.14) heranziehen, dieses 
besonders, wenn z.B. eine elektronische Rechenanlage zur Verfiigung steht, mit 
der das durch Kollokation erhaltene lineare Gleichungssystem gelést und die 


ELE a u—VU . 
minimalen und maximalen Werte von ——~ auf dem Rande ermittelt werden 
k6nnen. ? 


§ 2. Anfangswertaufgaben bei gewéhnlichen Differentialgleichungen 
Wir betrachten den Operator 
(2.1) Tu=w"+F(%,4,u) in w~n<x*45X Gra 
x 


mit der in einem Bereich S des (x, w, u’)-Raumes stetigen Funktion F. u(x)CK 
bedeutet: w(x) ist in x»<«*<X mindestens einmal und in Xy< xX mindestens 
zweimal stetig differenzierbar. 
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Satz 3. Fuiir eine festgewdhlte Funktion v(x)€K existiere eine in xySx<X 
stetig differenzierbare Funktion «(x), so daB 


(2.2) G(%, 2) = F(x, v4 2,0’ — az) — F(x, v, 0’) — («’ — a2) 220 


wird fiir 20 und xy< xX. Gilt Tu<Tvinxy< «SX, u<vundu'+au<v'+av 
fiir x=%X9, so gelten diese Ungleichungen in xySx<X. 


Beweis. Fiir die Funktionen z(x%)=«(x) —v(x) und Z (x) =2'(x)+a(x) - z(x) 
wird 2(%))<0O und Z(x))<0. Angenommen, es gilt nicht z<0 in x%<x<X, 
so existiert ein Wert x; mit x)< x,=X, so daB z(x,)=0 und z(x)<Oftirs,<%< x,, 
woraus Z (%,) =2'(x,) 20 folgen wiirde. Wegen Z(x))<0 und der Stetigkeit von 
7(x)wexistiert, dann, ein), Wert)(A_,mit 1% <%= 4,5 Z (%_)=0. und Z (x) <0) fiir 
Noa A< Myyod. h.tZ' (%)=0., Mit Z’=(u—v)" +a’ ta’z, 2 =Z—az, w=v-+2 
und w4’=v’+Z—az wird 


“(x)= Tu—Tv—F(4,0+2,0 +Z—a2z)+ F(x, 0,0’) +402 + (a’ — @) z. 


Hieraus Jolet fur x= wegen? fu— Jv<0 und 2(x)S04n +S %*5%,=%, 
Z(%)=0 und mit der durch (2.2) eingefithrten Funktion: Z’(%)<— G(x, z)|,_,, 0. 
Dieses steht im Widerspruch zu Z’(x,)=0, es existiert somit kein Wert x,, fiir 
den 20, wird. AlseZ.(y)-20-in 47, = 7S, eretbt wich danneZ(7,;) = (x,) <0, 


d.h. auch z(%) kann in %<%<X den Wert Null nicht annehmen. 


Zusatz. Ist (2.2) fiir alle Funktionen u aus einem abgeschlossenen Bereich S 
des (x, u, u’)-Raumes erfiillt, und ist F(x, u, u') in S stetig differenzerbar, so folgt 
aus TusTv im xm<xSX, uso und w’+ausv'’+av fiir x= x die Giiltigkett 
dieser Ungleichungen in Xx» xX. 

Der Satz 3 ist also auch giiltig bei Einschlu8 der Gleichheitszeichen in den 
Voraussetzungen. Zum Beweise setzen wir w=v+e-e?"—*” mit e>0. Fir 
x=, wird dann: uwSv=w—e<v,u'+audv'+av=w'+aw—e(d+a)<w'+aw 
fiir A> \ eae Fiir Tv erhalten wir nach einer kurzen Rechnung Tus Tv= 
Tw—ee’ *) (424 FF +1, ,), wobei die partiellen Ableitungen von F an einer 
Zwischenstelle zu nehmen sind. Wahlen wir A>—«(%») so groB, daB die eckige 
Klammer poe wird, so gilt Tw< Tw und nach Satz 3 folgt dann w<w und 
w+tau<w'+aw in %<x*<X. Hieraus ergeben sich fiir e—0 die Behauptungen 
des Zusatzes. 

Im folgenden wollen wir voraussetzen, daB F(x, u, wu’) in S die Voraussetzungen 
des Zusatzes erfiillen mége. — Wir betrachten die Anfangswertaufgabe 


(23) Tu=u"'+F(x,u,w)=r(x) fir %<%, u(x%)=u und Wi X;,) == Wee 
Bei der Konstruktion von Naherungslésungen v(x) wird man im allgemeinen 


von Funktionen ausgehen, die die Anfangsbedingungen exakt erfiillen. Gelingt 
es, zwei Funktionen v,;(x) (7=1, 2) so zu bestimmen, dab 
(24) Tv, TusTv fir <4 X,  — 0;(%) = uM und 0;(x)) = uM 
gilt, so erhalten wir bei Erfiillung der Bedingung (2.2) Fehlerabschatzungen fir 
die Funktion w(x) und ihre Ableitung: 

alo vy—a(v,—v) Sw Sv+a(vp—%) 

<u i 

(25)° tSwsv, und fir bee pas eal 


235 
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im Intervall 4,S%*%SX. Die Schranken fiir die Ableitung wv’ werden also am 
giinstigsten sein, wenn wir «(7)<0 wahlen kénnen. Es soll jetzt untersucht 
werden, wann dieses méglich ist. Nach (2.2) erhalten wir (die partiellen Ab- 
leitungen sind hier an einer Zwischenstelle zu nehmen) : 

G(«,2)=(E,— E,a+oe—a')z20 fir 250, 


uU 


wenn in S die Bedingung 
(2.6) Fo— F,-a+oe2—a'’S0 


UuU 


erfiillt ist. 

1. Fall. F.<0 in S, dann kann «=0 gewahlt werden, die Anwendung des 
Satzes 3 ist hier besonders einfach. Allerdings diirfte dieser Fall bei Anwendungen 
aus Physik und Technik (z.B. nichtlineare Schwingungen oder der auf Zug (bzw. 
Druck) und Biegung einseitig eingespannte Balken) keine allzu groBe Bedeutung 
besitzen. 

2. Fall. F.<w* in S mit w?>0. Hier kann « nicht mehr gleich Null gesetzt 
werden. Wegen (2.5) werden wir méglichst «(x)<0 wahlen. Wir wollen uns 
bei der Diskussion auf den Sonderfall F,,=0 beschranken. Dann ist (2.6) erfiillt, 
wenn « der Bedingung 


w(x — x) — F+G| 


(2.7) w*%+a?—«a’=0 geniigt, also a(x) =w-tan 
mit einer Integrationskonstanten C,>0. Diese Funktion ist nicht mehr durch- 
weg stetig fiir = x), d.h. wir erhalten zunachst nur Aussagen in einem gewissen 
Intervall 45.45 45 %+ =. (% —C,) mit 0<C\<a oder, wenn wir uns auf ein 
a) 
Intervall beschranken, in dem « (x) <0 bleibt, x,<%»)+ — (F = C,) mit O< Gy = y 
In diesem letzten Fall erhalten wir die Schranken v,SuSv,, vySu'<vg in 
X%pSxLx,. Jetzt kénnte man ein nachstes Intervall 4,5 xl x, x,+ £6 (F — C,) 
12) 2 
betrachten, in dem a, (*)=q@ - tan |w(x — x,) — = + q gewahlt wird und be- 


stimmen die Funktionen v,,(«) aus 


(2.8) los TVs Tt, in a, eee 
Vj40(%1) =9,;(%) und v,2=v; + (— 1)foy- (v, — v4) 


fiir *=%,, die dann die Schranken y,<u<y, und v3Su'/<y in 4<x<x, 
liefern. In dieser Weise kénnte man fortfahren und erhielte dann Schranken 
in einem Intervall xy< *< X. Praktisch jedoch diirfte diese Methode den meisten 
der anderen Methoden zur Fehlerabschatzung (z.B. Iterationsverfahren bei 
Anfangswertaufgaben) unterlegen sein. 

Eine physikalisch anschauliche Deutung erhalt (2.6) im Linearfall F=2bw!--cu 
mit 6 und c konstant. Dann setzen wir « ebenfalls konstant und erhalten 


(2.9) a*?— 2batcso. 


Die linke Seite dieser Ungleichung ist gerade die Frequenzgleichung, die durch 
den Ansatz w=e~** erhalten wird. (2.9) l4Bt sich fiir c+0 nicht mehr durch 
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reelle x-Werte erfiillen, falls b?— c?<0 ist. Dieses ist der bekannte Schwingungs- 
fall. Es ware méglich, daB sich aus (2.6) entsprechend ein Kriterium fiir das 
Eintreten eines Schwingungsvorganges gewinnen l4Bt, doch hieriiber sind — aus- 
gehend von diesen Entwicklungen in §2 — noch keine weiteren Untersuchungen 
angestellt worden. . 


Beispiel 4. Lineare Differentialgleichung 
(2.40)  Tu=(2— nu" -u=3— 2% fir 0%, u(0)=w'(0)=0. 


Tw ist hier anders definiert als in (2.1), wo der Koeffizient bei u’’ gleich 1 gesetzt 
wurde. Da aber 2—x?>0 in 0Xx< 2 ist, gilt der Satz 3 entsprechend in 
diesem Intervall fiir den hier definierten Operator. Wir wollen uns im folgenden 
auf das Intervall OX %<1 beschranken. Mit R,=>5- 3! wahlen wir «(x)= 


tan(~—1). Dann ist die Bedingung (2.6) erfiillt, und wir kénnen die in (2.5) 
angegebenen Abschatzungen fiir «<0 benutzen. Der Ansatz v (x) =ax?+bx8+c¢ x4 
+dx°-+ex® erfiillt die Anfangsbedingungen und liefert 


Tv — Tu= (4a — 3) + (120 -+ 2) x + (24c — a) x* + (40d — 50) x8 4+ 
+ (60e — 11c) x4 — 19dx® — 29e x® 


| ee oe Ate alien ee ln pee | Fle 
Caer init ¢ ==, b= — 35 = 55, b= — 7s: 


Tv —Tu= x [(60e — $5) + 2% — 29e x7]. 


fe 5 os 
Fiir a meee as (aus Tv —Tu=0 fiir x =1) folgt 
*  #f 1323" 5, pe tAD el 
Tv,—Tu re (- Sa 19x - a5 x) S0 
und fiir e, ee [bestimmt aus po gu O fiir x > 0): 
48 - 40 x } 


Z ee 
Tv,—Tu= oe (19 ene x} => 0. 


Insgesamt erhalten wir dann mit 


Jai! Bate pia ell tae = St Gy eS 
ee 5 a2) ang 60) Se 4g | 402 620 
in 0O<%<1 die Fehlerabschatzung 
441 , , 441 5 
= fe — AN ss = GO. 
|u(x) —v(x)|S CMTE und | «’(x) o'(*)| Seeazp% 


Speziell fiir x=1: |w(1) — 0,6958| <0,0037, | w’(1) — 1,0330| S$ 0,0223. 


Anmerkung. Uber Monotoniesaétze bei Anfangswertaufgaben hyperbolischer und 
iiber Fehlerabschatzungen aus den Defekten bei Randwertaufgaben elliptischer Dif- 
ferentialgleichungen soll in einer spateren Arbeit berichtet werden. 
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Stability and Boundedness of Systems 


TARO YOSHIZAWA 


Commumicated by L. CESARI 


Introduction 

Most problems in the theory of stability which have been discussed up to 
the present concern the behavior of the solutions relative to a fixed solution 
(cf. [2], [3], [5]), and most problems in the theory of boundedness concern the 
norms of solutions (cf. [6]). In this paper we discuss the same concepts relative 
to an arbitrary solution of the system and their relationships with the Liapunov 
function. 

Now we consider a system of differential equations 


(1) Pian 


In (1) x denotes an n-dimensional vector, and F(t, x) is a given vector field 
which is defined and continuous in the product space I x E’, where I is the 
interval 0Si<oo and E% is Euclidean u-space. Let x= x(t; %9, ty) be a solution 
of (4) through the initial point (¢), %)). We represent the Euclidean norm of % 
by ||| 
I. Definitions 

In this section we give the following definitions: 

(c) If corresponding to each «>0O and each 420 there exists a positive 
function B (fy, «) which for each «>0 is continuous in /) and is such that 


|| # (25 %. %) — % (4; Xo, t)|| SB (to, ) 
for all || x) —% || Sa and t=%), the system (1) is said to be equt-distance-bounded. 
(71) When B= (a) in (2) is independent of f), the system (4) is said to be 
uniform-distance-bounded. 


(271) If corresponding to each «> 0 and each f)= 0 there exist positive numbers 
Band T(t), «) such that 
|| ¥ (45 %0, fo) — %(¢5 Xo, fo)|| < B 
for all || % —%|| Sa and t>4+ T(t, «), the system (1) is said to be eqguiulti- 
mately distance-bounded for the bound B. 


(iv) When T= T(«) in (727) is independent of ¢), the system (1) is said to 
be uniform-ultimately distance-bounded for the bound B. 
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(v) If corresponding to each e>0 and each t,20 there exists a positive 
number 6(é), ¢) such that 


ll (C5: Monto HAE? Ne tg) 


for all || x9 —%o|| <6 and t2%, the system (1) is said to be egui-stable, where 6 
is continuous in fy. 

(vi) When 6=6(e) in (v) is independent of fo, the system (1) is said to be 
uniform-stable. 


(vii) If corresponding to each ¢>0 and each f= 0 there exist positive numbers 
6 (to) and T(t), €) such that 


|| ~(t; Xo, f) — x (t; Xo, te)|| <E 


for all || % —%||<5o (fo) and t>%+ T(t, €), the system (1) is said to be quasi- 
equiasymptotically stable. 

(viii) When every solution is stable in the sense of LIApuNov and the system 
(4) is equi-distance-bounded and moreover the system is quasi-equiasymptotically 
stable, the system (1) is said to be equiasymptotically stable. 

(ix) When 6, and T= Te) in (viz) are independent of f), the system (1) is 
said to be qguasi-uniform-asymptotically stable. 

(x) When the system (1) is uniform-stable and quasi-uniform-asymptotically 
stable, the system (1) is said to be uniform-asymptotically stable. 

To simplify the statements, we shall use ““D-bounded”’ in place of “‘distance- 
bounded’. 

We can also define the quasi-equiasymptotic stability of the system and 
the quasi-uniform-asymptotic stability of the system as follows: 

(xt) If for any e>0, any «>O0 and for any solution x=x(t; %, to) of (1) 
there exists a positive number T(t), e,«) such that if ||%)»—%9||Sa and 
t>ty+T(t, €,%), we have 


|| #3 Xo, to) — WUE; %, =e, 


the system (1) is said to be quasi-equiasymptotically stable. 

(x11) When T in (x7) depends only on ¢ and « but is independent of f), the 
system (1) is said to be guasi-uniform-asymptotically stable. 

Now we show that (x?) and (x77) are equivalent to (vii) and (ix), respectively. 

Lemma 1. Definition (vit) is equivalent to Definition (x1). 

Proof. It is clear that (x7) implies (vii), because it is sufficient that we 
consider some « as 6). We shall show that (vii) implies (a7). For any «>0, 
there exists a smallest integer N such that «/d9(¢))<.N. For any e«>0, we can 
find a positive number 7'(f, ¢/N) such that || x») — %p||< 6p (f) and >t) + T(ty, e/N) 
imply cal Cee to) — x(t; Xo, ty) ||<e/N. Then if ||%—%X||<a, by subdividing 
the line segment joining x) and X, we infer that || x(t; x9, to) — x(t; Xo, t) I|<e 
for t>%+ T(t, ¢/N). 

Lemma 2. Definition (ix) 1s equivalent to Definition (x12). 

Since the proof is the same as that of Lemma 1, it is omitted. 
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Remark 1. The above definitions of stability and D-boundedness imply that 
all the solutions of (1) are continuable to the right ; that is to say, all the solutions 
are defined for 4) <t<co. Conversely, when a solution is not continuable to 
the right, then the conditions in the definition are not fulfilled (cf. pp. 136-137 
in [6]). 

Remark 2. If we suppose that F(t, x) is continuous when 0<t<a, ||x||<co 
(aX oo), we can discuss the solutions when 0St<a. 


II. Relations 


In this section we shall discuss some relations between stability and D- 
boundedness. Of course, it is clear that the D-boundedness of a system does 
not imply the stability of the system (Ex. 8, p. 113 in [6}). 


Lemma 3. If the system (1) ts equi-stable, it is equi-D-bounded. 


Proof. By the equi-stability of the system, for some e>0, there exists a 
positive number 6(é)) such that if || *»—%||<d(f), we have for ‘=f, 


|| #5 %, to) — *(t5 %, b)|| <e. 


For any «>0, we can find the least integer N such that sy <N A henwit 
0 

is lear that if "|| 7px |[Soe “| 9 (GF Xo, %) — * (05%, %) |S Ne, where’ Ne is 

determined depending on « and ¢,. Therefore, the system is equi-D-bounded. 


Lemma 4. Jf the system (1) ts untform-stable, it 1s uniform-D-bounded. 


Proof. In the proof of Lemma 3, for any «>0, N is determined independent 
of tj, because we can choose 6 independent of f. 
To see the relations, we give an example. 


Example I. We consider a differential equation 


a2 (¢=0, x>0), 
\ diis2er stow pealldags4 4 SID ay ay tyity Vay v0) 
(2) dx 2 2isine t 
rE 38 (w#St, x>0), 
ae (0OSt<oo, <0). 


The general solution of (2) is x =C?sint+C(C=0) in the region 0OStSa, x20 
and is <= aC a (C=0) in the region w<t<oo, x20. Each solution of 
(2) is stable in the sense of LiapuNov, but the equation (2) is not equi-stable. 
Moreover the equation (2) is quasi-equiasymptotically stable, but it is not equi- 
D-bounded. 

From this, we can see that 

a) the stability of every solution in the sense of LiapuNov does not necessarily 
imply (v), 

b) (vit) does not necessarily imply (7), and 

c) (aii) does not necessarily imply (?). 
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The following example shows that the quasi-equiasymptotic stability of the 
system and the equi-D-boundedness of the system do not necessarily imply the 
equiasymptotic stability of the system. 


Example 2. Now we consider a differential equation 


2\x (OStS1,0<*SPh), 
2t (0S751,7=%), 
dx |(22—O){P—4t4+3+4V(P—414 3)? +44} (1<t<2,0<*S—P+4t—2), 
O) a= eae ely (1<t<2, —P+4t—2<x), 
Bs) (25t<00,0<%), 
x (0<t< 0, x<0). 


In the region 15/2, OS x —?#+ 4t—2, the solutions of this equation are 


a=—BP+ 4Bt+ (B?—3B) (B20). 


The equation (3) is quasi-equiasymptotically stable and is equi-D-bounded, but 
the solution x =0 is not stable in the sense of LIAPUNOV. 

However, in the case where the system (1) is quasi-equiasymptotically stable, 
for any solution x= x(t; %,¢) and any e>0, if ||%)»—%)||<a (a#>0: constant) 
and t>t)+ T(t), €), we have || «(¢; %, 4) — x(t; Xp, t)||<e. Then if every solution 
of the Cauchy-problem for (4) is unique in the interval 4; <?<#)+ T(t, €), there 
exists a positive number 6(<«) such that if || x —%9|| <6, we have || x(t; x», f) — 
x(t; Xp, t)||<e. Therefore if || x) — %||<6, we have || x(t; x9, 4) — x(t; %, t) || 
for t=ty); that is to say, every solution of (1) is stable in the sense of LIAPUNOV. 


Lemma 5. If every solution of the Cauchy-problem for (1) is unique, the quasi- 
equiasymptotic stability of the system implies the stability of every solution in the 
sense of LiapuNov. Therefore, in this case, if the system (1) is equi-D-bounded 
and 1s quasi-equiasymptotically stable, it is equiasymptotically stable. 

Example 2 shows the situation also for equiasymptotic stability in the 
large [6]. 

Example 1 shows also that 

a) (¢x) does not imply (x) (J in Example 4 is independent of f), 

b) (7x) does not imply (7) and 

c) (2v) does not imply (7). 

The following proposition is easily verified. 

Proposition 1. (a) (viii) implies the equiasymptotic stability in the large of 
every solution. 

(b) (~) implies the uniform-asymptotic stability in the large of every solution. 

From Example 4, it is clear that the converses are not true. 


As was stated before, the equiultimate D-boundedness of the system does 
not in general imply the equi-D-boundedness of the system. However, we have 
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Proposition 2. When the system (1) is linear, or, more generally, when 
F(t, x) satisfies a Lipschitz condition 
FU 2) — FWRI SLO |x — 3), 
where L(¢) is continuous and *€ E%, XC E%, then (ii) implies (7). 


Proof. By (7), for any «>0 and any solution x= x(t; %», tf), there exists a 
T(ty,) such that if ||%»—*x)||[Sa and t>f,+ T(t, a), we have || x(t; x9, f) — 
x(t; Xo, t)||<B, where B is a fixed constant. In the interval f,S¢<t,+T (ty, a) 
by the Lipschitz condition, 


4 


y+ T (toy) \ 
I ¥ (25 %o, to) — ¥(E; Xp, to) || S || % — %oll exp | fiapaaie £163) ds| 
to / 
tot T(to, a) 
& aexp| fe LAs) ds). 
Therefore, if 1=ty, é 


tot T (to, o) 
Wie (ts gs tol — a(t; Xt) || S max [2, exp | aga) as)) i 
fo 
This means that (777) implies (7). 


In the case where (1) is a linear system 


d 
(4) AT =AG x til, 
by the transformation 
(5) K=Y+ X(t; Xp, bo) 
the equation (4) is transformed into the equation 
dx 
(6) Ge =A (t) VE 


where x= x(t; %9,%) is a solution of (4). Therefore the behavior of solutions 
of (4) relative to the solution x(t; %, %) is the same as that of solutions of (6) 
relative to the solution y=0. Moreover, for any solution of (4), we obtain the 
same equation (6). For example, if a solution x= x(t; %9, %)) of (4) is stable in 
the sense of LIAPUNOV, then corresponding to each e>0 there is a positive number 
6 uch that if ||») —%9||<6, we have || x(t; x», to) —*(t; Xp, t)||<e for all th. 
Hence if || ¥o|| <6, we have || v(t; vo, f&)||<¢ for all ¢2¢). Therefore, the solution 
y(t) =0 of (6) is stable in the sense of LiapuNnovy, and then, from the stability 
of y=0, we can also see that every solution of (4) is stable in the sense of 
Lrapunov and that the system (4) is equi-stable, because we can suppose that 
ty=0. Therefore, we have 

Proposition 3. (a) The stability of y=0 of (6) in the sense of Lrapunov, 
the stability of one solution of (4) in the sense of LiapuNov and the equi-stability 
of the system (4) are equivalent. 

(b) The uniform-stability of y =0 of (6), the uniform-stability of one solution 
of (4) and the uniform-stability of the system (4) are equivalent. 

Corresponding to the D-boundedness of the system, when we discuss the 
behavior of solutions relative to some fixed solution, we shall say that the solution 
is D-bounded. 
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We can infer the following proposition in the same way as Proposition 3. 

Proposition 4. (a) The boundedness of (6), the D-boundedness of one solution 
of (4) and the equi-D-boundedness of the system (4) are equivalent. 

(b) The uniform-boundedness of (6), the uniform-D-boundedness of one 
solution of (4) and the uniform-D-boundedness of the system (4) are equivalent. 

Since the stability of (6) and the uniform-stability of (6) are equivalent to 
the equi-boundedness of (6) and the uniform-boundedness of (6), respectively 
(cf. p. 103 in [6]), we have 

Proposition 5. In a linear system, the stability in the sense of LiaPuNov, 
the equi-stability of the system, the D-boundedness and the equi-D-boundedness 
of the system are equivalent. Moreover, the uniform-stability, the uniform- 
stability of the system, the uniform-D-boundedness and the uniform-D-bounded- 
ness of the system are equivalent. 


Example 8. Now we consider a linear equation, 
(7) &* = (6teost — 21) x. 
The general solution is x= x,)exp(— 6¢psin ty — 6cost)+ t+ 6tsint+ 6cost — #). 
If we put —6fsint, — 6cost)+t=c (fy), we have |x| <|x%9| e°*°). From this, 
we can see that the solution x(t) =0 is stable in the sense of L1apunov. There- 
fore, the equation (7) is equi-D-bounded by Proposition 5. However, since 


«(2na+ os LONAOALL m)) 2° with u, the solution x(t) =0 1s not uniformly stable 


and hence the equation is not uniform-stable. This example shows that the 
equi-stability of the system and the uniform-stability of the system are different 
concepts. Of course, the equation (7) is equi-D-bounded, but it is not uniform- 
bounded. 

In the same way as the above we can infer the following proposition. 


Proposition 6. In a linear system, the asymptotic stability in the sense of 
LIAPUNOV, the equiasymptotic stability of the system, the ultimate D-bounded- 
ness and the equiultimate D-boundedness of the system are equivalent. More- 
over, the quasi-uniform-asymptotic stability (cf. [Z]), the quasi-uniform-asymp- 
totic stability of the system, the uniform-ultimate D-boundedness and the 
uniform-ultimate D-boundedness of the system are equivalent. The uniform- 


asymptotic stability is equivalent to the uniform-asymptotic stability of the 
system. 


Remark. In a linear system, the equiasymptotic stability of the system and 
the equiasymptotic stability in the large are equivalent, and also the uniform- 
asymptotic stability of the system and the uniform-asymptotic stability in the 


large are equivalent. 
III. Relations with the Liapunov function 


In this section we discuss the use of the Liapunov function V(t, x, y) in 


the present context. For this purpose, we consider the Liapunov function 
V(t, x,y) with respect to the system 


an Sa ‘ dy 
8 Gari), ares, 
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where the space considered is the product space I x E%xE%. Corresponding to 
the Liapunov function V(t, x, y), we define the function 


V'(t, x, 9) = lim {V(t +h, x +hF(t, x), y+hF(t, y)) — Vit. x, 9)}. 


Then we have theorems analogous to those for stability and boundedness [4], 
[5], [6]. 

Theorem 1. We suppose that there exists a Liapunov function V(t, x, y) defined 
in 0St<os, ||x—y||SH (H>0) in the product space I xE} XE’ and satisfying 
the following conditions; 

1° b(\|~—y||) SVE, x, y) Salt, ||x—y]|l), where a(t,r) and b(r) are positive 
for r>0 and are continuous increasing functions, and a(t,7)—>0 as r->0, 

2° V(t, x, y)€Co(x, y) (cf. p. 100 in [6]), 

Bo tex py) Ox 

Then the system (1) 1s equi-stable. 


Proof. For any e>0, if || ~ — y||=e, we have b(e)< V(t, x, y). Since a(t, 7) +0 
as r—>0, there exists a positive number d(é), «) such that if y<6d(t, ¢), then 
a(ty,7) <b(e). 

For some solution x= x(t; %, t)), we suppose that some solution «= x(t; Xp, ty) 
such that || %—%9||<d(%, €) satisfies || %(4; x, 4%) — %(t: Xp, %)||—e at some t, 
say t,. Then using the conditions 2° and 3°, we obtain an inequality 
b(e) SV (G,. x (4; Xs bo), (4; %o>%0)) S Vitor Xo. %o) S A(t, 0) < O(c). 
This statement is a contradiction. It follows that if ||%»—%||<6(t, «) and 
t=t), we have 
|| ¥ (45 %0, bo) — % (#5 Xo, b)|| < e. 

Theorem 2. If a(t,rv) in Theorem 1 1s independent of t, the system (1) is 
uniform-stable. 

Theorem 3. We suppose that there exists a Liapunov function V(t, x, y) 
defined in the product space IX EXE} and satisfying the following conditions; 

1° b(||x—y||) SV, x, y) Sa(t,||~—y||), where a(t,r) and b(r) are non-negative 
continuous functions and b(7r)—>oo as r—>oo, 

2° Vii, x, y) ECo (x, y), 

BY VG; %, V) SO. 

Then the system (1) ts equi-D-bounded. 

Theorem 4. If V(t, x, y) in Theorem 3 is defined when 0St<o, ||x—y||= Ro 
(Ry>0 may be sufficiently large) in IX EXE}, and if a(i,r) 1s independent of t, 
the system (1) 1s uniform-D-bounded. 

Next we discuss the equiultimate D-boundedness of the system and the 
equiasymptotic stability of the system. 


Theorem 5. We assume that F(t, x)€Co(x). In order that the system (1) be 
equi-D-bounded and equiultimately D-bownded, it 1s necessary and sufficient that 
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there exist a positive number B and a non-negative Liapunov function Vit, x, y) 
satisfying the following conditions in I x EZ XE}; 

1° b(||x—y||) SV, y) for ||x¥—y||ZB, where b(7) is a continuous function 
which is positive and increasing when r>B, and where b(r)—>0o as roo, 

2° Vit, x, y) sali, |l~>—y ), where a(t,r) is a continuous increasing function, 

3° Vit, x, y)ECo(t, x, y) and V'(t, x, y)S—cVit, x, y), where c is a positwe 
constant. 

Proof. The proof of sufficiency is the same as that of Theorem 11 in [6]. 
Now we shall show that the condition is necessary. Let 2, be the domain 


ostsy, |*llSa |lyllSa- 


If (f, %, Yo) €2,, we have ||x%)— yq||S2y. There exists a T(y) such that if 
t>T(n), we have || x(t; %, to) —¥ (ts Yo 40)||< 8, where B is the bound for the 
equi-D-boundedness of the system: for a solution z=2z(t; 29, t)) of the system 
(8), where z=(x, y), is continuable to t= and hence ||z(7; Zo, 4)|| <7, whence 
we have ||x||<7%, ||v||S7, and therefore it is sufficient that we take 7+7(27), 
where 7(27) is one for t=y, ||*—y||<27. We can assume that 7(y) is a non- 
negative continuous function defined for 720 and is monotone increasing in 7. 


For (f, %, Yo) €2, and 4) <¢<T(n), there is B(y) such that 
sup {|| «(¢; %q,t)||, Ily (2: Yo» ty) ||} S B(n), 
where #(7) is a positive continuous increasing function. Moreover, if we put 


x“ 4 = 24 + 
max ||F(¢, x)|| =2*(y), 
0S||x||<A(n) 


F*(y) is continuous and increasing, and we may assume that F'*(7)=1. More- 
over, let D, be the set of points (¢, x) such that 


OStST(y), 0<||x||<A(m). 
Since F(t, x) €Co(x), there is a positive constant L(y) such that 
|FG *) — F(t x) |] SL) lle — *', 


where (¢, x)€D, and (t, x’)€D,. We can also assume that L(y) is a positive 
continuous function and that it is increasing in y. 

Now we consider a function G(¢) as follows: 
Gam Bay - (Ceee) 


cO=| 0) (OSB )\e 


Moreover we put 

(9) 2F*(m) exp {(L(m) +1) T(m)} + 6 (28 (y)) exp {7() +} = Mn). 
Now we define V(¢, x, y) as follows: 

(10) VG, 2, y) = sup [G(llxGé+ 73 x, t) —y(t+7; y, t)\|)e"; OST]. 

It is clear that 

(11) G(||x — y||) < Vie, x, y). 
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Since the system (1) is equi-D-bounded and equiultimately D-bounded, there 


exist two positive constants f(t,«) and T(t,«) such that if ||~—y||<«, we have 
for T=0 


|x(t-+ 7; x,t) —y(t+T; y, t)||<B(t,«), 
and for t> T(E, «) 


llaté+ 7c; 4,2) —y(i+7; y,@)|| <B. 
Therefore 


2?) Vit, x, 9) SE (B(x — yl) eMle—M), 


and we can find a continuous function a(¢,7) such that V(t, x, y)<a(¢, ||*—y||). 
In the same way as in Theorem 11 in [6], we see that V(t, x, y) CC)(t, x, v) 


For instance, supposing that (?, v, y) €Q, and (t, x’, y’)€Q, and that V(t, x, y) = 
G(||~(¢+ 7; x,t) — y(t+7; ,8)|]) e&, we have 
Vit, x, 9) =Vit, x’, y’) SG(||x@+ 7; x4) —yt+; y,b||)e — 
—G(||x¢+r;«',t) —y(t+r: 9’, d||) 
ss oe {|| x ( (¢-+ 7; x,t) —x(t+ 7; x’,¢)|| + 
+\ly@+s ¥,8) —y@+7; yO} 
Zl el aL age 
etn) f 


In a similar way we obtain 


V(t, x, y) — Vit, x’, y') = — ek +T 


Therefore, by using (9), we have 

V(t, x,y) — Vee», 9") SMU) {lx — 2'l + ly — 9 I}. 
Finally, for two points (¢, x, y) and (¢’, x’, y’) which belong to 2,,, we have 
(13) Vt, x,y) — VU", y's M(m) [le (). 


/ 


 |le = 2] + [ly — 9 
Moreover, it is easy to show that 
(14) V(t, x, 9) S—VUex,9). 
Thus the theorem is proved. 

When the system (1) is uniform-D-bounded and uniform-ultimately D- 
bounded, there exists a B (a) such that if ||~—||Sa, we have for r20 


leé+rs x8) —yt+7; ¥, MS B(@) 
and in addition there exists a T(x) such that if ||x—y||<a, we have 


|n(¢+7t5 4,2 —y¢+t3 y,b||<B 
for t>T (a). Therefore 


V(t, x, 9) <G (B(x — yl) eT); 


that is, there exists a continuous function a(7) such that V(t, x, y)Sa({||x—y]]). 
Thus we have the following theorem: 
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Theorem 6. We assume that F(t, x)€Co(x). In order that the system (1) be 
uniform-D-bounded and uniform-ultimately D-bounded, it is necessary and sufficient 
that there exist a positive Liapunov function V(t, x, ¥) satisfying the following ) 
conditions in 0Xt<o, ||x—y||2Ry (Ro may be sufficiently large) : 

1° b(||x—y||) SV(t, x, vy) Sa (||~—y||), where a(r) and b(r) are positive con- 
tinuous increasing functions and b(r) +00 as r-> 00, 

2° V(t, , 9) CColts % 9), 

3° V(t, x, y)S—c(||x—Y||), where c(r) is a positive continuous function. 

Remark 1. In Theorems 5 and 6, for the sufficiency of the conditions we need 
not require that F(t, x) € Cy(x). 

Remark 2. When F(t, x) is autonomous or periodic in ¢, we can obtain a 
V(t, x, vy) which is autonomous or periodic in ¢. This is clear by the definition 
of V(t, x, y). 

By using Theorem 5, we can obtain the condition for the equiasymptotic 
stability of the system. 

Theorem 7. We assume that F(t, x)€Cy(x). In order that the system (1) be 
equiasymptotically stable, it 1s necessary and sufficient that there exists a Liapunov 
function V(t, x, y) satisfying the following conditions in I XE. XE}; 

1° V(t, x, y)=0 if ||x—y||=0 and Vit, x, y)>0 if ||x—y||+0, 

2° b(|\|x—y||) SVU x, y) Salt, ||” —y||), where a(t, r) and b(r) are continuous 
and b(r) is an increasing function such that b(r)>0 for r>0 and b(r) ce with r, 

3° Vit, x, y)ECy(t, x, y) and V't, x, y)S—cVit, x, y), where c 1s a positive 
constant. 


Proof. For the sufficiency of the conditions, we need not require F(t, x) € C9 (x). 
If we assume that I(t, x) €C)(x), every solution of the Cauchy problem for (1) 
is unique, and hence, by Lemma 5, every solution is stable in the sense of 
Liapunov. However, we can verify the stability without this assumption. Since 
V(t, x, y) is continuous and vanishes at x=y, for any e>0 and for fixed %, 
%q there is a positive constant 6 such that if || x) — v|| <6, we have V(ty, x», V9) < 
b(e). From this, by using the conditions 2° and 3°, we can see that if || %» — yo|| <0, 
|| ¥ (£5 X09, %) —¥(t3 Yo. to)||<e for tt). Therefore, x= x(t; xp, fo) is stable in the 
sense of LIAPUNOV. 

Since B in Theorem 5 is chosen arbitrarily, it is clear that the condition is 
sufficient. 

Now we show that the condition is necessary. Since the system is equi- 
ultimately D-bounded for any positive number ¢, we write T(y) in Theorem 5 
as Te, y), and we assume that T(e, 7) =T(1, y) for e>1. Then M(y) in Theorem 5 
depends also on ¢, and so we write it as M (e, 4); that is, 


M(e,n) = 2F*(e,n) exp{(L(e,y) + 1) T(e, n)} +G(e, 2B(e, n)) exp{T(e, y) + ny. 


M(e, 7) is a positive continuous function defined for e>0, 7=0. By MAssERA’s 
lemma in [4] (cf. p. 133 in [6]), there exist two functions g(e), h(n)EC, in 
[0, +00) such that h(y)>0, g(e)>0 when e>0, g(0)=0 and that 


(15) g(e) M(e,n) Sh(n). 
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To simplify the statements, we denote G(e,¢), B(e,7), T(e,7) and g(e) for 
e=1/n (n=1, 2,...) as G,(C), B,(n), (4) and g, respectively. 
Now we put 
(16) Vi (4%, ») = 8, sup [G, (llx(@+ 7; x 1) — y(E+ 73 9, 8)])) &; OS 7). 
It is clear that 


G, (lz — yl) SVG», 9) 
and 


V Atopy s Ol alii ws 
If, (é, x, y)€Q,, 


(17) V(t, %, 9) SB, G, (2B, (n)) e™ < h(y) (by (15). 
On the other hand, by the equi-D-boundedness of the system and the equi- 


ultimate D-boundedness of the system, there exist two functions B (t, yl) 
and 7, (t, ||~—y||) such that 


(18) Vis(b %, 9) S8n Gy (Bt le — yil)) ef le), 
If (¢, *, y)€Q, and_(’’, x, y')€Q,, 


V(t, x 9) —V, (0, xy’) Sa) [[t — | + || 2 
Moreover, 


y'|{I- 
Vs (t, x; y) SS V,(, x; y)- 


We are now going to obtain the desired function. We put 


(19) Vix, 9) = =e V,, (t, %, 9) 


Then V(é, x, y) is defined and continuous in the product space I x EY x E}, and 
if x=y, V(t, x, y)=0. In the same way as in Theorem 21 in [6], we can show 
that there exists a continuous function b(7) such that b(||«—v||)<V(é, x, y), 
b(r)>0 for r>0 and b(7) co with r. From (18) we have 


a ix a3 “l) ene nea 


On the other hand, if we put max (f, 
following inequality holds: 
T,, (%, 


elnltlle—vil) < eT) < M,, (1). 


(20) 


y||) =, we may consider that the 


Hence 


Therefore 
x, elnlt Ileal) < h(n) 
<h(n(tllx — yll)), 


whence by aid of (20) we can see that there exists a function a(t, 7) such that 


V,,(t, %, ¥) Salt, iz 
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Hence & 
eS 1 
V(t, x, y) = ar Vex») Salt, lle—yl) 2.3 
n=1 n= 
<a(t,| *—y|l). 


As in Theorem 21 in [6], we can see that 
Vit, x, v) E Colt, %, ¥) and V'(t,x,¥)s=—V(t,*,y). 


Next we shall discuss the uniform-asymptotic stability of the system. In 
this case, the system (1) is uniform-stable, and so by Lemma 4 it is uniform- 
D-bounded. Therefore there exist two continuous functions «(7) and B(7) such 
that for T20 


a(\|~ — yl) S|la@tr: x2) —y¢4+7; y, || SB (ll* — yl) 


and f(r) +0 as r-+0. Since the system is uniform-asymptotically stable, T is 
determined as a function of ||x—y|| only. Therefore, in place of (20) we have 


ie (t, x, y) = Sn & (B (Ix aos yll)) etallle— sil), 


Since J, is independent of ¢, we may assume that for 0St<oe 


elle) <M, (IJ — 9). 
Hence 
Va (t, %, 9) S Gq (B (ll* — yl)) Sa erm 


<B(\|x — y|l) g,.™@, (lx — yI)) 
SB (|x — y||) A(I|~ — 9) 


Therefore we obtain the following theorem. 


Theorem 8. We assume that F(t, x)€Co(x). In order that the system (1) be 
uniform-asymptotically stable, it is necessary and sufficient that there exists a 
Liapunov function satisfying the following conditions in IX EY XE}; 

1° V(t, x, y)=0 if ||x—y||=0 and Vie, x, y)>0 tf ||x—y|| +0, 

2° b(||x—y||)SV(t, x, y)Sa(||x—y||), where a(r) and b(r) are continuous in 
r and b(r) 1s an increasing function such that b(r)>0 for r>0 and b(r) >co with 
vy and that a(r) +0 as r->0, 

3° Vit, x, y) EC Colt, x,y) and V'(t, x, y)S—ce(||x—y 
tinuous and positive definite. 


), where c(r) ts con- 


Remark 1. When F(t, x) is autonomous or periodic in ¢, we can obtain a V 
which is autonomous or periodic in f. 


Remark 2. By Definition (x), it is sufficient that V(¢, x, y) in Theorem 8 be 
defined when 0<t<oe, ||x—y||<H, where H is some positive constant. 
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Fehlerabschatzungen 
fir die Ableitungen bei Randwertaufgaben 
mit linearer elliptischer Differentialgleichung 
A. G. MEYER 


Vorgelegt von L. COLLATZ 


Wie man Schranken fiir die Ableitungen der Lésung einer Randwertaufgabe 
aufstellen kann, wenn fiir die Lésung selbst bereits eine Abschatzung vorliegt, 
ist fiir die Potentialgleichung schon von O. D. KEettoce ([8], S. 227), fiir all- 
gemeinere elliptische lineare Differentialgleichungen von J. SCHAUDER [10], [11], 
fiir Systeme von elliptischen Differentialgleichungen von A. DouGLis und 
L. NIRENBERG [4] gezeigt worden. 

Dem Gedankengang der Arbeit [/0] von J. SCHAUDER folgend, werden in 
Nr. 1 Abschatzungsformeln mit zahlenmaBig bestimmten Konstanten angegeben, 
die im Inneren des betrachteten Gebietes gelten. In Spezialfallen lassen sich 
die Schranken auch auf den Rand ausdehnen. 

Die mathematische Behandlung von Strémungen fiihrt oft auf reine zweite 
Randwertaufgaben, bei denen die Konstante Lésung der homogenen Aufgabe 
ist. Von besonderem Interesse sind dann meist die Tangentialableitungen auf 
dem Rande. Fiir spezielle solche Aufgaben in der Ebene wird in Nr. 2 gezeigt, 
wie man Schranken fiir die Tangentialableitung gewinnen kann. 


Nr. 3 bringt einige Beispiele. 


1. Fehlerabschatzungen im Inneren des Grundgebietes 
1.1, Differentialgleichungen mit konstanten Koeffizienten 
Mit G sei ein offenes, zusammenhangendes, beschranktes Gebiet im 4%, ..., %,- 
Raum bezeichnet, mit J" der Rand von G. Fiir das n-tupel f,,..., £, wird auch 
kurz p oder P geschrieben (gelegentlich » als Spaltenvektor aufgefaBt). 
Der Abstand eines Punktes P vom Rande I" sei mit R(P) bezeichnet, d-h. 
Ba Pee pea |S e- qi)”. 


R sei der Radius der groéBten Hyperkugel in G+, d.h. R= sup R (P). Zunachst 


Noe, n 2 PEG 
sei mit dem Laplaceschen Operator A = a und konstantem c= 0 der spezielle 
lineare Differentialausdruck jai 9%) 
(1.1) L, (uw) = —Au+cu 


betrachtet. 
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Satz 1. Sez u(P) eine in G+T stetige Funktion, die in G stetige partielle Ab- 
leitungen bis zur zweiten Ordnung einschlieBlich hat. Sei L,{uj)—F(P) in G mit 
einer Funktion F(P), die in G stetig und beschrinkt ist und dort stetige partielle 
Ableitungen erster Ordnung hat. Sei Py ein Punkt aus G und fiir r< R(P,) 

Q, die Menge der Punkte Q mit | Q—P,| =r, 
K, die Menge der Punkte Q mit |Q—P,|<r, 
K, die Menge der Punkte Q mit |Q—P,|=r. 


Dann gilt 

(1.2) |as(Po)| S Ay-7- Max |F(Q)| + dae» Max | (Q)], 

(13) |ue(P)| SA 7 Max| PONTE) 4. 4... Max Ju (Q)], 
Hey (Pe)| S (|cos| By, + [sino] Bys) -7- Max Se ae se 

(1.4) + (|cos| Bo: +|sine| By») -— Max | (Q)| + 


2 eee Pie w eh). 


n 
Daber sind & und n beliebig vorgegebene Einhettsvektoren, uz, uz, die Ableitungen 
von u in den entsprechenden Richtungen, cosw = (&,n) = >) &; n; das innere Produkt 
j=1 
von € und yn. Die Konstanten A;, Af, B;; hingen noch von der Dimension n ab 
und sind in der Tabelle I angegeben. 


Bemerkungen zu Satz 1. Die Abschatzungen vereinfachen sich etwas, wenn 
man die auftretenden Maxima durch die entsprechenden oberen Grenzen iiber 
grad F(Q)| fiir ara setzt und r= R(P,) wahlt. AuBer- 

=e 


dem kann in (1.4) gesetzt werden: 


ganz G ersetzt, 


(|cosw| By, +|sinw| By.) =By= Viste Bis 
und entsprechend 


(|cos | Bo ae | sina | Bo») = By = VB3i + Bie. 


In Tabelle 2 sind fiir die Koeffizienten von (1.2) bis (1.4) im Falle c>0, 
m=2 und »=3 Formeln angegeben, die auBer von m auch noch von c und 7 
abhangen und insbesondere fiir groBe c oder groBe 7 bessere Schranken liefern. 
Fiir die Funktion uw selbst werden Schranken als bekannt vorausgesetzt. Diese 
seien etwa durch Anwendung der Methoden von [9] gewonnen. 

Trotz des Faktors 1/7 beim zweiten Glied der Formeln (1.2) und (1.3) laBt 
sich in gewissen Fallen auch noch am Rande eine Abschatzung mit diesen 
Formeln durchfiihren, nadmlich immer dann, wenn fiir «(Q) eine Abschatzung 
vorliegt, die geniigend stark verschwindet fiir R(Q)—0 (s. Beispiel 3.2). 


Beweis von Satz 1. P, ist ein beliebiger, fiir den Beweis fester Punkt aus G. 
Einfachheitshalber sei der Koordinatenanfangspunkt nach P, gelegt, also 
P,=(0, ..., 0). AuBerdem sei das Koordinatensystem so gedreht, daB der Vektor & 
in die Richtung der positiven x,-Achse weist und 7 in der x, x,-Ebene liegt, 
ile == (1, 07 s.., G); qalCos@; Sinw,;.0, <4.,.0). 
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Sei é fiir n=2, a fir n=3 
ai 2 
a = DA wane) Ss Wee 
ete | Pane ea a Yea -— fiir n>2 gerade 
Ix (n) [Jin o ¢ x) cu ae & 
Boo ag 2) A y a 
STR Gp laete fiir m >3 ungerade 
Tabelle 1 
n=2 N=3 allgemein 
4 3 an : 
cl igh if: 
Se ee ge (n—1) (ny 
4 2 2n : 
: 3 ae K 
A a 2 (n—1) a 
1 3 n £ 
Als || 5 — K( 
aE Dice 2000 (m—1) (n +2) =e) 
8 103 [m—iy=! nt2 . 
B os ae | Bde ARS de CLF 
sot Of a | n” nee wee 
Pe eR pee m+2 2 
52 Byes | nmt+3n 
eee Sk - 
8 20 3 n—1)"*-1 
Bag'| 5 \ a (n+2) 4K(n) 
8 10 ris eZ 
Be Bly {2} 
22 2 = (m+ 2) =< 
Tabelle 2 (mit e=Yc-7) 
n=2 n=3 
a ee a EN i ee ee | Pe 
@ 
1 “By Le Tey ete e).\ glo: . Sue )— (cosh gp—1— $e?) | 
PA {-§# té)+6& f) ee dé) =. oe ea ag 
sar AH E)FEATS) Ave) ie 8 30 (9 cosh 9 —sinh 0) 
0 
ie dirs Bee < ci P iaitsoveatal ae 
Se 21 (iQ) D 3 (g cosh g—sinh 0) 

At hee 40+8if(iQ) ee 10[9 cosh o—sinh 0 —4 0°) 

Se 20 072 J, (2 @) 20 Q (g cosh o—sinho) — 

; — < a 

B et : 373 
“a ; : ie (s. u.) melee sae n(Sae) 

Br — hie A Ary 2ONB <a Nena ae Soe 

8| Joli 0) PEAT G 0), 9 15{o?sinh e—3ecosh o+3sinh ro 
Q 

By» = 53," (8:'0,) als * By, (s.0.) 

a ee 
e 2 ; 

a neal pl BO bb 5) ; ok ea 1 H(t 9) —— A, (ig) 
maa: Brim oe | eyiMeis)—2 Mie) —[f8) +26 509]—=——2& lag 
: 0 Jolie) + ot Alee) 
n=3: Byj= _ : 1602 cosh g—4 489 sinh 9+48 cosh p— 2 0! +80?—48 


5 e® sinh g—15 0? cosh +15 e@sinh @ 
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Zu L, und K, existiert die Greensche Funktion G(P, Q), es ist in K, 


ay Fle P,Q)-4dVo fag u(Q)-Gre(P, Q)-dSo, 


wo dV, das eee in K,, % die auBere Normale und dS, das Ober- 
flachenelement auf Q, ist. Mit Beriicksichtigung des Verhaltens fiir P=Q 
(E. Hopr [7], §2) erhalt man durch Differentiation: 


(1.2*) 2 ee ) + Gey ( es Q) - dVo =af u(Q) ; Groth ke Q) -dS 
a u(P) = F(P)-| f G(P,Q) -4Vo]_ + 
gts bs i 
at [F(Q) —F(P)]-Ge,(P; Q)- dVo sate Gti (Py Q)- aSo, 
uz, (P) = F(P) | [olP, Q) aVo|. ue 
(1.4*) ~ tee 


=a [ F(Q) —F(P)| = Gast nlop Q)- )-dVo —fu(Q (Fo a kehoiie Q)- dSo. 


Die angehangten Indizes P, Q sollen angeben, beziiglich welcher Koordinaten 
die Ableitung zu nehmen ist. 


Nun ist g(P rails P, Q)-dVo BRE Serpe beziiglich P,, also 
9; (Po) = 9; x (Po) =0 fiir jk und 9,,(%)= Ag (P »); auBerdem gilt 


Ag —co—A, 
also 


| {G(R Q)-aVo), = "le SG, Q)-aVo— 
Ky 7] 


Mit GE, 4.(P, Q)'= Gepap(P, Q) — ——— ¢ G(P, Q) erbalt man 


tts (Pp) = { [F(Q) —F(P)]- Gbonp (Po, Q) dVo — 


= f 4(Q)- Gg tonp (Pos Q) : dSo+ — cu(L) =F (Pie 


Aus der Abschatzung der Darstellungen (1.2*), (1.3*) und (1.4**) ergeben sich 
fiir P= P, die Formeln (1.2) bis (1.4), wenn man die A,, Ay, B;; so wahlt, daB 
gilt: 
Bes {> iG Olav a= fr FG etre CO) as, 
Ky 2, 


—f£ 
AY A) ea |Gep(Po, Q)| -dVo, 
il MOE. 
|cosm| By, + |sinw| axe ee eo ehat paler: Q)|-4Vo, 


costo Bay + [sino] Baa J | Groepnp (Po, Q)| -dSq- 
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Zur Auswertung dieser Integrale werden zunachst die bendétigten Greenschen 
Funktionen und ihre Ableitungen aufgestellt. “ot 

Die sog. Hauptlésung y;,(@) hangt nur vom Abstand a=V>.(x;— 9;)? der 
Punkte P=(%,..., %,) und Q=(yy,---, Yn) ab: 


Sopa dbo ce? fiir n=2 | 
ae == 0 
Vg?) ni a fiir n> 2 | 
(Ga Dy oy oh) Ae = oP 
yi(a) =) gq Bore yea) Br eee 
: Prypa Neal seetsre ts Pa 
+(e aa 2 HD Gee) fiir n >2 gerade c= Ds 
4 22 — 
i ive a pe 1 <= . d 
ae (- +) 2 fh : At n=2 (‘jc a) fir n>2 ungerade 
Als regulare kugelsymmetrische Lésung bietet sich an: 
4 frees Or 


— Ne) 2 a Phere Mea) a. Greies O. 


Im folgenden wird auch bei der Greenschen Funktion immer noch die Di- 
mension als Index mit angegeben. Wenn der Index 1+ 2 auftritt, soll die ent- 


sprechende Funktion des R,,,. betrachtet werden fiir x 
Entsprechendes gilt fiir ~+ 4 und +6. Der Ansatz 


n+1— Xyn+e2 = Vat a =4,19=0- 


Gu (Py, Q) = Yolo) — 2g, (0) + (Ps Q) 


mit o0=|P,—Q| =\>\ y? fiihrt auf /(P,, Q)=0 in K, als eindeutige Lésung der 


Randwertaufgabe — Agh(Py, Q)+ch(Py, Q)=0 in K,, h(Py, Q)=0 auf Q,, so 
daB gilt 


G, (Py, Q) = (0) — ZH ¢¢ (9). 


Vnép (4) = — 20 (%,— %)- -Vn42(a), 


Vn Ep np (a) = COS@- < ; Vin (a) SF 


- 4n* cos ON Mts ae a t SIN @ + (%, — 4) (%2 — ys) -Yn+4 (a) 
erhalt man entsprechend 
Gnep(Ly, Q) = 22 yy (CP aEY WT BOA, 


9 2 A 
Greer eos Q) = 4*|cos GQ): (i ra ) + SIN @ - y;- ys| CES Q) 
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und durch Ableitung nach @ und 9-7 auf dem Rande Q,: 


‘ 0 fit cs= 0 
G,, nNQéP (Py, Q) = 47 Vi he (7) si | Vi (7) | fa 6 
7 ur ¢ F 


Gin inp Po Q) == 828 1 ‘cos : (vi _ an + sinw 4, ys| x 


Ww 


‘i 0 fir c=0 
x |¥n+6(7) + asalr) : 


Grey) |. file = 0. 
‘Sed 14 (7 y) 


Damit lassen sich alle Integranden explizit angeben. Ihre Integration ergibt 
fiir c=O die Tabelle 1. Fir c>0 ist die Integration teilweise nur angendhert 
méglich. In Tabelle 2 sind fiir ~=2 und n=3 die Ergebnisse, soweit elementar 
angebbar, zusammengestellt. Die Schranken werden kleiner bei wachsendem c. 
Einfacher kann man auch fiir c>0 durch die Formeln (1.2) bis (1.4) unter 
Benutzung von Tabelle 1 abschatzen, wie sich aus den Ungleichungen 


| Gai (By Q) Ls s has Pas Q)| 
é 
jp Cvri(Po QO), S |g, Crai(Fo. Q) 


ergibt, wobei der naecens Index c bzw. 0 angeben soll, daB die Greensche 
Funktion fiir den Differentialausdruck —AG-+-cG zu nehmen ist mit c>0 bzw. c=0. 

Die zweite Ungleichung ist eine direkte Folge der ersten, die erste folgt aus 
der Monotonie der Randwertaufgabe ({9], Monotoniesatz 1) mit dem Laplace- 
Operator, denn es ist 


— Ao Gisi(Py, Q) =— ¢Ga4i(Po, Q) SO= — Ap Gr4i(H,Q) in K,, 


G? (Py, O)=G), , (P,, O)=0 auf QO, und fir kleines 9 ist Gh, ;(P,, Q)=Giu;(Ay, Q), 
wie man leicht aus der Foi nentwieklane der Zylinderfunktionen [3], (Bd. 1, 
S. 420, 432) sieht. 

Damit ist Satz 1 bewiesen. 


i=2,4 
—< 


o=7 


Nun sei allgemeiner in G ein linearer Differentialausdruck 


(1.5) L [uw] = — D) ag) Mp1 — D0; 4, + cu 

kis al 
vorgelegt mit konstanten reellen Koeffizienten a,,, 0;, c. Bei der Funktion w 
geben die Indizes die betreffenden partiellen Ableitungen an. Der Differential- 
ausdruck erfiille die Voraussetzungen : 


(1.6) Die Matrix A=(a,,) ist symmetrisch und positiv definit. Es ist c20. 
Die Eigenwerte 2; von A (d.h. die Wurzeln von det(A —AE)= 0) sind also 


positiv, es gibt zwei Zahlen /, und J), so daB 
(1.7) Oral dak geist flit elon ty. 
Mit den Koeffizienten }; von (1.5) und J,, /) wird gebildet 


(1.8) es 


| VA. 
21, 1,’ l= |) Sp 
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Die Differentialgleichung L[{w]=/(x) fiir w(x) 1aBt sich reduzieren in eine Diffe- 
rentialgleichung L,{[v] =F(%) fiir eine Hilfsfunktion v (¥ ) ({3], Bd. 2, S. 137), auf 
welche die Abschatzungen nach Satz 1 angewandt werden kénnen: 

Sei T die orthogonale Matrix der Transformation von A auf Hauptachsen, 
D die eae ta mit den Elementen d,;—|/A;. Durch die affine Transforma- 


f(x) tiber in L[@] =f (%) mit L[@] = Aa — > 5,a,4 cf, 
j=1 


Im *-Raum wird die Hilfsfunktion v(% =A (0,2) % (z) 


i (3) = u(x), {R= (x 


“|= 
? 
eingefiihrt mit (b, %) = 2 b,x;. Aus E [wz] =f (%) wird dann die Differentialgleichung 


Lv) =—Av+ev=F(%) mit €=c+2/b|2, FH) =f (HZ) Abe”) 

Die Formeln (1.2) bis (1.4) von Satz 1 liefern Schranken fiir die Ableitungen 
von v(%). Riicktransformation auf die Funktion u(x) im x-Raum ergibt den 

Satz 2. Sei u(P) eine in G+T stetige Funktion, die in G noch stetige partielle 
Ableitungen bis zur zweiten Ordnung einschlieBlich hat. Set L{u\)=f(P) im G mit 
einer Funktion f{(P), die in G stetig und beschrankt ist und dort stetige partielle 
Ableitungen erster Ordnung hat. 

Dann gelten im Punkt P,€G die Abschatzungen: 


(1.9) | ws (Py) | a Aihantn Mex I/O) ae pe or Mae ake Q)| +0. 5° -|u(B) 


u u 


|w:(R)| S Afr ae a oe +a |f(P)|| + 


1(Q) —f (Lo) 
Me OE 


(1.10) ; ou . 
+ Se bE" Max |u(Q)] +a 4 |W(R) 
He (P| S (By + 20 AB) [Max | OMA) 4 a 17 (|| + 
ed) ne -} 2a <4) = Mae lel |w(Q)| + 


+ Se [lf (Pa)| + (e+ (+ 1) 22) | (B)I] 
fiir r= R(B). 

Daber sind € und y wieder beliebig vorgegebene Einheitsvektoren. A,, A*, B; 
sind die Konstanten der Tabelle 1 mit B,=\B?,+B?,, B B;=/) B24. B%.. 


Die anschlieBend an Satz 1 gemachten Bemerkungen kénnen sinngemaB hier- 
her tibertragen werden. 


1.2. Differentialgleichungen mit verdnderlichen K oeffizienten 
In G sei ein linearer Differentialausdruck 


Hee Dn 2 ay P\ ae =P (P) +.c(P) -u(P) 


j,k=1 
vorgelegt. Er erfiille die VoLustal ey 
(1.13) Die Koeffizienten ayn, b;, c sind stetig und mit stetigen Ableitungen erster 


Ordnung versehen in G-- J. Die Matrix A der a;,(P) ist symmetrisch 
und positiv definit in jedem Punkt PEG4+T und es ist c(P)=O dort. 
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Die Eigenwerte j;(P) der Matrix A(P) sind also positiv fiir PEG+I. Wegen 
Stetigkeit im abgeschlossenen Bereich G+" gibt es analog zu (1.7) zwei Zahlen 
1, und J), so daB gilt 


(1.14) 0<1,<V4,(P) <1, fir i=1,....n und PCGHLT. 
Analog zu (1.8) wird gebildet 


1 4 . W n is ar 
(1.15) a= 5) Max |b(P)| mit o(P)| =| Dye), 
Zu den Koeffizienten von (1.12) werden Konstanten c, A, B, C gebildet mit 
G2oMaxre(P); 
PEG 
A= reer piles o= etwa 2a Mer |gra ays), 
(1.16) as b 
. b;(Q)—); ) aad - << ah 
B= pilax.. 2, Se " [etwa ae Max grad bil), 
O<G ial = 
c(Q)—e(P) — Hie 
C2 Max | Aes | (etwa C = Max |gradc|). 
QeG4r 


Seien weiter U, i, 4, Konstanten mit 


Ui eax | 2 ( Py ioe | lye wee \E(P)|., 


PE€G+r 
(Q)—F(P) 
ne alps, || |=] 21 


QE€G4r 


(1.17) 


Mit diesen Konstanten und mit den in Tabelle1 angegebenen Ay AS AA 
B, =\Bi,+ Be., By= \/B2,+ B3. werden Hilfsfunktionen gebildet: 


hyo (8) = det AE (H+ CU) +a ace 
oe bn ai 
d2etRO AF. R? ek ee Res A 
hy, (0) == a ly, ip Be hy »(0) (1— 6)? Li, fe , 


Hle+ (m+1) 0? G ) Ul 
4n lz 


hy (0) = 6 e*®°(B,+2a Af RO) a (RaeC Uy ke? [£o- 


(1.18) e%ko (ks 


hy (0) = 2. e*®?(B, + 20 AERO) = 


f= 6: eB ‘ 


) * 
Neg (0) = CE) e*RS(B + 2a AT RO) aA. 
Alle diese Funktionen sind positiv und stetig in offenen Intervall 0<d0<1. Fiir 
6—>1 wachsen die Funktionen h,;(6) mit 7,7=1, 2 iiber alle Grenzen. Fiir 60 
on ( feet 2 


Aare 
lim A,;() =0 fur | fact 
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Zu vorgegebenem k>0 gibt es also ein e>0, so dab h;,,(0)<h ist fiir 0<d<e 
undiae/ ae dy 2: 
Damit kann der Abschatzungssatz formuliert werden: 


Satz 3. Sei u(P) eine in G+I" stetige Funktion mit stetigen partiellen Ab- 
leitungen bis zur zweiten Ordnung einschlieplich in G. Mit dem in (A 12) defimerten 
Differentialausdruck L sei L{w\=F(P) in G beschrankt, stetig und mit stetigen 
partiellen Ableitungen erster Ordnung versehen. Unter Benutzung der Bezeichnungen 
von (1.14) bis (1.18) seien 6, und 6, zwei Zahlen mit 0O<6,, 09<1 und 


hy (04) . hy 1 (02) 


— <1. 
1—hy,(0,)  1—hg2(0s) 


Iy1(0,) <1, Mtg (be) <1, 


Dann ist fiir beliebige Richtungen &, 4 


R(B,) - |tte(B,)| <a) ata eo __ 
(1.19) : (1 fy 1) (1 hye) —Iy 9 hey 


(1 —hy 1) hoo thei Mio 
2(P\ lu. a : 
B (Fo) IMs (Fa)| ~~ (1—fyy) (1 he2) Me her 
Dabei ist P\EG beliebig. Fiiri=1, 2;7=0, 1, 2 1st kurz h;;(0;) =h,; gesetzt worden. 
Fiir eine méglichst gute Abschaétzung kann es niitzlich sein, als G nicht das 
Gebiet der vorliegenden Randwertaufgabe, sondern das Innere einer Hyperkugel 
um P, in diesem Gebiet zu nehmen. 


Beweis von Satz 3. P, ist ein beliebiger, fiir den Beweis fester Punkt von G. 
G* sei das abgeschlossene Teilgebiet |P—P,|SdR(Pj) mit 0<d<1 von G, 
R*(P) der Abstand des Punktes P vom Rande von G*. Die Ableitungen von u 
sind in G* beschrankt. Sei 


2) N* = Max R* mall? 1k Spas x R*2 lic 
(1.20) A Noe GRY 2G Peal en ek (P) -|%,,.,(P)|, 
lor|=1 \Q2|=1 
|Q3|=1 
wobei zu beachten ist, daB die GroBtwerte auch iiber alle méglichen Richtungen 
(Einheitsvektoren) 0,, 02, @3 Zu nehmen sind. 
Damit gibt es zu y>1 im Inneren von G* Punkte PB und P, sowie Rich- 
tungen 01, 02, 03, So daB gilt 


Nit <y--R*(R)-|u,,(B)| und NF Sy-R*2(B)-|u,, 


Zu den F, 1=1, 2 werden Differentialgleichungen L*{w]—/,;(P) mit konstanten 
Koeffizienten gebildet: 


Mit ihrer Hilfe lassen sich die rechten Seiten obiger Ungleichungen weiter ab- 
schatzen tiber die Formeln (1.10) bzw. (1.11). Unter Benutzung der Konstanten 
von (1.16) und (1.17) und der Schranken N* fiir die Ableitungen von w erhalt 
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man fiir y—1 die Ungleichungen 
Ny S hy 9 (01) rip hy, (0;) ; Ne ai hy 2 (04) : Ny 
NF = hy 9 (de) 3 hy ; (02) ; Nv at hy » (0s) i N+ 


mit den Funktionen /;; von (1.18) und 7,=6; R*(B), 0<6;<1. Nach Voraus- 
setzung wurden 6, und 6, so gewahlt, daB diese Ungleichungen aufgelést werden 
k6énnen, man erhalt 


Nes (1 —hz2) Moth, hyo ; Ne a (1 —hy 1) hag the hyo 


Pe( hy) (1 hy») hain : (1 hy) (1 ho») bale 
Die rechten Seiten sind unabhangig von d. d—>1 liefert die Behauptung des 
Satzes. 


Aus Satz 3 folgt insbesondere, daB R(P)-|u;(P)| und R?(P) -|u-,,(P)| be- 
schrankt sind in G+ TI. 


2. Fehlerabschatzungen fiir die Tangentialableitung 
auf dem Rande des Gebietes bei ebenen Aufgaben 
Sei G ein offenes, beschranktes, einfach zusammenhangendes Gebiet in der 
x, y-Ebene mit stiickweise glattem Rand J”. In G sei der selbstadjungierte Diffe- 
rentialausdruck 


(2-4) Ewe (abu) a (bp cu,), 
vorgelegt mit Funktionen a(x, y), b(x, y) und c(x, y), die inG+T/'stetig sind und 


dort stetige partielle Ableitungen erster Ordnung haben. Die Matrix A -(; , 
sei positiv definit in G+ J’, d.h. es sei a>0 und D=ac—l?>0 in G+I. G 
Von den Funktionen, auf die der Differentialausdruck (2.1) angewandt wird, 
sei vorausgesetzt, daB sie in G stetige partielle Ableitungen zweiter Ordnung 
haben, inG+/J/'stetig sind und stetige partielle Ableitungen erster Ordnung haben. 
I" habe die Darstellung («(s), y(s)), OSs</ mit der Bogenlange s als Para- 
meter und stiickweise stetig differenzierbaren Funktionen x(s), y(s) und %?-++ y?>0 
in O0<s</. I" werde von s im mathematisch positiven Sinn durchlaufen. 
Durch ¢=(%,¥) und v=(—¥y,%) sind jedem Punkt von J’ ein Tangenten- 
vektor ¢ und ein innerer Normalvektor y zugeordnet (eventuell je zwei fiir dasselbe 
s als Grenzwerte von kleineren und gréBeren s her). 
Mit der Koeffizientenmatrix A von (2.1) la8t sich auf J” der sog. Konormalen- 
vektor o bilden: 


(2.2) C— Mmitt0* = Aye 


Da A positiv definit auf I" ist, gilt |o*|=-0 und (o,v)>0. Ableitungen von 
Funktionen in Richtung der Tangente oder der Konormale werden kurz durch 
‘die Indizes ¢ bzw. o gekennzeichnet. 

Bei der zweiten Randwertaufgabe werden oft auf dem Rande I’ die Werte 
vorgeschrieben fiir den Differentialausdruck 


(2.3) L*(u] = —|o*| u, = — (au,+ bu,)y — (bu, +¢4,) ae 
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a 


Satz 4. Sei p eine Losung der Differentialgleichung L{[p|=0 in G. I’ bestehe 
aus zwei Kurvenstiicken K,, K, (UU =K,+K,), so daB 
L*[p]S0 auf K, 2st, 
L*[p]} 20 auf K, ist. 
Dann gut: 
Wenn I glatt ist im Endpunkt P, von Ky, dann ist p, 20 in B,, 
Wenn I glatt ist im Endpunkt P, von K,, dann ist y,S0 im Fy. 
Dabei ist unter Kurvenstiick das stetige Bild eines Intervalles verstanden. 


Die Bezeichnung Endpunkt ist im Sinne der Orientierung von J’ durch den Bogen- 
parameter gemeint. 


Beweis. Die Hilfsfunktionen 
f(x,y) =ag,+y, und g(x,y) =—b9,—cyy 


sind eindeutig und stetig in G+J/’, haben stetige partielle Ableitungen erster 
Ordnung in G und aus L[y|=0 folgt #,=g, in G. Nach dem Hauptsatz tiber 
Kurvenintegrale ist die Funktion 


HC Sek asc 


mit festem Punkt P*€G und vorgegebenem Wert y(P*) in G+J’ eindeutig und 
hat stetige partielle Ableitungen erster Ordnung dort mit y,=—g, y,=/, unab- 
hangig vom gewahlten Integrationsweg, solange dieser in G+ J’ verlauft. w(x, y) 
gentigt in G einer elliptischen Differentialgleichung, denn y hat auch noch stetige 
partielle Ableitungen zweiter Ordnung in G und aus y,,—q,, in G folgt 


aD,+bD, piped bd 
aa 5 <= Vx — D . ae, 


4 perk lO e)ean (Ome +f Om) yan py = 0. 
Nach dem Randmaximumssatz nimmt also y seinen GroBtwert und seinen Kleinst- 
wert bezitiglich G+ J’ auf I’ an. 


Auf I’ gilt y,=L*[@], L*[y] =—D,, auf K, ist also y,=0, d.h. y(P)—y(B)=0, 
auf A, ist yO, d.h. p(P)—yp(B)=0, insgesamt p(P)=Sy(R) auf I. Ent- 
sprechend erhalt man y(P)<y(R) auf 7. yw nimmt also seinen Kleinstwert 
(GréBtwert) beziiglich G+J" in R (PR) an. Wenn I’ glatt ist in B, weist die 
Konormale ins Innere von G bei P.. Es ist dann also y, =0 in B, wegen D>0 
und |o*|>0 also 


1 Nee 


= —— ———— * [, 
a ee 
Entsprechend folgt 


pe 0) Puna, 


Yt D Wo = 0 nN 2? 
wenn /’ glatt ist dort. 


Fir die Durchfithrung einer Fehlerabschatzung bei der Randwertaufgabe 
[ee] = 7 (%) in NG= 
L*(u) = y(s) auf I 


Fehlerabschatzung bei elliptischen Differentialgleichungen 433 


sei v(x, y) eine Naherung mit L[v)=r(x, y) in G, L*(v]=y(s)+e(s) auf I. Sei 
p 


weiter q=> a; 4; (x, y) eine Lésung von L{g]=0 mit noch freien Parametern a,, 
j=1 p 
j—1,..., 0. Auf I ist L*[q|=> a,-y,(s)mit y,(s) = L*[q,]. 
j=1 
Nun sei & ein Punkt von J’, bei welchem J" glatt ist. Zu P gehore der 


Wert s, der Bogenlange von J’. Es miissen Zahlen a‘" so bestimmt werden, daB 


b 
fiir gQU= >) a. g, auf F gilt 
ABE 


2 
e(s) + Dial-y,(s) 20 fir s,Ss<sy, 
(2.4) f 
&(s) + dal. y;(s) SPOS Tar sf SSS, =F TS 
j=l 
wobei e(/-+s)=e(s) und y;(/+s)=y;(s) zu setzen ist und s, eine Zahl mit 
S;<S3<s,+/ ist. Dann ist nach Satz 4 


_ 


gi) <= vy, bel Bs 


Um eine méglichst gute Schranke zu erhalten, wird man zu den Ungleichungen 
(2.4) noch die Extremalbedingung 


(2.5) qr (R) = Min 


hinzunehmen bei der Bestimmung der a}. 


Entsprechend miissen Zahlen a\*) so bestimmt werden, daB gilt 


p 
2(S)2e ray (6) 0 Far SS FS sy) 
i qi (PR) = Max, 
e(s) + a)-y,(s) 20 fir s,Ss<s,+47, 
j=1 
denn damit erhalt man v,— u, < — gq!” bei B. 

Die Abschatzung mittels Satz 4 kann nur punktweise erfolgen. 

In den Anwendungen besonders wichtig sind zweite Randwertaufgaben bei 
unbeschranktem Grundgebiet. 

Sei B ein nichtbeschrankter Bereich der x, y-Ebene, der alle Punkte P mit 
v=r, fiir ein passend gewahltes 7) enthalt, wo r= /x?+ y? ist. Der Koordinaten- 
anfangspunkt liege auBerhalb B. Der Rand J’ von B bilde eine stiickweise glatte, 
geschlossene, doppelpunktfreie Kurve. Auf J” sei wieder die Bogenlange s einge- 
fiihrt, B liege zur Linken beim Durchlaufen von J’ mit wachsendem s. 


Satz 5. Sei p eine Lésung der Potentialgleichung Ap=O0 1m B. Es gebe positive 
Konstanten c und ¢, so daB in B gilt 


c 
vite © 


lgrad 9| < 


T bestehe aus zwei Kurvenstiicken K,, K,, so dab —y,S0 auf K,, — 9,20 auf 
Ky tst. 
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Wenn I glatt ist im Endpunkt P, von K,, dann ist py, 20 is As 
Wenn I glatt ist im Endpunkt P, von Ky, dann ist g,S0 im ae 


Der Beweis ist sehr einfach. Er wird hier nur angedeutet: 


y(P) ist beschrankt fiir r—>oo wegen | grad y| < 


C 
yite 


Die Transformation durch reziproke Radien 


x E 
f= 9 Me ees 
ye 0? ; 
mit. <%= 474)? 
wey: B53) 
if cae eee 


fiihrt den Bereich B in einen beschrankten Bereich B* der &, y-Ebene iiber. 
Es ist Ag*=0 in B* fiir p*(&, n)=9(x, y). Da p(x, y) beschrankt ist fiir 70, 
gilt dasselbe fiir m*(E,) bei g+0. Dann hat y* im Nullpunkt eine hebbare 
Singularitat, kann dort also so definiert werden, daB A p*=0 ist in dem einfach 
zusammenhangenden Bereich G*, der aus B* durch Hinzunahme des Null- 
punktes entsteht ((6], Sh 453). Nun 1la8t sich Satz 4 anwenden. Die Riick- 
transformation liefert die Aussage von Satz 5. 

SchlieBlich wird noch der Fall der raumlichen achsensymmetrischen Potential- 
gleichung bei nichtbeschranktem Bereich betrachtet. 

Sei B ein nichtbeschrankter, einfach zusammenhangender Bereich der oberen 
Halbebene y>0 in der x, y-Ebene. Der Rand I’ von B bestehe aus den beiden 
Abschnitten x<a und x=b der x-Achse und einem stiickweise glatten Kurven- 
stiick K, das ganz in der Halbebene y>0 liegt, ausgenommen Anfangspunkt P, 
und Endpunkt P,. K habe die Darstellung (x(s), y(s)), O<s</ mit der Bogen- 
lange s. K sei von s im mathematisch positiven Sinn durchlaufen (B zur Linken). 
t und y sind Tangenten- und innerer Normalvektor auf I. 


Satz 6. Sei m eine Lésung der Differentialgleichung = [(y . Jet(y y,),]=0 


mm B. Es gebe positive Konstanten c und e, so dap |gradg| S =< 
r=\x2+y2. Auf K sei (mit einer Zahl sy aus OSs, Sl) 


wr in B ist mit 


= Ys ON Ur, One GS Se 
= 0 0 Jie Sy SS Se 
P, sev der dem Parameterwert s, entsprechende Punkt von K. Wenn K glatt ist 
ber B,, dann gilt 
j= 0 1A. 


Fir g=const in B+ ist der Satz richtig. Fiir nicht konstantes y lauft 
der Beweis analog zum Beweis von Satz 4. Wie dort ist 


p(P) =y(P*) +f (gdx + f dy) 


mit f= y py, = — y@, nach Vorgabe von w(P*) eine eindeutige, stetige Funktion 
mit stetigen partiellen Ableitungen erster Ordnung in B+T° solange der Inte- 


’ 


grationsweg in B+ J" verlauft. In B hat p noch stetige zweite Ableitungen und 
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erfiillt die elliptische Differentialgleichung 
Pa PRE, Sie 
ey eadegs ay a 


y nimmt also seine untere Grenze beziiglich B in keinem inneren Punkt von B an. 
Auf K ist »(P) —y(R)=0 wegen 


0, Tuts; = 57 


wan 2S THEO S59 Ss; . 


Auf den Abschnitten x<a@ und «= der x-Achse gilt y(P) =y(B,) bzw. p(P) =y(B) 
wegen y=0. Also ist 
p(P)=y(B) fir Per. 


wie man leicht 


Fir 7(P) oo wird p(P)>y(B)=y(B) wegen |grad 7 se aie : 
P 


einsieht, wenn man in y(P) —p(B) =f (fdx+gdy) =f y (p, dx — —, dy) zunachst 
Py 
ein Stiick langst der x-Achse und dann iiber ein Rnoceanrack integriert. 


Es ist also y,=0 in A und wegen y,= yy, und y>0 folgt 
Cpe Oe Wier 


Bisher wurde immer gefordert, daB der Rand J" durch das Vorzeichen des Rand- 
ausdruckes L*[{q] in zwei Kurvenstiicke zerlegt wird. Eine Verallgemeinerung 
des Satzes 4 ist auch fiir den Fall méglich, da8 J’ durch das Vorzeichen von 
L*[q] in eine gerade Anzahl Kurvenstiicke mit jeweils wechselndem Vorzeichen 
zerlegt wird. Auch auf mehrfach zusammenhangende Gebiete kann der Satz 
erweitert werden, wenn man noch verlangt, daB $L*[p] ds=0 ist iiber jeden 
Randteil, der eine geschlossene Kurve bildet. Der praktischen Anwendbarkeit 
dieser Verallgemeinerungen steht allerdings die Tatsache entgegen, da auBer 
iiber L*[q] auch noch Kenntnisse iiber die Tangentialableitung y, an gewissen 
Randstellen erforderlich sind. 


3. Beispiele 
3.1. Randwertaufgabe mit nichtlinearer Randbedingung 
In Beispiel 4.3 von [9] wurde fiir die Lésung der Randwertaufgabe 


=Au—=1 inG: |*|<4, |gradu| =u auf I 


eine Naherung 
= y? 


Lara + Ay + a Vy 1 Ay Ve 
mit v,=Re[(x+7y)**], a9=0,848044, a= — 0,019562, a4,=0,000529 aufgestellt 
und die Fehlerabschatzung 

|%¢—v|S0,008 nG+Tr 


durchgefiihrt. Mit den Formeln von Satz1 laB8t sich damit eine Fehlerab- 
schatzung fiir die Ableitungen durchfiihren, die allerdings nur im Inneren von G 
endliche Schranken liefert. 
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Sei R(x, y) der Abstand des Punktes P= (x, y) von I’, d.h. 
R(x, y) =Min (1 —|«|,1—|y]). 


Dann gilt nach (1.2) fiir die Ableitung in einer beliebigen Richtung: 


epee 0,00637_ 
(x,y) — R(4, y) 


4 
| He — ¥2| S—- 0,005 - > 


also z.B. 
0,247 <4 — Ux, (z, 0) S 0,273 : 


3.2. Temperaturverteilung in einer Platte 
Eine Halfte einer quadratischen Platte werde erwarmt, die andere Halfte 
gekiihlt. Am Rande wird eine feste Temperatur eingehalten. Unter verein- 
fachenden Annahmen erhalt man fiir die Temperatur uw die Randwertaufgabe: 


—Au=6«x(i1—x) inG: Paige .eaNp) Ete © 
= = 

u=0 auf I’: Ms 1, |ylst, 
l=4, lalst. 


Der Ansatz 
sinj 7% -COSJ UY 
coshj x 


v=0,7% x 0,348 > a, Dts ag== 


p 
= 


Vo d 


erfiillt bereits die Differentialgleichung in G und die Randbedingungen auf den 
Randteilen x= +1, |y|<1. Eine Fehlerabschatzung fiir die Lésung selbst liefert 
der Randmaximumssatz. Es ist 


|u — v| = Max |u — | in ganz G+TJI". 
Fiir die Ableitung in beliebiger Richtung folgt damit nach (1.2): 


4 
Me Ye prey Mae eee 


wo R(x, y)=Min (1—|x 


,1—|y|) der Abstand vom Rand ist. 
Je nachdem, wieviele Ansatzfunktionen man beriicksichtigt, erhalt man: 
a, = — 0,23 53 


a, = — 0,2353 &,= +0,00735 
a, = —0,2365 a,= +0,0076 ag= —0,00108 


|ju—v|S 0,0080 0,0011 0,0003 


R(x, y) |ug—ve|S 0,0102 0,0014 0,0004 


also 0,6360S u, (0, 0) <0,6368. 


Gegen den Rand wachsen die Schranken iiber alle Grenzen. Bei den Rand- 
teilen «= +1, |y|S1 laBt sich Abhilfe schaffen, wenn man beriicksichtigt, daB 
der Fehler |~—v| auf diesen Randteilen verschwindet. Man kann etwa folgender- 
maBen verfahren: 
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Als Gebiet G* mit Rand J’** wird wegen des ungeraden Verhaltens aller 
Funktionen beziiglich des Nullpunktes der x-Achse das Teilgebiet 0<x<1, 


|y|<41 von G eingefiihrt. 
Nun werden zwei Niaherungen v® und v® bestimmt mit w—v™=0 auf I'*, 
v (1, 0) méglichst klein und u—v®<0 auf I’*, v®) (1, 0) méglichst ero8. Dann 


i : : (1) 2 
ist VVSusv® in G*+I"* oder mit Geren te 
(2) y(1) 
| ES eg 
2 


v0)(Q) —v®)(Q) 


Sei P=(x, y) mit festem y, |y|<1 und x=1—¢, 0<e<1—|y|. Nach (1.2) ist 
; 


Jus(P) —%9(P)| SF, Max |u(Q) — v(Q)] SA Max 


Nun ist, wenn man P = (1, y) einfiihrt, 
ax |/(Q)|S Max |f(Q)| 


M 
\O-Plse |O-B\s2e 
QcG* 


und, nach Mittelwertsatz und Schwartzscher Ungleichung 


|v (Q) — v (Q)| S|Q — B] - grad (& — o™)[p (0-2). 


Damit erhalt man 
4 1 
a aleg (7) SG 28 aa Mes, leretlel (0) (0) 
QEG* 


also fiir e—+0 die Abschatzung auf dem Rand: 
4 
|me(4, ) — vg (4, y)] S — [grad (v (1, y) — vo (1, y))]. 


Fiir den Ansatz mit zwei Parametern v =y,+ av, + a?) Op t= est 
|v — vp + 0,23655 v1 — 0,007525 vg| S 0,00225 v, — 0,000425 v2, 


also 
|g — Upe + 0,23655 ry ¢ — 0,007525 vo¢| S ee. |0,00225 vy — 0,000425 r,,| 


in den Randpunkten R= (1, y). Fir y=0 z.B. erhalt man 
0,7342 S — u, (1, 0) S$ 0,7373. 


3.3. Umstrémung eines Rotationsellipsords 
Das Geschwindigkeitspotential ® einer achsensymmetrischen raumlichen Stré- 
mung um ein Rotationsellipsoid geniigt unter vereinfachenden Annahmen der 


Randwertaufgabe 
po : : x 2 2 
—A®-() _in B: er ge fe + 22>1, q (20 
1m 71) =o 
r—>oo\ Ox 


a? eee 
Saamoan aig S41, 
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Dabei ist » die innere Normale auf I” und 9?=y?+ 2%, 7?=x?+ 02. Durch Wahl 
gecigneter Koeffizienten a; laBt sich eine Naherung 
p 
V=x4+)4,V; 
j=l 


mit V;=—, Vz= * (4x%— 602), ... aufstellen. 
We y@ 
Mit Hilfe des Randmaximumssatzes und der Transformation durch reziproke 
Radien hat L. Cottatz [2] eine Abschatzung angegeben von der Form 


|v—@|=<=* 
Vs 


mit einer Konstanten a, die vom Defekt in der Randbedingung abhangt. Damit 
laBt sich auch gleich wieder eine Fehlerabschatzung fiir die Ableitungen nach 
Nr.1 angeben. Man erhalt 
3a 
Vee Pl rae 

wobei d= Min[R(P), 7/2] ist mit dem Abstand R(P) von I’. 

Diese Schranken lassen sich nicht auf den Rand ausdehnen. Andererseits 
interessiert hier besonders die Tangentialgeschwindigkeit am Profilrand. 

Wenn man Zylinderkoordinaten x, 0, # einfiihrt und die Achsensymmetrie 
beriicksichtigt, erhalt man fiir ® eine Randwertaufgabe mit den in Satz 6 von 
Nr. 2 geforderten Eigenschaften. 


7 (le Pe + (ee) = 0 ih (BPA O-.0; + et>t, 


Pa : x = - (é@ e 
®,=0 aufK: 920, F+e=1, Jim(Z>—1)=0. 


Sei (®, — 1) =O (1/7*) und ©, =O (1/r°) fiir y+ 00, Wegen grad V;=0 (1/r°) ist auch 
b 

grad p=O (1/73) fir p=V —®=x+ > a,V;—@. Zur Durchfiihrung einer Ab- 
1 

schatzung fir die Tangentialableitung im Punkt PR. mit «=0, e=1 werden zwei 

Naherungen (i) (i) (i : 

VO = 4 a Ve av. (== 452) 

aufgestellt, so daB fir po =V" —®@ auf I gilt 


—p20, —gpS<0 fir x<0, 
— Gr 0, =O 20 time O: 

Dann ist 

Gla late 


VO<G,< VO in P. 


Man erhalt 2.B. VY=x%+0,35V,+0,017V,, V®=x+1,08V,+0,088V, und da- 
mit die Schranken 
1,248 <@,<1,552 im Punkt BR (x=0, o=1). 


Diese Arbeit ist ein Teilauszug aus der Dissertation des Verfassers. Ein anderer 
Teilauszug wurde in [9] gegeben. 
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ni Bemerkung 


: zur Arbeit H. EHRMANN, Iterationsverfahren mit verdnderlichen Operatoren, 
Vol. 4, 45 —64 (1959). 


a Herr J. SCHRODER aus Hamburg machte mich freundlicherweise darauf auf- 

merksam, daB das Zusatzaxiom (2.1g) auf S. 47 eine haufig nicht erfiillte Forde- 
rung darstellt. 

" Unter der Voraussetzung B 1a geniigt es jedoch, stets gleichmaBige Konver- 

_-genz, zu fordern, also (2.14g) wie folgt abzuschwachen: 

_ Aus lim e,=0 und 0Se, fiir »=—1, 2,... folgt die Existenz eines (festen) 
Elementes t derart, daB fiir jede positive Zahl k ein ,(k) existiert mit 0, <t/k 
fiir alle n=n,(). f 


y Der Beweis bleibt dabei unveradndert. 
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